Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,324,959 C→T 5.9% intergenic (‑48/+42) folP ← / ← ftsH 7,8‑dihydropteroate synthase/protease, ATP‑dependent zinc‑metallo

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,324,9590CT 5.9% 546.0 / 13.5 167intergenic (‑48/+42)folP/ftsH7,8‑dihydropteroate synthase/protease, ATP‑dependent zinc‑metallo
Reads supporting (aligned to +/- strand):  ref base C (71/86);  new base T (2/8);  total (73/94)
Fisher's exact test for biased strand distribution p-value = 1.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 170 total aligned reads displayed.

CGTGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTA  >  NC_000913/3324813‑3325106
                                                                                                                                                  |                                                                                                                                                   
   gCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCgg                                                                                                                                                 <  3:393478/148‑1 (MQ=255)
   gCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCgg                                                                                                                                                 >  3:539620/1‑148 (MQ=255)
   gCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCgg                                                                                                                                                 <  3:660251/148‑1 (MQ=255)
     caccaTCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGGGCGTGCTCCGGGGt                                                                                                                                              <  2:368639/149‑1 (MQ=255)
     caccaTCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGt                                                                                                                                              >  2:638854/1‑149 (MQ=255)
      accaTCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGGGCGTGCTACGGGGtt                                                                                                                                             <  2:588377/149‑1 (MQ=255)
        caTCCGAAAAGGAATCAGGCGTGACGTTGAGAATCCCCAATACGGGGGGATGGCTAATGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTCATAAGAAAAAACCCCGGGGCGTGCTCCGGGGtttt                                                                                                                                           <  1:93238/149‑1 (MQ=255)
          tCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGGGCGTGCTCCGGGGttttt                                                                                                                                          <  3:146798/148‑2 (MQ=255)
          tCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTc                                                                                                                                          <  4:558402/148‑1 (MQ=255)
             gAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGGGCGTGCTCCGGGGtttt                                                                                                                                           <  1:980084/144‑1 (MQ=255)
                 aGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACaaa                                                                                                                                  >  2:661274/1‑149 (MQ=255)
                     aTCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATAc                                                                                                                               <  4:912377/148‑1 (MQ=255)
                     aTCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACa                                                                                                                              <  1:661274/149‑1 (MQ=255)
                     aTCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACa                                                                                                                              >  2:891822/1‑149 (MQ=255)
                        aGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTc                                                                                                                           <  2:538940/149‑1 (MQ=255)
                         ggCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTc                                                                                                                           >  3:277151/1‑148 (MQ=255)
                                  ttGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCg                                                                                                                 <  1:792270/149‑1 (MQ=255)
                                    gAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCgg                                                                                                                <  4:568810/148‑1 (MQ=255)
                                      ggATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGaa                                                                                                             <  2:970232/149‑1 (MQ=255)
                                      ggATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGaa                                                                                                             <  1:199223/149‑1 (MQ=255)
                                      ggATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGaa                                                                                                             >  2:1059983/1‑149 (MQ=255)
                                      ggATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGa                                                                                                              <  4:835134/148‑1 (MQ=255)
                                         tCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAactt                                                                                                          <  1:105884/149‑1 (MQ=255)
                                         tCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAact                                                                                                           >  4:144203/1‑148 (MQ=255)
                                          ccccATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAactt                                                                                                          <  4:247412/148‑1 (MQ=255)
                                             cATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCg                                                                                                                 <  3:821651/138‑1 (MQ=255)
                                             cATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAacttact                                                                                                       >  4:302501/1‑148 (MQ=255)
                                                tACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTc                                                                                                   <  1:276473/149‑1 (MQ=255)
                                                      aGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTaa                                                                                             <  2:663579/149‑1 (MQ=255)
                                                      aGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTa                                                                                              <  3:1050212/148‑1 (MQ=255)
                                                       ggATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAAc                                                                                            <  2:822532/149‑1 (MQ=255)
                                                               aaGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAactt                                                                                                          >  4:331860/1‑127 (MQ=255)
                                                               aaGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGa                                                                                    <  2:414294/149‑1 (MQ=255)
                                                               aaGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTAGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAactt                                                                                                          <  3:331860/127‑1 (MQ=255)
                                                                   tCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATg                                                                                <  1:974564/149‑1 (MQ=255)
                                                                    ccAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATgg                                                                               <  2:703256/149‑1 (MQ=255)
                                                                     cAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGgt                                                                              <  1:143329/149‑1 (MQ=255)
                                                                          aaGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTAc                                                                         >  2:324864/1‑149 (MQ=255)
                                                                           aGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTAc                                                                         >  4:229080/1‑148 (MQ=255)
                                                                            gTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTAccc                                                                       >  2:286518/1‑149 (MQ=255)
                                                                             tACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGt                                                                                                    >  3:595391/1‑119 (MQ=255)
                                                                             tACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATgg                                                                               <  3:71192/140‑1 (MQ=255)
                                                                             tACCCTGCGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTAccc                                                                       >  4:768927/1‑148 (MQ=255)
                                                                              aCCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCg                                                                      >  4:954211/1‑148 (MQ=255)
                                                                               cccTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCggg                                                                    >  1:655768/1‑149 (MQ=255)
                                                                                     gCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTaa                                                                                             >  1:605152/1‑118 (MQ=255)
                                                                                     gCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCgg                                                               <  4:418868/148‑1 (MQ=255)
                                                                                     gCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCgg                                                               >  3:731537/1‑148 (MQ=255)
                                                                                     gCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGc                                                              <  1:324864/149‑1 (MQ=255)
                                                                                     gCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGc                                                              <  1:781818/149‑1 (MQ=255)
                                                                                        aaGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGt                                                            >  4:340397/1‑148 (MQ=255)
                                                                                         agagTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGgtgt                                                                            <  4:81935/131‑1 (MQ=255)
                                                                                         agagTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGgtgt                                                                            >  3:81935/1‑131 (MQ=255)
                                                                                           agTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGct                                                                                        >  3:871942/1‑117 (MQ=255)
                                                                                              ttCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCgg                                                      >  3:304506/1‑148 (MQ=255)
                                                                                               tCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGt                                                     <  4:435674/148‑1 (MQ=255)
                                                                                                     tgTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTtcatca                                               >  3:542326/1‑148 (MQ=255)
                                                                                                     tgTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCaa                                              <  1:286518/149‑1 (MQ=255)
                                                                                                      gTTATGTCTGGTATGATTTGATAAGAAAAAACTCCTGAGCACGCTCCGGGGTTTTCGGTCCAAATACATTCATCTGATGCGGGAACTTACTCGTTCCCCAACTGCTCTGACACGCTTTTACCCTGTTTCCCCGCACGCGGTTCATCaa                                              <  4:111531/148‑1 (MQ=255)
                                                                                                          tCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCgg                                          >  3:874369/1‑148 (MQ=255)
                                                                                                          tCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCgg                                          <  3:302501/148‑1 (MQ=255)
                                                                                                          tCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCgg                                          <  4:255099/148‑1 (MQ=255)
                                                                                                          tCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGa                                         <  2:981386/149‑1 (MQ=255)
                                                                                                             cTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACga                                      >  2:154343/1‑149 (MQ=255)
                                                                                                             cTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTCACATGGTGTTACCCGGGTTCGCGTACGCGGTTCATCAACCGGACga                                      >  3:12438/1‑148 (MQ=255)
                                                                                                              tGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACgag                                     >  1:858718/1‑149 (MQ=255)
                                                                                                                gTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACgagga                                   <  2:918800/149‑1 (MQ=255)
                                                                                                                gTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACgagga                                   <  1:154343/149‑1 (MQ=255)
                                                                                                                gTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACgagg                                    <  3:144203/148‑1 (MQ=255)
                                                                                                                  aTGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACgaggag                                  <  3:768927/148‑1 (MQ=255)
                                                                                                                  aTGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACgaggag                                  <  3:945742/148‑1 (MQ=255)
                                                                                                                      aTTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATAc                                                                                                                               >  4:907442/1‑51 (MQ=255)
                                                                                                                      aTTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATAc                                                                                                                               <  3:907442/51‑1 (MQ=255)
                                                                                                                       ttGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCttt                             >  4:489445/1‑148 (MQ=255)
                                                                                                                         gATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTgg                           <  4:624270/148‑1 (MQ=255)
                                                                                                                         gATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATAGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTgg                           <  4:731537/148‑1 (MQ=255)
                                                                                                                            aaGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACt                        >  4:398400/1‑148 (MQ=255)
                                                                                                                            aaGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTa                       >  2:517839/1‑149 (MQ=255)
                                                                                                                             aGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTa                       >  3:376496/1‑148 (MQ=255)
                                                                                                                             aGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTAc                      <  1:1046255/149‑1 (MQ=255)
                                                                                                                              gAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTAcc                     >  2:618724/1‑149 (MQ=255)
                                                                                                                                aaaaaCCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCa                    >  3:915693/1‑148 (MQ=255)
                                                                                                                                 aaaaCCCCGGAGCACGCCCCGAGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCAtt                  >  2:258086/1‑149 (MQ=255)
                                                                                                                                    aCCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGt                <  4:824175/148‑1 (MQ=255)
                                                                                                                                    aCCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGt                <  3:973299/148‑1 (MQ=255)
                                                                                                                                     ccccGGGGCGTGCTCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCg              >  2:766207/1‑149 (MQ=255)
                                                                                                                                     ccccGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTc               >  4:706633/1‑148 (MQ=255)
                                                                                                                                      cccGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCg              <  4:542326/148‑1 (MQ=255)
                                                                                                                                       ccGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGc             >  4:902065/1‑148 (MQ=255)
                                                                                                                                       ccGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGc             >  4:923199/1‑148 (MQ=255)
                                                                                                                                        cGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGcc            >  3:245987/1‑148 (MQ=255)
                                                                                                                                        cGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGcc            <  3:157032/148‑1 (MQ=255)
                                                                                                                                              aCGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTg     <  2:295968/149‑1 (MQ=255)
                                                                                                                                               tgCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTg     >  4:203363/2‑148 (MQ=255)
                                                                                                                                                gCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGtt   <  1:109404/149‑1 (MQ=255)
                                                                                                                                                gCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGt    <  4:277151/148‑1 (MQ=255)
                                                                                                                                                gCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGt    <  3:416247/148‑1 (MQ=255)
                                                                                                                                                gCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGt    <  3:229080/148‑1 (MQ=255)
                                                                                                                                                 ccccGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGtt   <  3:489445/148‑1 (MQ=255)
                                                                                                                                                  cccGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTa  <  3:398400/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CGTGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTA  >  NC_000913/3324813‑3325106

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: