Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,918,022 A→G 5.9% intergenic (+209/+23) barA → / ← gudD hybrid sensory histidine kinase, in two‑component regulatory system with UvrY/D‑glucarate dehydratase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,918,0220AG 5.9% 220.1 / 3.7 68intergenic (+209/+23)barA/gudDhybrid sensory histidine kinase, in two‑component regulatory system with UvrY/D‑glucarate dehydratase 1
Reads supporting (aligned to +/- strand):  ref base A (38/26);  new base G (2/2);  total (40/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.91e-01

ACTGTTGGACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACTTGATCCATATCGATTTCTACACCCAGCCCCGGTTTTTCTGGCACCTGTACCAGCCCGCCTTTGATCTCAAACGGTTCTTTGGTCAGGCGCTGATTGCCTTCCTGCCAAATCCAGTGCGTATCAATAGC  >  NC_000913/2917759‑2918296
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                   
aCTGTTGGACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggtttt                                                                                                                                                                                                                                                                          >  2:224284/1‑275 (MQ=255)
                                                                ttgtcggtggaaaTCGGCTTTTTTTGACGCTGGTTTAAGATTGGTTACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGa                                                                                                                                                                                                          <  4:293274/263‑1 (MQ=255)
                                                                                                                          tgtTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAGATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGTCGGTTTTGGTTTTCTGTCTTAACGCCCCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCCTTCCCATCCCATCGTCCCGCCCGCCCATCCCGTGTTTCTTATACAGTTTATg                                                                                                                                               >  2:93117/1‑275 (MQ=255)
                                                                                                                             ttATGCCGAATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCAGAACAGGCTATTTTGTGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttt                                                                                                                                                                                                                                                                           <  6:388427/149‑1 (MQ=255)
                                                                                                                               aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttg                                                                                                                                                                                                                                                                         >  7:194670/1‑149 (MQ=255)
                                                                                                                                     gataaggcgTTAACGCCGCATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTc                                                                                                                                                                                                                                                                   >  5:90477/6‑149 (MQ=255)
                                                                                                                                     gATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTc                                                                                                                                                                                                                                                                   <  5:54402/149‑1 (MQ=255)
                                                                                                                                      ataaggcgTTAACGCCGCATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCt                                                                                                                                                                                                                                                                  <  6:90477/145‑1 (MQ=255)
                                                                                                                                          ggcgTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTACTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCt                                                                                                                                                                                                                                                              >  5:340947/1‑149 (MQ=255)
                                                                                                                                               tAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCt                                                                                                                                                            >  2:236670/1‑241 (MQ=255)
                                                                                                                                                   gCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATa                                                                                                                                                        <  4:236666/241‑1 (MQ=255)
                                                                                                                                                     cTTATCCGGCCTAGAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCAt                                                                                                                                                                                                                                                   <  8:301776/149‑1 (MQ=255)
                                                                                                                                                             gCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGACTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGgc                                                                                                                                                                                                                                           <  8:380885/149‑1 (MQ=255)
                                                                                                                                                              ccTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCCCATGGGCTTTCATCACTTTATCCATATCGATTTCTACAcc                                                                                                           >  2:492612/1‑275 (MQ=255)
                                                                                                                                                                          aTTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTGCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGACCAGCCAGGAATCAGATAc                                                                                                                                                                                                       >  2:515898/1‑171 (MQ=255)
                                                                                                                                                                          aTTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAAt                                                                                                                                                                                                              >  2:343932/1‑164 (MQ=255)
                                                                                                                                                                              aaTTCAATGAATANAAAANTGATGTAGACCTGATACATAGCACAACATGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCGCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGAACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCa                                                                                                                                                                                                   <  4:515893/171‑1 (MQ=255)
                                                                                                                                                                               aTTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGaa                                                                                                                                                                                                                             >  6:199510/1‑144 (MQ=255)
                                                                                                                                                                               aTTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGaa                                                                                                                                                                                                                             <  5:199510/144‑1 (MQ=255)
                                                                                                                                                                                     tGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAAc                                                                                                                                                                                                                                                          <  8:22545/110‑1 (MQ=255)
                                                                                                                                                                                     tGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAAc                                                                                                                                                                                                                                                          >  7:22545/1‑110 (MQ=255)
                                                                                                                                                                                     tGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTccagccag                                                                                                                                                                                                                  >  6:106466/1‑149 (MQ=255)
                                                                                                                                                                                            aaaaaTGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAg                                                                                                                                                                                                           >  7:429570/1‑149 (MQ=255)
                                                                                                                                                                                                   atgTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCCTCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACTTGGTCCATATCGATTTCTACACCCAGCCCCGTTTTTTCTGGCCACTTTACCAGCACCCCt                                                                      >  4:347279/1‑275 (MQ=255)
                                                                                                                                                                                                         aCCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATTCTGCATTccc                                                                                                                                                                                              >  7:10590/1‑149 (MQ=255)
                                                                                                                                                                                                         aCCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGc                                                                                                                                                                      >  2:120120/1‑173 (MQ=255)
                                                                                                                                                                                                            tGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTACAGCCAGGAATCAGATACTGCATTCCcatc                                                                                                                                                                                           >  7:23030/1‑149 (MQ=255)
                                                                                                                                                                                                             gATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGaa                                                                                                                                                                                                               >  1:71443/1‑128 (MQ=255)
                                                                                                                                                                                                             gATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCgt                                                                                                                                                                  <  4:120119/173‑1 (MQ=255)
                                                                                                                                                                                                              aTACATATCACAACAGACTATTTTGCGTTTAGCATCTGTCTCCAACCGGCTCCAGATAGAGCCGGTTTTGGTTTGCTGTCTTAACTCACCATGCACGGGC‑GCTTGTTATCGCAAGTCCAGCCAGCAATCCAATACTGCATTCCCATCCAATCGTCACTCCCGCCCCGCCCGGTTTTCTGATAGACCTCATGTGCTTCAATACCTTAATCCTTATCGTTTTCTACACCTACCCCGGGTTTTTACGGCACTGGTACCATCCCCCCCTTgccaccaa                                                            >  4:42479/1‑267 (MQ=255)
                                                                                                                                                                                                                aCATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttggtttt                                                                                                                                                                                                                                                                    >  6:424938/1‑73 (MQ=255)
                                                                                                                                                                                                                aCATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttggtttt                                                                                                                                                                                                                                                                    <  5:424938/73‑1 (MQ=255)
                                                                                                                                                                                                                 cATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAg                                                                                                                                                                                                           <  3:71443/128‑1 (MQ=255)
                                                                                                                                                                                                                    aGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGc                                                                                                                                                                      >  2:137499/1‑162 (MQ=255)
                                                                                                                                                                                                                     gCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACTTGATCCATATCGATTTCTACACCCGGCCCCGGTTTTTCTGGCACCTGTACCAGCCCGCCTTTGATTCCAAACGGtttt                                                    >  2:400630/1‑273 (MQ=255)
                                                                                                                                                                                                                     gCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTAGTTTTCTGGCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAAGGTCCAGCCCGGAATCAGATACTGCATTCCCATGGCAACGTCACGACCGCCAAGCCCGTGTTTCTGCTACCGCTAATGGCCTTTCATTACTTTCTCCATTTCGATTTCCCCAACCAACCCCGGTTTTCGTGGCCACTGCAACAGCCCCCCTTTGACCTCAACCGGTTCt                                                    >  4:399589/1‑275 (MQ=255)
                                                                                                                                                                                                                        caacaGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCgt                                                                                                                                                                  <  4:137498/162‑1 (MQ=255)
                                                                                                                                                                                                                        caacaGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTc                                                                                                                                                             >  4:413247/1‑167 (MQ=255)
                                                                                                                                                                                                                         aacaGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTa                                                                                                                                                                                                                                  >  6:210602/1‑97 (MQ=255)
                                                                                                                                                                                                                         aacaGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTa                                                                                                                                                                                                                                  <  5:210602/97‑1 (MQ=255)
                                                                                                                                                                                                                         aacaGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACATCTCATGGGCTTTCATCACTTCATCCATATCGATTGCTACACCCAGCCCCGGTTTTTCTGGCATCTGTCCCCGCCCCACTTTTATCCCAAACGGTTATTTgg                                                >  1:75256/1‑275 (MQ=255)
                                                                                                                                                                                                                          acaGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCAcgcgc                                                                                                                                                                             <  8:194670/149‑1 (MQ=255)
                                                                                                                                                                                                                            aGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGAt                                                                                                                                                         <  2:413252/167‑1 (MQ=255)
                                                                                                                                                                                                                              gCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttgg                                                                                                                                                                                                                                                                        >  2:468818/1‑55 (MQ=255)
                                                                                                                                                                                                                                  ttttGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttggtttt                                                                                                                                                                                                                                                                    <  4:468813/55‑1 (MQ=255)
                                                                                                                                                                                                                                   tttGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGccc                                                                                                                                                                    >  7:255775/1‑149 (MQ=255)
                                                                                                                                                                                                                                      gCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCgtg                                                                                                                                                                 <  8:429570/149‑1 (MQ=255)
                                                                                                                                                                                                                                               aTCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCCTTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATa                                                                                                                                                        >  7:420823/1‑149 (MQ=255)
                                                                                                                                                                                                                                               aTCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATa                                                                                                                                                        >  6:89257/1‑149 (MQ=255)
                                                                                                                                                                                                                                                  aGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAAc                                                                                                                                                                                                                            >  2:361820/1‑78 (MQ=255)
                                                                                                                                                                                                                                                   gTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCgtg                                                                                                                                                                 >  1:116324/1‑136 (MQ=255)
                                                                                                                                                                                                                                                      tcAAACCGGCTCCAGAGAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTcc                                                                                                                                                                                                                        <  4:361815/78‑1 (MQ=255)
                                                                                                                                                                                                                                                      caaaaCCGGCTCTATCTGGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCcatcgca                                                                                                                                                                                        >  8:145998/3‑110 (MQ=39)
                                                                                                                                                                                                                                                      caaaaCCGGCTCTATCTGGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCcatcgca                                                                                                                                                                                        <  7:145998/108‑1 (MQ=39)
                                                                                                                                                                                                                                                       cAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATTCTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCAt                                                                                                                                                <  8:10590/149‑1 (MQ=255)
                                                                                                                                                                                                                                                       cAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTc                                                                                                                                                             <  3:116324/136‑1 (MQ=255)
                                                                                                                                                                                                                                                       cAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGa                                                                                                                                                          >  5:94855/1‑139 (MQ=255)
                                                                                                                                                                                                                                                       cAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGa                                                                                                                                                          <  6:94855/139‑1 (MQ=255)
                                                                                                                                                                                                                                                        aaaCCGGCTCTATCTGGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTc                                                                                                                                                             >  8:313483/1‑135 (MQ=255)
                                                                                                                                                                                                                                                        aaaCCGGCTCTATCTGGAGCCGGTTATGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTc                                                                                                                                                             <  7:313483/135‑1 (MQ=255)
                                                                                                                                                                                                                                                        aaaCCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACTTGATCCATATCGATTTCTACACCCAGCCCCGGTTTTTCTGGCACCTGTACCAGCCCGCCTTTGATCCCAAACGGTTCTTTGGTCAGGCGCTGATTGCCTTTCTGCCAAATCCa                 >  2:466826/1‑275 (MQ=255)
                                                                                                                                                                                                                                                              gCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCgtg                                                                                                                                                                 >  2:298031/1‑125 (MQ=255)
                                                                                                                                                                                                                                                              gCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACTTGATCCATATCGATTTCTACACCCAGCCCCGGTTTTTCTGGCACCTGTACCAGCCCGCCTTTGATCTCAAACGGTTCTTTGGTCAGGCGCTGATTTCCTTCCTGCCAAATCCAGTGcct           >  2:269083/1‑273 (MQ=255)
                                                                                                                                                                                                                                                               cTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTtc                                                                                                                                        <  5:106466/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                tCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTtca                                                                                                                                       >  5:220426/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  cAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTc                                                                                                                                                             <  4:298027/125‑1 (MQ=255)
                                                                                                                                                                                                                                                                  cAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACGGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACTTGATCCATATCGATTTCTACACCCAGCCCCGGTTTTTCTGGCACCTGTACCAGCCCGCCTTTGATCTCAAACg                                                         <  2:157974/225‑1 (MQ=255)
                                                                                                                                                                                                                                                                   aGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTtcatca                                                                                                                                    >  6:249322/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                       agagTCGGTTTATGTTTTCTGTCTTAACTACCCATACACGGGCGGCTTG‑TCTCGAACGGCGAGCCAGGAACCAGAGTCTGAATTCCCATCGCATCGTCACGAGCGCCAAGCCCGTGTTTCTTATACAACTCATGGGCTTTCATCACTTGTTCCATATCGCTTGCTACACCCAGCCCCGGTTTTTCTGGCACCTGTACGAGCCCGCCTTTGATCTCAAACGGTTCTTTGGTCAGGCGCTGATTGCCTTCCTGCCAAAACCAGTGCGTATCAATAGc  <  4:400625/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                   
ACTGTTGGACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGC‑GCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACTTGATCCATATCGATTTCTACACCCAGCCCCGGTTTTTCTGGCACCTGTACCAGCCCGCCTTTGATCTCAAACGGTTCTTTGGTCAGGCGCTGATTGCCTTCCTGCCAAATCCAGTGCGTATCAATAGC  >  NC_000913/2917759‑2918296

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑