Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,212,477 G→A 5.9% intergenic (‑173/+63) ymgK ← / ← ymgL protein YmgK/protein YmgL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,212,4770GA 5.9% 482.1 / 11.1 151intergenic (‑173/+63)ymgK/ymgLprotein YmgK/protein YmgL
Reads supporting (aligned to +/- strand):  ref base G (66/76);  new base A (3/6);  total (69/82)
Fisher's exact test for biased strand distribution p-value = 5.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.02e-01

Only 100 of 153 total aligned reads displayed.

AATAATGAATATAAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGA  >  NC_000913/1212330‑1212626
                                                                                                                                                   |                                                                                                                                                     
aataatGAATATAAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCa                                                                                                                                                     <  2:730697/150‑1 (MQ=255)
aataatGAATATAAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCa                                                                                                                                                     <  2:685049/150‑1 (MQ=255)
  taatGAATATAAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACg                                                                                                                                                   >  1:197812/1‑150 (MQ=255)
     tGAATATAAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGcc                                                                                                                                                <  2:1310251/150‑1 (MQ=255)
     tGAATATAAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGcc                                                                                                                                                <  2:521796/150‑1 (MQ=255)
       aaTATAAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCgt                                                                                                                                              >  2:1259520/1‑150 (MQ=255)
         tataAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATTATTTAAAATGTTCCTGGGTAAAATTATTTTTACCGGATTTTTAATAAATAAAAATAAATGTGCGTTAAGTATTTAAACGGCCGTGCCACACGTCCgttt                                                                                                                                            <  2:669798/150‑2 (MQ=255)
           taAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGCCGTGCCACACGTCCgtttca                                                                                                                                          <  1:597445/150‑4 (MQ=255)
           cttgCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCgttgtt                                                                                                                                          <  2:1342032/147‑1 (MQ=255)
            aaGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGttt                                                                                                                                         >  1:1365869/1‑150 (MQ=255)
            aaGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGttt                                                                                                                                         <  2:1018519/150‑1 (MQ=255)
             aGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGtttt                                                                                                                                        <  2:589899/150‑1 (MQ=255)
             aGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGtttt                                                                                                                                        <  1:665681/150‑1 (MQ=255)
                atgacatgaGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTtat                                                                                                                                     <  1:1225558/150‑1 (MQ=255)
                 tgacatgaGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTtata                                                                                                                                    <  1:990349/150‑1 (MQ=255)
                  gacatgaGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATaa                                                                                                                                   >  2:1068208/1‑150 (MQ=255)
                   acatgaGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATaaa                                                                                                                                  <  2:813177/150‑1 (MQ=255)
                    catgaGAATAAGGTTTCAATTTTTGAGTTATATAGTAATGATTTAAGCTGTTCCTGGCTAAAATAAATATAACCGGATTTTGACTAAGCGAAATTAGATGTGCGTTAAGTATTGAAACGGACGTGTGGGACGGCCGTTGTTTTTATAAat                                                                                                                                 <  2:1071304/150‑1 (MQ=255)
                     atgaGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTGACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAata                                                                                                                                <  1:985044/150‑1 (MQ=255)
                        agaATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCATAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGt                                                                                                                             <  1:733298/150‑1 (MQ=255)
                         gaATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGtt                                                                                                                            >  2:258029/1‑150 (MQ=255)
                           aTAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTaa                                                                                                                          <  1:948744/150‑1 (MQ=255)
                            tAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAAc                                                                                                                         >  2:872921/1‑150 (MQ=255)
                              aGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCg                                                                                                                       <  2:19133/150‑1 (MQ=255)
                               ggTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTTTAAATATGTTAACCGt                                                                                                                      <  2:1081184/150‑1 (MQ=255)
                                 tttCAATTTTTGAGTTATATAGGAATGATTTAACCTGTGCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTta                                                                                                                    <  1:1007339/150‑1 (MQ=255)
                                    cAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATaa                                                                                                                 <  1:48803/150‑1 (MQ=255)
                                    cAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATaa                                                                                                                 <  1:886288/150‑1 (MQ=255)
                                     aaTTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATaaa                                                                                                                >  1:859524/1‑150 (MQ=255)
                                           tGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACg                                                                                                          <  2:708311/150‑1 (MQ=255)
                                           tGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACg                                                                                                          <  2:39954/150‑1 (MQ=255)
                                               ttATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATTTTATAACggtta                                                                                                      <  1:328309/150‑5 (MQ=255)
                                                tatataGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCa                                                                                                     <  1:413744/150‑1 (MQ=255)
                                                  tataGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATTTTATAACggttaaca                                                                                                   <  1:511694/150‑8 (MQ=255)
                                                      ggAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTc                                                                                               <  2:1258880/150‑1 (MQ=255)
                                                      ggAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTc                                                                                               >  1:627109/1‑150 (MQ=255)
                                                      ggAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATAAATGTGCGTTAAGTATTGAAACGGACTTGTGGCAAGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGta                                                                                               >  1:558093/1‑149 (MQ=255)
                                                       gAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGAGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCa                                                                                              <  1:918847/150‑1 (MQ=255)
                                                           gATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGGTGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTg                                                                                          >  2:1338303/1‑150 (MQ=255)
                                                           gATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTg                                                                                          >  1:14694/1‑150 (MQ=255)
                                                              ttAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACTTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATc                                                                                       <  2:825724/150‑1 (MQ=255)
                                                              ttAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATc                                                                                       <  1:1258214/150‑1 (MQ=255)
                                                              ttAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATc                                                                                       >  2:647780/1‑150 (MQ=255)
                                                               tAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATca                                                                                      <  1:1096682/150‑1 (MQ=255)
                                                               tAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATca                                                                                      >  2:1300256/1‑150 (MQ=255)
                                                                        ccTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCaaa                                                                             >  1:55161/1‑150 (MQ=255)
                                                                               aaaaTACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAg                                                                      <  1:470130/150‑1 (MQ=255)
                                                                               aaaaTACATATAACCGGATGATGACTAAACCAAAATAAATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAg                                                                      <  1:382216/150‑1 (MQ=255)
                                                                                     cATATAACCGGATGATGACTAAACAAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGATTTTATAAATATGTTAACAGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTCATa                                                                >  1:444883/1‑150 (MQ=255)
                                                                                      atatAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAg                                                               >  1:1420226/1‑150 (MQ=255)
                                                                                       tataACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGt                                                              >  2:37227/1‑150 (MQ=255)
                                                                                         taACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGt                                                                                                         >  2:224953/1‑105 (MQ=255)
                                                                                          aaCCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTAtt                                                           <  2:1113748/150‑1 (MQ=255)
                                                                                          aaCCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTAtt                                                           <  1:258029/150‑1 (MQ=255)
                                                                                           aCCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAATAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTa                                                          >  2:1239815/1‑150 (MQ=255)
                                                                                           aCCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAata                                                                           <  2:319668/133‑1 (MQ=255)
                                                                                            ccGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAg                                                         >  1:946381/1‑150 (MQ=255)
                                                                                             cGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTAGCAAAAGTCAAGTGCTCCCAGTTCACATATCAAGTGGATAGTATTAgg                                                        >  1:1038032/1‑149 (MQ=255)
                                                                                                atgatgACTAAACCAAAATACATGTTCTTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTg                                                     <  1:1388586/150‑1 (MQ=255)
                                                                                                atgatgACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTg                                                     <  1:608600/150‑1 (MQ=255)
                                                                                                 tgatgaCTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTgg                                                    <  1:872921/150‑1 (MQ=255)
                                                                                                   atgaCTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTg                                                  <  2:875030/150‑1 (MQ=255)
                                                                                                    tgaCTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGa                                                 >  1:923035/1‑150 (MQ=255)
                                                                                                     gaCTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGat                                                >  2:1169436/1‑150 (MQ=255)
                                                                                                     gaCTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGat                                                >  1:281141/1‑150 (MQ=255)
                                                                                                      aCTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGatt                                               >  1:392526/1‑150 (MQ=255)
                                                                                                         aaaCCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGattatt                                            >  2:614947/1‑150 (MQ=255)
                                                                                                           aCCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGAttattta                                          >  1:240873/1‑150 (MQ=255)
                                                                                                            ccAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGAttatttat                                         <  2:627109/150‑1 (MQ=255)
                                                                                                            ccAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGAttatttat                                         >  2:677222/1‑150 (MQ=255)
                                                                                                                 aTACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTg                                    <  2:285906/150‑1 (MQ=255)
                                                                                                                  tACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTta                                                                                                                                      <  1:1374909/51‑1 (MQ=255)
                                                                                                                  tACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTta                                                                                                                                      >  2:1374909/1‑51 (MQ=255)
                                                                                                                  tACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATgg                                       >  2:286361/1‑146 (MQ=255)
                                                                                                                  tACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATgg                                       <  1:286361/146‑1 (MQ=255)
                                                                                                                     aTGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACaa                                >  1:889239/1‑150 (MQ=255)
                                                                                                                         gCGTTAAGTATTGAAACGGACGTGTGACACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATaaa                            >  2:1370076/1‑150 (MQ=255)
                                                                                                                            ttAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAg                         <  2:474537/150‑1 (MQ=255)
                                                                                                                                gTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGt                                                                                                                      <  1:663665/53‑1 (MQ=255)
                                                                                                                                  aaaacaaCGGCCGTGCCACACGTCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTa                   >  1:800519/6‑150 (MQ=255)
                                                                                                                                  aTTGAAACGGACGTGTGGCACTGCCGATGTTTTTATAAATATGTTAAACGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGCTGATTATTTATGGGTGACAATAAAAAGACAGTa                   >  2:393848/1‑150 (MQ=255)
                                                                                                                                  aTTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTa                   >  2:240490/1‑150 (MQ=255)
                                                                                                                                   ttGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTAt                  >  2:1008124/1‑150 (MQ=255)
                                                                                                                                   ttGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTAt                  >  2:634213/1‑150 (MQ=255)
                                                                                                                                     gAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTa                >  1:544360/1‑150 (MQ=255)
                                                                                                                                     gAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTa                >  2:377004/1‑150 (MQ=255)
                                                                                                                                     gAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATAAGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTa                >  2:252607/1‑150 (MQ=255)
                                                                                                                                     acaaCGGCCGTGCCACACGTCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTa                >  1:20921/3‑150 (MQ=255)
                                                                                                                                       aaCGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAat              <  2:55161/150‑1 (MQ=255)
                                                                                                                                        aCGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAatc             >  2:517337/1‑150 (MQ=255)
                                                                                                                                          ggACGTGTGGCCCGGCCGTTGTTTTTATAAATATCTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAatcat           >  2:451233/1‑150 (MQ=255)
                                                                                                                                          ggACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAatcat           <  2:1222890/150‑1 (MQ=255)
                                                                                                                                            aCGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATcc         <  1:1338303/150‑1 (MQ=255)
                                                                                                                                             cGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCa        <  2:765795/150‑1 (MQ=255)
                                                                                                                                             cGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCa        <  1:57365/150‑1 (MQ=255)
                                                                                                                                               tgtgGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATa      >  2:309436/1‑150 (MQ=255)
                                                                                                                                                gtgGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATag     <  2:923035/150‑1 (MQ=255)
                                                                                                                                                gtgGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTATGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATag     >  2:547969/1‑150 (MQ=255)
                                                                                                                                                  ggCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATagag   >  1:1002634/1‑150 (MQ=255)
                                                                                                                                                   gCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATAAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAgaga  >  1:309020/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AATAATGAATATAAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGA  >  NC_000913/1212330‑1212626

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