Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 494,830 Δ1 bp 6.2% coding (591/768 nt) djlA → DnaJ‑like protein DjlA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947494,8300G. 6.2% 202.9 / 6.6 64coding (591/768 nt)djlADnaJ‑like protein DjlA
Reads supporting (aligned to +/- strand):  ref base G (33/27);  new base . (2/2);  total (35/29)
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.28e-01

TGAGTTGCTGCTGGACTGGGGGCAGAAGCTGGGCATGAGCCGGCGCCAGGTGCAGGCCATGTCCCTGGAGTACGAGCCGCGCAAGGCCACGGGAGCAGAAAGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCAGGTACCGGTGCCAGCGAGGCGCAGGTGCGTGAAGCAACCGCGCGTACCCGTGAGTTGCACCAGGCCTATG  >  NC_002947/494686‑494978
                                                                                                                                                |                                                                                                                                                    
tGAGTTGCTGCTGGACTGGGGGCAGAAGCTGGGCATGAGCCGGCGCCAGGTGCAGGCCATGTCCCTGGAGTACGAGCCGCGCAAGGCCACGGGAGCAGAAAGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACg                                                                                                                                                  >  1:69764/1‑149 (MQ=255)
  aGTTGCTGCTGGACTGGGGGCAGAAGCTGGGCATGAGCCGGCGCCAGGTGCAGGCCATGTCCCTGGAGTACGAGCCGCGCAAGGCCACGGGAGCAGAAAGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGc                                                                                                                                                 <  4:259020/148‑1 (MQ=255)
        tgctgGACTGGGGGTAGAAGCTGGGCATGAGCCGGCGCCAGGTGCAGGCCATGTCCCTGGAGTACGAGCCGCGCAAGGCCACGGGAGCAGAAAGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAgacacc                                                                                                                                          >  7:273379/1‑149 (MQ=255)
        tgctgGACTGGGGGCAGAAGCTGGGCATGACCCGGCGCCAGGTGCAGGCCATGTCCCTGGAGTACGAGCCGCGCAAGGCCACGGGAGCAGAAAGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAgacacc                                                                                                                                          >  1:26428/1‑149 (MQ=255)
            ggACTGGGGGCAGAAGCTGGGCATGAGCCGGCGCCAGGTGCAGGCCATGTCCCTGGAGTACGAGCCGCGCAAGGCCACGGGAGCAGAAAGCGTGCCCATGACCTACGCGGCAGCGTTGTGCTTGCTGGCGGTGGACGCAgacaccgaca                                                                                                                                      >  5:199303/1‑149 (MQ=255)
             gACTGGGGGCAGAAGCTGGGCATGAGCCGGCGCCAGGTGCAGGCCATGTCCCTGGAGTACGAGCCGCGCAAGGCCACGGGAGCAGAAAGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAgacaccgacac                                                                                                                                     >  7:93160/1‑149 (MQ=255)
                tGGGGGCAGAAGCTGGGCATGAGCCGGCGCCAGGTGCAGGCCATGTCCCTGGAGTACGAGCCGCGCAAGGCCACGGGAGCAGAAAGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGAcaccgacaccga                                                                                                                                  >  8:146349/1‑149 (MQ=255)
                  ggggCAGAAGCTGGGCATGAGCCGGCGCCAGGTGCAGGCCATGTCCCTGGAGTACGAGCCGCGCAAGGCCACGGGAGCAGAAAGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGAt                                                                                                                                 <  7:92572/148‑1 (MQ=255)
                     gCAGAAGCTGGGCATGAGCCGGCGCCAGGTGCAGGCCATGTCCCTGGAGTACGAGCCGCGCAAGGCCACGGGAGCAGAAAGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAgg                                                                                                                             >  5:209230/1‑149 (MQ=255)
                       aGAAGCTGGGCATGAGCCGGCGCCAGGTGCAGGCCATGTCCCTGGAGTACGAGCCGCGCAAGGCCACGGGAGCAGAAAGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTc                                                                                                                           <  5:312388/149‑1 (MQ=255)
                               ggCATGAGCCGGCGCCAGGTGCAGGCCATGTCCCTGGAGTACGAGCCGCGCAAGGCCACGGGAGCAGAAAGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGc                                                                                                                   <  8:97497/149‑1 (MQ=255)
                               ggCATGAGCCGGCGCCAGGTGCAGGCCATGTCCCTGGAGTACGAGCCGCGCAAGGCCACGGGAGCAGAAAGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGAAAAGGTCAAGCAgg                                                                                                                    <  7:119767/148‑1 (MQ=255)
                                gCATGAGCCGGCGCCAGGTGCAGGCCATGTCCCTGGAGTACGAGCCGCGCAAGGCCACGGGAGCAGAAAGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCg                                                                                                                  <  3:260751/149‑1 (MQ=255)
                                  aTGAGCCGGCGCCAGGTGCAGGCCATGTCCCTGGAGTACGAGCCGCGCAAGGCCACGGGAGCAGAAAGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGc                                                                                                                   >  4:64890/1‑146 (MQ=255)
                                  aTGAGCCGGCGCCAGGTGCAGGCCATGTCCCTGGAGTACGAGCCGCGCAAGGCCACGGGAGCAGAAAGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGc                                                                                                                   <  3:64890/146‑1 (MQ=255)
                                   tGAGCCGGCGCCAGGTGCAGGCCATGTCGCTGGAGTACGAGCCGCGCAAGGCCACGGGAGCAGAAAGCGGGCCCATGACCTACGCGGCAGCGTTGCGCTTTCTGGCGGTGGACGCAGACACCGACACCGAGAAGGTCAAGCAGGCGTAt                                                                                                               >  7:273752/1‑149 (MQ=255)
                                        cGGCGCCAGGTGCAGGCCATGTCCCTGGAGTACGAGCCGCGCAAGGCCACGGGAGCAGAAAGCGTGCCCATGACCTACGCGGCAGCGTTGTGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTAtcgtcg                                                                                                          <  6:199303/149‑1 (MQ=255)
                                          gcgcCAGGTGCAGGCCATGTCCCTGGAGTACGAGCCGCGCAAGGCCACGGGAGCAGAAAGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGAc                                                                                                                                                   >  6:259916/1‑106 (MQ=255)
                                          gcgcCAGGTGCAGGCCATGTCCCTGGAGTACGAGCCGCGCAAGGCCACGGGAGCAGAAAGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGAc                                                                                                                                                   <  5:259916/106‑1 (MQ=255)
                                            gcCAGGTGCAGGCCATGTCCCTGGAGTACGAGCCGCGCAAGGCCACGGGAGCAGAAAGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCtt                                                                                                      >  1:124921/1‑149 (MQ=255)
                                                   gCAGGCCATGTCCCTGGAGTACGAGCCGCGCAAGGCCACGGGAGCAGAAAGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGcc                                                                                               >  5:258588/1‑149 (MQ=255)
                                                    cAGGCCATGTCCCTGGAGTACGAGCCGCGCAAGGCCACGGGAGCAGAAAGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCg                                                                                              >  8:30572/1‑149 (MQ=255)
                                                      ggCCATGTCCCTGGAGTACGAGCCGCGCAAGGCCACGGGAGCAGAAAGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGc                                                                                            <  2:26428/149‑1 (MQ=255)
                                                             tCCCTGGAGTACGAGCCGCGCAAGGCCACGGGAGCAGAAAGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCAccc                                                                                     >  4:293596/1‑149 (MQ=255)
                                                              cccTGGAGTACGAGCCGCGCAAGGCCACGGGAGCAGAAAGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGt                                                                                                    <  2:79281/133‑1 (MQ=255)
                                                              cccTGGAGTACGAGCCGCGCAAGGCCACGGGAGCAGAAAGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGt                                                                                                    >  1:79281/1‑133 (MQ=255)
                                                               ccTGGAGTACGAGCCGCGCAAGGCCACGGGAGCAGAAAGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGCCACCGACACCGCTAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCgg                                                                                   >  7:305516/1‑149 (MQ=255)
                                                                     gTACGAGCCGCGCAAGGCCACGGGAGCAGAAAGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCGGACAAGc                                                                             <  7:30572/149‑1 (MQ=255)
                                                                          aGCCGCGCAAGGCCACGGGAGCAGAAAGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCa                                                                        <  2:124921/149‑1 (MQ=255)
                                                                           gCCGCGCAAGGCCACGGGAGCAGAAAGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACTGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCAg                                                                       >  8:219381/1‑149 (MQ=255)
                                                                                  aaGGCCACGGGAGCAGAAAGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCAGGTACCgg                                                                <  5:99088/149‑1 (MQ=255)
                                                                                    ggCCACGGGAGCAGAAAGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCAGGTACCGGTg                                                              >  6:132488/1‑149 (MQ=255)
                                                                                     gCCACGGGAGCAGAAAGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCAGGTACCGGTGc                                                             >  3:251567/1‑149 (MQ=255)
                                                                                           ggAGCAGAAAGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCAGGTACCGGTGCCAGc                                                         >  1:263139/1‑147 (MQ=255)
                                                                                            gatcAGAAAGCGTGCCCATGTCCTTGGCGGGAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCTTCGGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCAGGTACCGGTGCCAGCGAg                                                      <  8:93160/146‑1 (MQ=255)
                                                                                             aGCAGAAAGCGTGCCCATGACCTAGGCGGCAGCGTTGCGCTTGGTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCAGGTACCGGTGCCAGCGAgg                                                     <  8:273752/149‑1 (MQ=255)
                                                                                                aGAAAGCGTGCCCATGACCTACGCGGCAGCGTTGAGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGACGGCAGCACCCGGACCACCTGGCAGGTAACGGTGCCAGCGAGgcg                                                   >  6:234705/1‑148 (MQ=255)
                                                                                                 gAAAGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGTCGTATCGTCAGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCAGGTACCGGTGCCAGCGAGGCGCa                                                 <  8:273379/149‑1 (MQ=255)
                                                                                                  aaaGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCAGGTACCGGTGCCAGCGAGGCGCAg                                                >  6:66194/1‑149 (MQ=255)
                                                                                                  aaaGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGT‑GACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCAGGTACCGGTGCCAGCGAGGCGCAgg                                               >  4:154872/1‑149 (MQ=255)
                                                                                                  aaaGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGT‑GACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCAGGTACCGGTGCCAGCGAGGCGCAgg                                               >  4:149894/1‑149 (MQ=255)
                                                                                                   aaGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCAGGTACCGGTGCCAGCGAGGCGCAg                                                <  7:146349/148‑1 (MQ=255)
                                                                                                     gCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCAGGTACCGGTGCCAGCGAGGCGCAGGTg                                             >  8:117465/1‑149 (MQ=255)
                                                                                                       gTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCAGGTACCGGTGCCAGCGAGGCGCAGGTGCg                                           >  8:41194/1‑149 (MQ=255)
                                                                                                        tGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCAGGTACCGGTGCCAGCGAGGCGCAGGTGCGt                                          >  5:228639/1‑149 (MQ=255)
                                                                                                          cccATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCAGGTACCGGTGCCAGCGAGGCGCAGGTGCGTg                                         <  6:228639/148‑1 (MQ=255)
                                                                                                           ccATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACTGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCAGGTACCGGTGCCAGCGAGGCGCAGGTGCGTGaa                                       <  7:219381/149‑1 (MQ=255)
                                                                                                           ccATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCCCCCGGACAAGCTGGCAGGTACCGGTGCCAGCGAGGCGCAGGTGCGTg                                         >  7:18542/1‑147 (MQ=255)
                                                                                                           ccATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCAGGTACCGGTGCCAGCGAGGCGCAGGTGCGTGaa                                       <  5:132488/149‑1 (MQ=255)
                                                                                                           ccATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCAGGTACCGGTGCCAGCGAGGCGCAGGTGCGTGaa                                       <  2:263139/149‑1 (MQ=255)
                                                                                                                   tACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCAGGTACCGGTGCCAGCGAGGCGCAGGTGCGTGAAGCAACcgc                               >  3:209446/1‑149 (MQ=255)
                                                                                                                   tACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCAGGTACCGGTGCCAGCGAGGCGCAGGTGCGTGAAGCAACcgc                               >  5:41920/1‑149 (MQ=255)
                                                                                                                     cgcgGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCAGGTACCGGTGCCAGCGAGGCGCAGGTGCGTGAAGCAACcgcgc                             <  2:69764/149‑1 (MQ=255)
                                                                                                                      gcgGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCAGGTACCGGTGCCAGCGAGGCGCAGGTGCGTGAAGCAACcgcgcg                            >  6:8991/1‑149 (MQ=255)
                                                                                                                           aGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCAGGTACCGGTGCCAGCGAGGCGCAGGTGCGTGAAGCAACCGCGCGTAccc                       <  3:293596/149‑1 (MQ=255)
                                                                                                                             cGTTGCGCTTGCTGGCGGT‑GACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCAGGTACCGGTGCCAGCGAGGCGCAGGTGCGTGAAGCAACCGCGCGTACCCGTg                    <  3:149894/149‑1 (MQ=255)
                                                                                                                             cGTTGCGCTTGCTGGCGGT‑GACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCAGGTACCGGTGCCAGCGAGGCGCAGGTGCGTGAAGCAACCGCGCGTACCCGTg                    <  3:154872/149‑1 (MQ=255)
                                                                                                                               ttGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCAGGTACCGGTGCCAGCGAGGCGCAGGTGCGTg                                         <  6:39488/127‑1 (MQ=255)
                                                                                                                               ttGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCAGGTACCGGTGCCAGCGAGGCGCAGGTGCGTg                                         >  5:39488/1‑127 (MQ=255)
                                                                                                                               ttGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCAGGTACCGGTGCCAGCGAGGCGCAGGTGCGTGAAGCAACCGCGCGTACCCGTGa                   <  5:8991/149‑1 (MQ=255)
                                                                                                                                 gcgcTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCAGGTACCGGTGCCAGCGAGGCGCAGGTGCGTGAAGCAACCGCGCGTACCCGTGAg                  <  2:300148/148‑1 (MQ=255)
                                                                                                                                        cTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCAGGTACCGGTGCCAGCGAGGCGCAGGTGCGTGAAGCAACCGCGCGTACCCGTGAGTTGCAc            >  3:133051/1‑147 (MQ=255)
                                                                                                                                             ggTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCCGGTACCGGTGCCAGCGCGGCGCCGGTGCGTGCTGCAACCGCGCGTACCCGGGAGGTGCAGCCGGcc      >  8:117527/1‑148 (MQ=255)
                                                                                                                                             ggTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCAGGTACCGGTGCCGGCGAGGCGCAGGTGCGTGAAGCAACCGCGCGTACCCGTGAGTTGCACCAGGcc      >  5:140508/1‑148 (MQ=255)
                                                                                                                                                ggACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCAGGTACCGGTGCCAGCGAGGCGCAGGTGCGTGAAGCAACCGCGCGTACCCGTGAGTTGCACCAGGCCTATg  >  3:208281/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TGAGTTGCTGCTGGACTGGGGGCAGAAGCTGGGCATGAGCCGGCGCCAGGTGCAGGCCATGTCCCTGGAGTACGAGCCGCGCAAGGCCACGGGAGCAGAAAGCGTGCCCATGACCTACGCGGCAGCGTTGCGCTTGCTGGCGGTGGACGCAGACACCGACACCGATAAGGTCAAGCAGGCGTATCGTCGGCTTGTCAGCCGGCATCACCCGGACAAGCTGGCAGGTACCGGTGCCAGCGAGGCGCAGGTGCGTGAAGCAACCGCGCGTACCCGTGAGTTGCACCAGGCCTATG  >  NC_002947/494686‑494978

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: ‑