Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,431,811 G→T 5.4% A51E (GCG→GAG)  yrdD ← putative DNA topoisomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,431,8110GT 5.4% 306.9 / 3.9 91A51E (GCG→GAG) yrdDputative DNA topoisomerase
Reads supporting (aligned to +/- strand):  ref base G (43/43);  new base T (2/3);  total (45/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

GTCCGGTTTATCGATAAGTTCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCACCGTGAACAGTGCTGATTTCGCC  >  NC_000913/3431663‑3431960
                                                                                                                                                    |                                                                                                                                                     
gTCCGGTTTATCGATAAGTTCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGc                                                                                                                                                      <  2:1408469/150‑1 (MQ=255)
        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:441164/1‑150 (MQ=255)
                                            tAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACtgtg                                                                                                          <  1:476528/150‑1 (MQ=255)
                                            tAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACtgtg                                                                                                          <  2:745400/150‑1 (MQ=255)
                                            tAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACtgtg                                                                                                          <  1:1451845/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                                                                                      >  2:218531/1‑150 (MQ=255)
                                                                ccAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCg                                                                                      <  1:1616558/150‑1 (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                                                                        <  2:271666/150‑1 (MQ=255)
                                                                              gCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCg                                                                        <  1:1098961/150‑1 (MQ=255)
                                                                                    tACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCacac                                                                                                                       <  2:795419/97‑1 (MQ=255)
                                                                                    tACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCacac                                                                                                                       >  1:795419/1‑97 (MQ=255)
                                                                                     aCCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAAcc                                                                 >  1:280391/1‑150 (MQ=255)
                                                                                     aCCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAAcc                                                                 <  1:1185140/150‑1 (MQ=255)
                                                                                     aCCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAAcc                                                                 >  1:822135/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                                                  >  2:472725/1‑150 (MQ=255)
                                                                                                    caTGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCAc                                                  >  1:523046/1‑150 (MQ=255)
                                                                                                      tGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACtt                                                >  1:841907/1‑150 (MQ=255)
                                                                                                           ggCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGc                                           <  2:1161122/150‑1 (MQ=255)
                                                                                                           ggCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGc                                           <  2:441164/150‑1 (MQ=255)
                                                                                                                   tGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCa                                   <  2:520615/150‑1 (MQ=255)
                                                                                                                          ccAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTAc                            <  2:1603015/150‑1 (MQ=255)
                                                                                                                          ccAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTAc                            <  2:318350/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:773081/1‑150 (MQ=255)
                                                                                                                                           tGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCACCGTGAACAGTGCt           >  1:507848/1‑150 (MQ=255)
                                                                                                                                           tGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCACCGTGAACAGTGCt           >  2:231864/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)
                                                                                                                                                  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)
                                                                                                                                                    |                                                                                                                                                     
GTCCGGTTTATCGATAAGTTCGGTATGTTCGCATTCAGGGTAGTTAATGCAACCAATAAACATACCAAAGCGTCCCTGGCGTAATACCAGATTTGCGCCACATGCAGGGCAAACCTGCCCCTCCAGAACTTTGACGATATGTCCATCCGCTGAAGATTTCAGAGGACGGACGTAGTCACACGCCGGATACTGTGAGCATCCAAGAAACGGACCGTGTTTCCCGGATCGAATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCACCGTGAACAGTGCTGATTTCGCC  >  NC_000913/3431663‑3431960

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: