Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 532,945 C→G 7.2% Y70* (TAC→TAG) ‡ tufA → elongation factor Tu‑A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947532,9450CG 7.2% 197.5 / 6.0 69Y70* (TAC→TAG) ‡tufAelongation factor Tu‑A
Reads supporting (aligned to +/- strand):  ref base C (33/31);  new base G (3/2);  total (36/33)
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.22e-01

ACCACGGTAAGACCACTCTGACCGCAGCTCTGACTCGCGTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCACAAACCCGTGA  >  NC_002947/532800‑533088
                                                                                                                                                  |                                                                                                                                               
aCCACGGTAAGACCACTCTGACCGCAGCTCTGAGCTCGCGTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACaa                                                                                                                                               >  1:225293/1‑149 (MQ=21)
aCCACGGTAAGACCACTCTGACCGCAGCTCTGACTCGCGTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACaa                                                                                                                                               <  7:88695/148‑1 (MQ=34)
 ccACGGTAAGACCACTCTGACCGCAGCTCTGACTCGCGTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAAc                                                                                                                                              <  1:106617/148‑1 (MQ=34)
 ccACGGTAAGACCACTCTGACCGCAGCTCTGACTCGCGTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACt                                                                                                                                             >  2:351328/1‑149 (MQ=34)
    cGGTAAGACCACTCTGACCGCAGCTCTGACTCGCGTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGa                                                                                                                                          >  5:88620/1‑149 (MQ=34)
     ggTTAGAGCACTCTGTCGGCAGCGCTGACTCGCGTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGTGTTCGACAAGATCGTCTCGGCTCCGGAAGAAAAAGCGCGCGGTAGCACCATCACCACCGCTCACGTCGAGTTGTACTCGAc                                                                                                                                         <  5:17229/149‑1 (MQ=11)
      gTAAGACCACTCTGACCGCAGCTCTGACTCGCGTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTTGTACTCGAcg                                                                                                                                        >  3:30258/1‑148 (MQ=255)
          gACCACTCTGACCGCAGCTCTGACTCGCGTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTTGTACTCGAcgtgag                                                                                                                                    >  6:173651/1‑144 (MQ=255)
          gACCACTCTGACCGCAGCTCTGACTCGCGTTTGATCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTc                                                                                                                                    <  2:225293/149‑1 (MQ=37)
           aCCACTCTGACCGCAGCTCTGACTCGCGTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTTGTACTCGAcgtgagc                                                                                                                                   <  4:30258/149‑7 (MQ=11)
            ccACTCTGACCGCAGCTCTGACTCGCGTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGt                                                                                                                                  <  8:138387/149‑1 (MQ=37)
             cACTCTGACCGCAGCTCTGACTCGCGTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTc                                                                                                                                 <  4:140940/149‑1 (MQ=37)
                   gACCGCAGCTCTGACTCGCGTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACg                                                                                                                           >  4:132268/1‑149 (MQ=37)
                      cGCAGCTCTGACTCGCGTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTc                                                                                                                        >  8:296379/1‑149 (MQ=37)
                            tctGACTCGCGTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGtt                                                                                                                   >  4:125378/1‑148 (MQ=37)
                                aCTCGCGTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGAc                                                                                                                <  1:203480/147‑1 (MQ=37)
                                aCTCGCGTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTg                                                                                                              >  6:326396/1‑149 (MQ=37)
                                   tCGCGTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCa                                                                                                                                <  1:36675/129‑1 (MQ=38)
                                   tCGCGTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCa                                                                                                                                >  2:36675/1‑129 (MQ=38)
                                   tCGCGTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGcc                                                                                                            >  1:92506/1‑149 (MQ=37)
                                   tCGCGTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGcc                                                                                                            >  8:282467/1‑149 (MQ=37)
                                     gcgTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGcc                                                                                                            >  6:19225/1‑147 (MQ=37)
                                          tGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCa                                                                                                     >  5:352845/1‑149 (MQ=34)
                                             tCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTAGCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGc                                                                                                  >  6:27574/1‑149 (MQ=21)
                                                    ttttCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTAc                                                                                           <  6:352845/149‑1 (MQ=34)
                                                     tttCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTAc                                                                                           <  7:296379/148‑1 (MQ=34)
                                                      ttcggttcggCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGt                                                                                         >  1:81871/1‑149 (MQ=34)
                                                          gttcggCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTgaag                                                                                     <  5:19225/149‑1 (MQ=34)
                                                          gttcggCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTgaag                                                                                     <  8:167236/149‑1 (MQ=34)
                                                           ttcggCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTgaag                                                                                     >  4:223044/1‑148 (MQ=34)
                                                               gCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACAt                                                                                <  1:351328/149‑1 (MQ=34)
                                                                  aTCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGAt                                                                             <  2:81871/149‑1 (MQ=32)
                                                                         aGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGt                                                                      >  8:59855/1‑149 (MQ=32)
                                                                           ttCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGtgc                                                                    >  1:19943/1‑149 (MQ=32)
                                                                               aCAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTc                                                                                                      >  5:349966/1‑111 (MQ=32)
                                                                               aCAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTc                                                                                                      <  6:349966/111‑1 (MQ=32)
                                                                                    aTCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGAt                                                           >  5:101649/1‑149 (MQ=32)
                                                                                     tCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGtgct                                                                   <  2:40030/140‑1 (MQ=32)
                                                                                     tCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGtgct                                                                   >  1:40030/1‑140 (MQ=32)
                                                                                          tCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCAGGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTgaat                                                                                     <  2:187709/117‑2 (MQ=17)
                                                                                          tCGCCTCCGGAAGAAAAAGCGCGCGGTATCAACAACAACACCGCTCACGTCGAGTACAACTCGACCAATCCTCACAACGCTCACGTAGACTGCCCAGGTCACCCTGACTACGTgaag                                                                                     >  1:187709/1‑117 (MQ=11)
                                                                                             gCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGtgctgc                                                                 <  7:99819/134‑1 (MQ=32)
                                                                                             gCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGtgctgc                                                                 >  8:99819/1‑134 (MQ=32)
                                                                                               tCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTgaag                                                                                     <  4:40115/112‑1 (MQ=32)
                                                                                               tCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTgaag                                                                                     >  3:40115/1‑112 (MQ=32)
                                                                                                 cGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATc                                              <  2:19943/149‑1 (MQ=32)
                                                                                                 cGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATc                                              >  3:76886/1‑149 (MQ=32)
                                                                                                      gaaAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGt                                         >  1:17015/1‑149 (MQ=32)
                                                                                                      gaaAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGt                                         >  2:296685/1‑149 (MQ=32)
                                                                                                        aaaaGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGttt                                       <  2:92506/149‑1 (MQ=32)
                                                                                                          aaGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGtt                                        <  1:284496/146‑1 (MQ=32)
                                                                                                          aaGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGtt                                        >  2:284496/1‑146 (MQ=32)
                                                                                                           aGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCt                                    <  7:282467/149‑1 (MQ=32)
                                                                                                                      tCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGa                          <  7:59855/148‑1 (MQ=32)
                                                                                                                        aCCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCg                                  >  4:130183/1‑138 (MQ=32)
                                                                                                                        aCCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCg                                  <  3:130183/138‑1 (MQ=32)
                                                                                                                        aCCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATgg                       <  5:326396/149‑1 (MQ=32)
                                                                                                                        aCCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATgg                       >  3:144621/1‑149 (MQ=32)
                                                                                                                        aCCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATg                        <  3:125378/148‑1 (MQ=32)
                                                                                                                           aTCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTcc                    >  3:185566/1‑149 (MQ=32)
                                                                                                                              aaCACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGAt                 >  5:366871/1‑149 (MQ=32)
                                                                                                                                    gCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTc                     <  7:141136/139‑1 (MQ=32)
                                                                                                                                    gCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTc                     >  8:141136/1‑139 (MQ=32)
                                                                                                                                      tCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCACaaa         >  6:184320/1‑149 (MQ=34)
                                                                                                                                        aCGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCACaaa         <  5:27574/147‑1 (MQ=34)
                                                                                                                                          gTCGAGTTGTACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCACAAACCCg     >  5:321530/1‑149 (MQ=255)
                                                                                                                                             gAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCACAAACCCGt    <  3:223044/147‑1 (MQ=34)
                                                                                                                                             gAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCACAAACCCGTGa  <  4:76886/149‑1 (MQ=34)
                                                                                                                                             gAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCACAAACCCGTGa  <  5:184320/149‑1 (MQ=34)
                                                                                                                                                  |                                                                                                                                               
ACCACGGTAAGACCACTCTGACCGCAGCTCTGACTCGCGTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGAGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCACAAACCCGTGA  >  NC_002947/532800‑533088

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