Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,431,817 G→A 5.1% S49F (TCT→TTT) ‡ yrdD ← putative DNA topoisomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,431,8170GA 5.1% 317.4 / 3.9 98S49F (TCT→TTT) ‡yrdDputative DNA topoisomerase
Reads supporting (aligned to +/- strand):  ref base G (48/44);  new base A (2/3);  total (50/48)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01

TATCGATAAGTTCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCACCGTGAACAGTGCTGATTTCGCCATAAC  >  NC_000913/3431671‑3431965
                                                                                                                                                  |                                                                                                                                                    
tATCGATAAGTTCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAAAAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGAtt                                                                                                                                                   <  1:1233632/150‑1 (MQ=255)
 aTCGATAAGTTCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACATTGACGATATGTCCATCCGCTGAAGAttt                                                                                                                                                  >  2:1083846/1‑150 (MQ=255)
  tCGATACGTTCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACAAAAGCGTCCCTGGCCTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCATCCAGACCTTTGACGATATGTCCATCCGCTGAAGATTTc                                                                                                                                                 <  1:326603/150‑1 (MQ=255)
    gATAAGTTCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCag                                                                                                                                               >  1:789193/1‑150 (MQ=255)
    gATAAGTTCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCag                                                                                                                                               >  1:1603015/1‑150 (MQ=255)
       aaGTTCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAgg                                                                                                                                            <  2:450191/150‑1 (MQ=255)
             ggTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTCACGATATGTCCATCCGCTGAAGATTTCAGAggacggac                                                                                                                                      >  2:263964/1‑150 (MQ=255)
              gTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCGTCCTCTGAAATCTTCAGCggatggaca                                                                                                                                     <  2:1602316/150‑2 (MQ=255)
               tATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGt                                                                                                                                    <  2:914447/150‑1 (MQ=255)
                 tGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAg                                                                                                                                  <  2:789193/150‑1 (MQ=255)
                 tGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGGAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCGTCCTCTGAAATCTTCAGCGGATGGACATat                                                                                                                                  <  2:442536/150‑2 (MQ=255)
                   ttCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTc                                                                                                                                >  1:747228/1‑150 (MQ=255)
                    tCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTca                                                                                                                               <  2:209971/150‑1 (MQ=255)
                     cGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCGTCCTCTGAAATCTTCAGCGGATGGACATatcgtc                                                                                                                              <  1:1457390/150‑6 (MQ=38)
                     cGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTcac                                                                                                                              <  2:835257/150‑1 (MQ=255)
                          tCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGcc                                                                                                                         <  2:432384/150‑1 (MQ=255)
                             gggTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGa                                                                                                                      >  1:1465248/1‑150 (MQ=255)
                              ggTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGAt                                                                                                                     <  2:935044/150‑1 (MQ=255)
                                 aGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACt                                                                                                                  >  2:1521267/1‑150 (MQ=255)
                                  gTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACtg                                                                                                                 >  1:1161122/1‑150 (MQ=255)
                                    tAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACtgtg                                                                                                               <  1:1451845/150‑1 (MQ=255)
                                    tAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACtgtg                                                                                                               <  2:745400/150‑1 (MQ=255)
                                    tAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACtgtg                                                                                                               >  1:441164/1‑150 (MQ=255)
                                    tAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACtgtg                                                                                                               <  1:476528/150‑1 (MQ=255)
                                     aaTGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGa                                                                                                              <  2:729993/150‑1 (MQ=255)
                                        gCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCa                                                                                                           >  2:484920/1‑150 (MQ=255)
                                        gCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCa                                                                                                           <  1:224319/150‑1 (MQ=255)
                                         cAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAg                                                                                                             >  1:437164/1‑147 (MQ=255)
                                         cAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAg                                                                                                             <  2:437164/147‑1 (MQ=255)
                                         cAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCAt                                                                                                          >  1:271666/1‑150 (MQ=255)
                                           aCCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATcc                                                                                                        >  1:1239965/1‑150 (MQ=255)
                                             cAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGTGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCaa                                                                                                      >  2:1590221/1‑150 (MQ=255)
                                                tAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGaa                                                                                                   >  1:71192/1‑150 (MQ=255)
                                                    cATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACgg                                                                                               <  1:100650/150‑1 (MQ=255)
                                                       aCCAAAGCGTCTCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGAcc                                                                                            >  2:1507600/1‑150 (MQ=255)
                                                        ccAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCg                                                                                           <  1:1616558/150‑1 (MQ=255)
                                                        ccAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCg                                                                                           >  2:218531/1‑150 (MQ=255)
                                                        ccAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCg                                                                                           >  2:866032/1‑150 (MQ=255)
                                                             gCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGttt                                                                                      >  2:1290024/1‑150 (MQ=255)
                                                                  ccTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCgg                                                                                 <  1:484920/150‑1 (MQ=255)
                                                                     ggCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATc                                                                              <  2:1465248/150‑1 (MQ=255)
                                                                     ggCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATc                                                                              <  2:1239965/150‑1 (MQ=255)
                                                                      gCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCg                                                                             <  2:71192/150‑1 (MQ=255)
                                                                      gCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCg                                                                             <  1:1098961/150‑1 (MQ=255)
                                                                      gCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCg                                                                             <  2:271666/150‑1 (MQ=255)
                                                                            tACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCacac                                                                                                                            >  1:795419/1‑97 (MQ=255)
                                                                            tACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCacac                                                                                                                            <  2:795419/97‑1 (MQ=255)
                                                                             aCCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAAcc                                                                      <  1:1185140/150‑1 (MQ=255)
                                                                             aCCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAAcc                                                                      >  1:822135/1‑150 (MQ=255)
                                                                             aCCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAAcc                                                                      >  1:280391/1‑150 (MQ=255)
                                                                              ccAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCa                                                                     >  2:209761/1‑150 (MQ=255)
                                                                                  aTTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTc                                                                 <  2:559162/150‑1 (MQ=255)
                                                                                   tttGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGTCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCa                                                                >  2:895395/1‑150 (MQ=255)
                                                                                       cgcCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGccc                                                            >  1:318350/1‑150 (MQ=255)
                                                                                        gcCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGcccc                                                           <  2:1308254/150‑1 (MQ=255)
                                                                                          cacaTGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGc                                                         >  1:1591886/1‑150 (MQ=255)
                                                                                            caTGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCAc                                                       >  1:523046/1‑150 (MQ=255)
                                                                                            caTGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCAc                                                       >  2:472725/1‑150 (MQ=255)
                                                                                              tGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACtt                                                     >  1:841907/1‑150 (MQ=255)
                                                                                                   ggCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGc                                                <  2:441164/150‑1 (MQ=255)
                                                                                                   ggCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGc                                                <  2:1161122/150‑1 (MQ=255)
                                                                                                           tGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCa                                        <  2:520615/150‑1 (MQ=255)
                                                                                                                  ccAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTAc                                 <  2:318350/150‑1 (MQ=255)
                                                                                                                  ccAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTAc                                 <  2:1603015/150‑1 (MQ=255)
                                                                                                                      aaCTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCAc                             >  2:38420/1‑150 (MQ=255)
                                                                                                                        cTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCACCg                           >  2:542825/1‑150 (MQ=255)
                                                                                                                          ttGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCACCGTg                         >  2:566567/1‑150 (MQ=255)
                                                                                                                            gACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGTACTCATTATTACGCAACGTGaa                       >  2:614365/1‑150 (MQ=255)
                                                                                                                             aCGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCACCGTGAAc                      <  2:841907/150‑1 (MQ=255)
                                                                                                                               gATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATAACGCACCGTGAACAg                    <  2:265346/150‑1 (MQ=255)
                                                                                                                                 tatGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCACCGTGAACAGTg                  >  1:572658/1‑150 (MQ=255)
                                                                                                                                 tatGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCACCGTGAACAGTg                  >  2:1447341/1‑150 (MQ=255)
                                                                                                                                  atGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGAGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCACCGTGAACAGTGc                 <  1:866032/150‑1 (MQ=255)
                                                                                                                                   tGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCACCGTGAACAGTGCt                >  2:231864/1‑150 (MQ=255)
                                                                                                                                   tGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCACCGTGAACAGTGCt                >  2:773081/1‑150 (MQ=255)
                                                                                                                                   tGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCACCGTGAACAGTGCt                >  1:507848/1‑150 (MQ=255)
                                                                                                                                   tGTCCATCAGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTATGGGCAGGACTCATTATTACGCACCGGGAACAGTGCt                >  2:71865/1‑150 (MQ=255)
                                                                                                                                    gtccgtccTCTGAAATCTTCAGCGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCACCGTGAACAGTGCTg               >  2:1131172/1‑150 (MQ=255)
                                                                                                                                    gTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCACCGTGAACAGTGCTg               <  2:572658/150‑1 (MQ=255)
                                                                                                                                    gTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCACCGTGAACAGTGCTg               <  2:747228/150‑1 (MQ=255)
                                                                                                                                     tCCGTCCTCTGAAATCTTCAGCGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCACCGTGAACAGTGCTGa              >  2:222218/1‑150 (MQ=255)
                                                                                                                                     gccATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCACCGTGAACAGTGCTGa              <  1:263964/149‑1 (MQ=255)
                                                                                                                                      ccATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCACCGTGAACAGTGCTGAt             >  2:1133842/1‑150 (MQ=255)
                                                                                                                                      ccATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACAAGCCGGAAACTGTGAGCATCAAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGAACTTTGGGCAGGACACATTATTACGCAACGTGAACAGTGCTGaa             >  2:932368/1‑149 (MQ=255)
                                                                                                                                       cATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCACCGTGAACAGTGCTGAtt            >  1:197150/1‑150 (MQ=255)
                                                                                                                                       cATCCGCTGAAGATTTCAGAGGACGGACCTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCCTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGAAGGACTCATTATAACTCAACGTGAAAAGATCTGAtt            >  2:164759/1‑150 (MQ=255)
                                                                                                                                         tCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCACCGTGAACAGTGCTGATTTc          <  1:1290024/150‑1 (MQ=255)
                                                                                                                                         ctctcTGAATATTTCAGAGCACCGCCCTAGTCACACCCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCACCGTGAACAGTGCTGATTTc          <  1:359666/146‑1 (MQ=255)
                                                                                                                                          ccGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCACCGTGAACAGTGCTGATTTCg         <  1:780852/150‑1 (MQ=255)
                                                                                                                                           cGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCACCGTGAACAGTGCTGATTTCGc        <  2:146571/150‑1 (MQ=255)
                                                                                                                                           cGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCACCGTGAACAGTGCTGATTTCGc        <  1:1083846/150‑1 (MQ=255)
                                                                                                                                            gCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCACCGTGAACAGTGCTGATTTCGcc       >  2:128944/1‑150 (MQ=255)
                                                                                                                                               gAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCACCGTGAACAGTGCTGATTTCGCCATa    >  2:656578/1‑150 (MQ=255)
                                                                                                                                               gAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCACCGTGAACAGTGCTGATTTCGCCATa    >  1:1318305/1‑150 (MQ=255)
                                                                                                                                               gAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCACCGTGAACAGTGCTGATTTCGCCATa    >  2:1263627/1‑150 (MQ=255)
                                                                                                                                                 aGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCACCGTGAACGGTGCTGATTTCGCCATaac  <  2:822135/150‑1 (MQ=255)
                                                                                                                                                 aGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCACCGTGAACAGTGCTGATTTCGCCATaac  >  2:717889/1‑150 (MQ=255)
                                                                                                                                                 aGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCACCGTGAACAGTGCTGATTTCGCCATaac  <  1:566567/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TATCGATAAGTTCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCACCGTGAACAGTGCTGATTTCGCCATAAC  >  NC_000913/3431671‑3431965

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: