Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,297,194 A→T 5.0% intergenic (‑69/+25) yjcO ← / ← fdhF Sel1 repeat‑containing protein YjcO/formate dehydrogenase H

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,297,1940AT 5.0% 921.9 / 20.8 273intergenic (‑69/+25)yjcO/fdhFSel1 repeat‑containing protein YjcO/formate dehydrogenase H
Reads supporting (aligned to +/- strand):  ref base A (135/124);  new base T (6/8);  total (141/132)
Fisher's exact test for biased strand distribution p-value = 5.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 278 total aligned reads displayed.

AAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCG  >  NC_000913/4297046‑4297343
                                                                                                                                                    |                                                                                                                                                     
     cTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAaggag                                                                                                                                                 >  1:1760780/1‑150 (MQ=255)
       ggCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTC‑AAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCt                                                                                                                                              >  2:624481/1‑150 (MQ=255)
         cTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTg                                                                                                                                             >  1:1982844/1‑150 (MQ=255)
           gCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCTTTCGGAGGCTGtt                                                                                                                                           <  2:951760/150‑2 (MQ=255)
           gCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTa                                                                                                                                           >  1:85930/1‑150 (MQ=255)
            ccAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCTTTCGGAGGCTGttt                                                                                                                                          <  1:105970/150‑3 (MQ=255)
            ccAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAg                                                                                                                                          <  2:32472/150‑1 (MQ=255)
            ccAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAg                                                                                                                                          >  1:388343/1‑150 (MQ=255)
              aGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCTTTCGGAGGCTGttttt                                                                                                                                        <  1:326600/150‑5 (MQ=255)
               ggCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGACTCCGAAAGGAGGCTGTAGaaa                                                                                                                                       >  2:237934/1‑150 (MQ=255)
                gCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAg                                                                                                                                      <  2:85930/150‑1 (MQ=255)
                 cTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCTTTCGGAGGCTGttttttta                                                                                                                                     <  1:1001797/150‑8 (MQ=255)
                 cTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAgg                                                                                                                                     <  1:1152593/150‑1 (MQ=255)
                 cTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAgg                                                                                                                                     <  2:106102/150‑1 (MQ=255)
                  tCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGa                                                                                                                                    >  2:1913920/1‑150 (MQ=255)
                   cGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGAc                                                                                                                                   >  2:1619683/1‑150 (MQ=255)
                      aGTCGTTGGCGTGTGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGt                                                                                                                                <  2:712620/150‑1 (MQ=255)
                       gTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTa                                                                                                                               >  1:243728/1‑150 (MQ=255)
                         cGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTAtt                                                                                                                             >  2:1953632/1‑150 (MQ=255)
                         cGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTAtt                                                                                                                             >  1:1884119/1‑150 (MQ=255)
                            tGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACg                                                                                                                          <  1:293124/150‑1 (MQ=255)
                             ggCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGc                                                                                                                         >  2:1410048/1‑150 (MQ=255)
                             ggCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGc                                                                                                                         >  2:801020/1‑150 (MQ=255)
                              gCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGcc                                                                                                                        <  1:33463/150‑1 (MQ=255)
                              gCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGcc                                                                                                                        <  1:492764/150‑1 (MQ=255)
                               cGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCa                                                                                                                       >  2:198516/1‑150 (MQ=255)
                                   ggCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGc                                                                                                                   >  2:1073897/1‑150 (MQ=255)
                                    gCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGcc                                                                                                                  >  1:196781/1‑150 (MQ=255)
                                         gAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCCCCTTTCGGAGGCTGt                                                                                                                                            <  2:551498/119‑1 (MQ=255)
                                           aTGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTcgc                                                                                                           <  2:1520211/150‑1 (MQ=255)
                                             gTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCgcgc                                                                                                         >  2:447288/1‑150 (MQ=255)
                                             gTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCgcgc                                                                                                         >  1:461035/1‑150 (MQ=255)
                                              tCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTa                                                                                                                               >  1:903186/1‑127 (MQ=255)
                                              tCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCa                                                                                                        >  2:1586737/1‑150 (MQ=255)
                                               cAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACg                                                                                                                                  >  1:1480537/1‑123 (MQ=255)
                                               cAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAg                                                                                                       >  1:1373220/1‑150 (MQ=255)
                                               cAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAg                                                                                                       >  2:514273/1‑150 (MQ=255)
                                               cAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAg                                                                                                       >  1:183457/1‑150 (MQ=255)
                                                 aaaaCAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGt                                                                                                                                <  2:913898/123‑1 (MQ=255)
                                                 aaaaCAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGt                                                                                                                                >  1:913898/1‑123 (MQ=255)
                                                  aaaCAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCg                                                                                                    <  1:823528/150‑1 (MQ=255)
                                                      aacaTCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGtt                                                                                                <  2:393359/150‑1 (MQ=255)
                                                        caTCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGtttt                                                                                              <  1:1619683/150‑1 (MQ=255)
                                                          tCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCa                                                                                            <  1:750285/150‑1 (MQ=255)
                                                          tCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCa                                                                                            <  1:510222/150‑1 (MQ=255)
                                                            aaTGCGATAATTTTTTTCATTTTGATAATACGATTGGGTTGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCGGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAAc                                                                                          <  2:1476322/150‑1 (MQ=255)
                                                               gCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAAc                                                                                          >  2:181590/1‑147 (MQ=255)
                                                                cGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACttgt                                                                                      >  2:1459688/1‑150 (MQ=255)
                                                                cGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACttgt                                                                                      <  2:432290/150‑1 (MQ=255)
                                                                   tAATTTTTTTCATTTTGATAATACGATT‑GGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAg                                                                                                                      >  1:1781206/1‑114 (MQ=255)
                                                                     aTTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTAc                                                                                 >  1:888981/1‑150 (MQ=255)
                                                                     aTTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGATGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTTTAc                                                                                 >  1:1103926/1‑150 (MQ=255)
                                                                            tCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCaa                                                                          >  2:360959/1‑150 (MQ=255)
                                                                             cATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAAt                                                                         <  2:1760780/150‑1 (MQ=255)
                                                                              aTTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATc                                                                        <  2:820270/150‑1 (MQ=255)
                                                                              aTTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATc                                                                        >  2:1010476/1‑150 (MQ=255)
                                                                                 ttGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACg                                                                     >  1:700226/1‑150 (MQ=255)
                                                                                  tGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACttgttg                                                                                    <  2:976513/134‑1 (MQ=255)
                                                                                  tGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACttgttg                                                                                    >  1:976513/1‑134 (MQ=255)
                                                                                  tGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGt                                                                    >  2:1857237/1‑150 (MQ=255)
                                                                                  tGATAATACGATTGGGTAGGGAGT‑ACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGcc                                                                                                                        >  1:1538847/1‑97 (MQ=255)
                                                                                   gATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTa                                                                   <  1:1913920/150‑1 (MQ=255)
                                                                                   aataataCGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAg                                                                                                                                          >  1:144392/2‑79 (MQ=255)
                                                                                    ataataCGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGt                                                                    <  2:255593/148‑1 (MQ=255)
                                                                                    ataataCGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTAc                                                                  <  1:1773324/150‑1 (MQ=255)
                                                                                    ataataCGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTAc                                                                  <  2:1667299/150‑1 (MQ=255)
                                                                                       ataCGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGCCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGc                                                               >  1:1409052/1‑150 (MQ=255)
                                                                                        taCGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCt                                                              >  1:1560693/1‑150 (MQ=255)
                                                                                          cGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTAc                                                                                 >  2:894769/1‑129 (MQ=255)
                                                                                              tGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTcggcgg                                                        <  1:1023459/150‑1 (MQ=255)
                                                                                              tGGGTAGGGAGTAACCAGGATAAAGACTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTcggcgg                                                        <  2:1741303/150‑1 (MQ=255)
                                                                                                   aGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCATCCGAAAGGAGGCTGTAGAAAGGACGGTATTACCCCAGTCCCCCTACGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGgcgcgc                                                   >  2:1386212/1‑150 (MQ=255)
                                                                                                    gggAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCt                                                  >  2:1798765/1‑150 (MQ=255)
                                                                                                      gAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTga                                                >  1:1689063/1‑150 (MQ=255)
                                                                                                       aGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTgat                                               >  1:557227/1‑150 (MQ=255)
                                                                                                           aCCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTgatcggc                                           <  2:1689063/150‑1 (MQ=255)
                                                                                                            ccAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCCgatcgccg                                          >  1:1607415/1‑150 (MQ=255)
                                                                                                                tataAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTgatcggcgatcg                                      <  2:1560693/150‑1 (MQ=255)
                                                                                                                tataAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTgatcggcgatcg                                      <  1:1238029/150‑1 (MQ=255)
                                                                                                                tataAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTgatcggcgatcg                                      >  1:181005/1‑150 (MQ=255)
                                                                                                                 ataAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTgatcggcgatcgg                                     >  1:596830/1‑150 (MQ=255)
                                                                                                                        gTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGAcg                              <  1:1598979/150‑1 (MQ=255)
                                                                                                                           cGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGa                           >  2:1431151/1‑150 (MQ=255)
                                                                                                                            caaTGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGTCGAGCCGATCGCCGATCAGCGcgccgccga                          <  1:1705365/149‑4 (MQ=255)
                                                                                                                                ggATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGgc                      >  1:316009/1‑150 (MQ=255)
                                                                                                                                 gATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGgcg                     >  2:136377/1‑150 (MQ=255)
                                                                                                                                  aTAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGgcgc                    <  1:419132/150‑1 (MQ=255)
                                                                                                                                   tAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCa                   >  1:581040/1‑150 (MQ=255)
                                                                                                                                   tAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCa                   >  2:1273409/1‑150 (MQ=255)
                                                                                                                                     aaaaaaCAGCCTCCGAAAGGAGGGTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTGAATCACGTACTGCTCGGCGGCGCGCTGATCGGCTATCGGCTCGACGCGAACGGCGCAGt                 <  1:1369273/150‑1 (MQ=255)
                                                                                                                                        aaaCAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTAtt              >  2:1778682/1‑150 (MQ=255)
                                                                                                                                        aaaCAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTAtt              >  1:669620/1‑150 (MQ=255)
                                                                                                                                             gCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATTGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTAc         <  1:801020/150‑1 (MQ=255)
                                                                                                                                              ccTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACt        >  2:496251/1‑150 (MQ=255)
                                                                                                                                               cTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTg                                                 <  2:1495735/108‑1 (MQ=255)
                                                                                                                                                tCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTcc      >  1:1666912/1‑150 (MQ=255)
                                                                                                                                                 ccGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACttgttg                                                                                    >  2:1128111/1‑71 (MQ=255)
                                                                                                                                                  cGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCgg    >  1:787455/1‑150 (MQ=255)
                                                                                                                                                    aaaGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCg  <  1:1717210/150‑1 (MQ=255)
                                                                                                                                                    aaaGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCg  <  2:1251980/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCG  >  NC_000913/4297046‑4297343

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑