Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 902,622 C→G 7.7% intergenic (+48/+28) trxB → / ← nadC thioredoxin reductase/quinolinate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947902,6220CG 7.7% 196.0 / 6.4 65intergenic (+48/+28)trxB/nadCthioredoxin reductase/quinolinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (26/34);  new base G (2/3);  total (28/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.54e-01

ACCACGTTTACCGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAGGGTCGTTTCGTT  >  NC_002947/902474‑902760
                                                                                                                                                    |                                                                                                                                          
aCCACGTTTACCGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAgc                                                                                                                                            <  6:222577/149‑1 (MQ=255)
  cACGTTTACCGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAgcgc                                                                                                                                          >  7:90667/1‑149 (MQ=255)
    cGTTTACCGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCg                                                                                                                                        <  7:279056/149‑1 (MQ=255)
    cGTTTACCGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCg                                                                                                                                        <  2:3189/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                                                                                                                                    >  4:351359/1‑128 (MQ=255)
                             gCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGttt                                                                                                                                    <  3:351359/128‑1 (MQ=255)
                               cGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTcac                                                                                                             <  1:186167/149‑1 (MQ=255)
                                ggTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTcaca                                                                                                            <  8:90667/149‑1 (MQ=255)
                                ggTGCGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAg                                                                                                           <  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                                                                                                                              <  6:177573/78‑1 (MQ=25)
                                                                                     cAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCGCTTGGCCGGTTTTTTAtt                                                                                                                              >  5:177573/1‑78 (MQ=25)
                                                                                       aCGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAAt                                                     >  2:214656/1‑149 (MQ=255)
                                                                                          ccTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCga                                                  <  7:104458/149‑1 (MQ=255)
                                                                                             cGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGAc                                                                                           <  1:312367/105‑1 (MQ=255)
                                                                                             cGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGAc                                                                                           >  2:312367/1‑105 (MQ=255)
                                                                                             cGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGAc                                                                                           <  5:201355/105‑1 (MQ=255)
                                                                                             cGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGAc                                                                                           >  6:201355/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)
                                                                                                                                                    |                                                                                                                                          
ACCACGTTTACCGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAGGGTCGTTTCGTT  >  NC_002947/902474‑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: