Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 902,618 C→G 7.2% intergenic (+44/+32) trxB → / ← nadC thioredoxin reductase/quinolinate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947902,6180CG 7.2% 212.7 / 6.3 70intergenic (+44/+32)trxB/nadCthioredoxin reductase/quinolinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (28/36);  new base G (2/3);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.34e-01

GCTGACCACGTTTACCGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAGGGTCGTTTCGTT  >  NC_002947/902470‑902760
                                                                                                                                                    |                                                                                                                                              
gCTGACCACGTTTACCGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGcc                                                                                                                                                <  2:201656/149‑1 (MQ=255)
  tGACCACGTTTACCGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGcc                                                                                                                                                >  3:83189/1‑147 (MQ=255)
   gACCACGTTTACCGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAg                                                                                                                                             >  6:90379/1‑149 (MQ=255)
    aCCACGTTTACCGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAgc                                                                                                                                            <  6:222577/149‑1 (MQ=255)
      cACGTTTACCGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAgcgc                                                                                                                                          >  7:90667/1‑149 (MQ=255)
        cGTTTACCGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCg                                                                                                                                        <  2:3189/149‑1 (MQ=255)
        cGTTTACCGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCg                                                                                                                                        <  7:279056/149‑1 (MQ=255)
         gTTTACCGCCAGGCCATCACTACGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCAAAAAAAAAACCGCCAAACCCCCgg                                                                                                                                       >  2:322540/1‑149 (MQ=255)
             aCCGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGttt                                                                                                                                    <  2:25757/148‑1 (MQ=255)
              ccGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGGCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGttttt                                                                                                                                  <  4:69032/149‑1 (MQ=255)
               cGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGtttttt                                                                                                                                 <  4:68750/149‑1 (MQ=255)
                  cAGGCCATCACTTCGGCCGGTGCCGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGAGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTAtt                                                                                                                              <  7:273837/149‑1 (MQ=255)
                    ggCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCGCTTGGCCGGTTTTTTATTgg                                                                                                                            <  3:1568/149‑2 (MQ=39)
                    ggCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCGCTTGGCCGGTTTTTTATTgg                                                                                                                            <  7:277636/149‑2 (MQ=39)
                      ccATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCt                                                                                                                          <  6:13159/149‑1 (MQ=255)
                       cATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGACGGTTGTTTATTGCCTg                                                                                                                         >  6:300240/1‑149 (MQ=255)
                        aTCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGtttt                                                                                                                                   >  1:276383/1‑138 (MQ=255)
                        aTCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGtttt                                                                                                                                   <  2:276383/138‑1 (MQ=255)
                            cTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAAt                                                                                                                     >  7:309206/1‑148 (MQ=255)
                             ttCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAAt                                                                                                                     >  7:28198/1‑147 (MQ=255)
                               cGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAg                                                                                                                  <  5:90379/148‑1 (MQ=255)
                                ggCCGGTGCGGGCTGCATGGCGGCGCTTGATGTCGAGCGCTACCTGGAAGGCTTGGCCAACGCCTCGTGCTGATCTGAGAAAAAGGGGGGCAAAGCGCCCCCAATAAAAAACCCGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCt                                                                                                                <  3:224867/149‑1 (MQ=255)
                                 gCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGttt                                                                                                                                    <  3:351359/128‑1 (MQ=255)
                                 gCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGttt                                                                                                                                    >  4:351359/1‑128 (MQ=255)
                                   cGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTcac                                                                                                             <  1:186167/149‑1 (MQ=255)
                                    ggTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCaa                                                                                                                                              >  7:42320/1‑115 (MQ=255)
                                    ggTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCaa                                                                                                                                              <  8:42320/115‑1 (MQ=255)
                                    ggTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTcaca                                                                                                            <  8:90667/149‑1 (MQ=255)
                                    ggTGC‑GGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAg                                                                                                           <  7:67802/149‑1 (MQ=255)
                                     gTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAg                                                                                                           >  3:169133/1‑149 (MQ=255)
                                             tGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTc                                                                                                   <  2:142021/149‑1 (MQ=255)
                                               cATGGCGGCGCTGGATGTCGCTCGCTACCTGGACGTTTTGGCCTACGCATACTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCGCCAATAATAAACAGGCCAAGCGCCTGTTATTTATTGCCTGCAATCAGCTCACAGGCGCAGTc                                                                                                   <  4:26612/147‑1 (MQ=255)
                                               cATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGc                                                                                                 <  4:83189/149‑1 (MQ=255)
                                               cATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTGTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCg                                                                                                  >  6:181932/1‑148 (MQ=255)
                                                    cggcgCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCTTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGa                                                                                            <  3:201703/149‑1 (MQ=255)
                                                        gcTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAgg                                                                                        >  5:265412/1‑149 (MQ=255)
                                                        gcTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGAAATCAGCTCACAGGCTCAGTCGCATCGACAgg                                                                                        >  4:217300/1‑149 (MQ=255)
                                                         cTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGt                                                                                       <  3:377002/149‑1 (MQ=255)
                                                             aTGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCAc                                                                                   <  7:233926/149‑1 (MQ=255)
                                                                       cTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTcac                                                                          >  8:104458/1‑148 (MQ=255)
                                                                        tACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACAt                                                                        >  6:299449/1‑149 (MQ=255)
                                                                             ggACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTg                                                                   >  5:176880/1‑149 (MQ=255)
                                                                              gACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTgg                                                                  >  4:83411/1‑149 (MQ=255)
                                                                              gACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTg                                                                   <  6:345304/148‑1 (MQ=255)
                                                                                cGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTc                                                                <  8:28198/149‑1 (MQ=255)
                                                                                  gTTTGGCCAACGGCTTGGTCTGATCTGATAAAAACGGGGGTAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCAt                                                              <  1:322540/149‑1 (MQ=255)
                                                                                         cAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCGCTTGGCCGGTTTTTTAtt                                                                                                                              >  5:177573/1‑78 (MQ=25)
                                                                                         cAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCGCTTGGCCGGTTTTTTAtt                                                                                                                              <  6:177573/78‑1 (MQ=25)
                                                                                           aCGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAAt                                                     >  2:214656/1‑149 (MQ=255)
                                                                                              ccTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCga                                                  <  7:104458/149‑1 (MQ=255)
                                                                                                 cGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGAc                                                                                           >  6:201355/1‑105 (MQ=255)
                                                                                                 cGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGAc                                                                                           <  1:312367/105‑1 (MQ=255)
                                                                                                 cGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGAc                                                                                           <  5:201355/105‑1 (MQ=255)
                                                                                                 cGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGAc                                                                                           >  2:312367/1‑105 (MQ=255)
                                                                                                  gTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCgaga                                                <  5:97497/147‑1 (MQ=255)
                                                                                                  gTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGAt                                               <  5:299449/148‑1 (MQ=255)
                                                                                                   ttCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATg                                              <  4:169133/148‑1 (MQ=255)
                                                                                                      tgatctgaGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGACACCAATCGAGATGTAGt                                          <  3:83411/149‑1 (MQ=255)
                                                                                                       gatctgaGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAAGCGAGATGTAGTc                                         >  6:223363/1‑149 (MQ=255)
                                                                                                        atctgaGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTc                                         >  2:287133/1‑148 (MQ=255)
                                                                                                          ctgaGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTTGTCCAc                                      >  1:282498/1‑149 (MQ=255)
                                                                                                           tgaGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCAcc                                     >  8:152486/1‑149 (MQ=255)
                                                                                                             agaAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATAAGATCACAGGCTCAGTCGCATCGAGAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCAcccc                                   >  5:274853/1‑149 (MQ=255)
                                                                                                                        tGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCaa                        >  2:382275/1‑149 (MQ=255)
                                                                                                                              gcgcCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCAccc                  <  7:152486/149‑1 (MQ=255)
                                                                                                                                  ccccAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAg              >  2:278035/1‑149 (MQ=255)
                                                                                                                                  aggcAATAAAAAACCGGCGCTTGGCCGGTTTGTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAg              >  1:380344/4‑149 (MQ=39)
                                                                                                                                          aaaaaCCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAGGGtcgttt      >  4:43669/1‑149 (MQ=255)
                                                                                                                                           aaaaCCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAGGGtcgtttc     <  5:97449/149‑1 (MQ=255)
                                                                                                                                            aaaCCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAGGGtcgtttcg    >  6:20695/1‑149 (MQ=255)
                                                                                                                                              aCCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCATGGtcgtttcgtt  <  3:153050/149‑1 (MQ=255)
                                                                                                                                              aCCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAGGGtcgtttcgtt  <  1:382275/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
GCTGACCACGTTTACCGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAGGGTCGTTTCGTT  >  NC_002947/902470‑902760

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