Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,793,708 T→A 6.3% intergenic (+17/+110) PP_2445 → / ← PP_2446 membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,793,7080TA 6.3% 178.2 / 2.5 63intergenic (+17/+110)PP_2445/PP_2446membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (27/32);  new base A (2/2);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

ATCGCCGAAGGCTTTGGCGCCCATGTACCCAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTAC  >  NC_002947/2793560‑2793853
                                                                                                                                                     |                                                                                                                                                  
aTCGCCGAAGGCTTTGGCGCCCATGTACCCAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCt                                                                                                                                                    >  4:115239/1‑149 (MQ=255)
 tCGCCGAAGGCTTTGGCGCCCATGTACCCAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTc                                                                                                                                                   >  5:5412/1‑149 (MQ=255)
  cGCCGAAGGCTTTGGCGCCCATGTACCCAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTc                                                                                                                                                   <  3:80877/148‑1 (MQ=255)
  cGCCGAAGGCTTTGGCGCCCATGTACCCAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTAAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTcg                                                                                                                                                  >  5:267386/1‑149 (MQ=255)
   gCCGAAGGCTTTGGCGCCCATGTACCCAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTcgc                                                                                                                                                 <  5:271641/149‑1 (MQ=255)
   gCCGAAGGCTTTGGCGCCCATGTACCCAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTcgc                                                                                                                                                 <  1:252961/149‑1 (MQ=255)
    ccGAAGGCTTTGGCGCCCATGTACCCAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTcgcg                                                                                                                                                >  6:113731/1‑149 (MQ=255)
    cagAAGGCTTTGGCGCCCATGTACCCAAGGGCTATGTGTATGCGTCAATGGCGTTCTCGACGGCTATCGAGCTCATCAACATACTGGCCCATCGAGAGCGTTTGAAGCGAGAGGAGCAAGAAGGTTGATTTAAAAAGGGGCACCTcgcg                                                                                                                                                >  4:65759/3‑149 (MQ=255)
           cTTTGGCGCCCATGTACCCAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCctct                                                                                                                                         >  2:258613/1‑149 (MQ=255)
            tttGGCGCCCATGTACCCAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCtt                                                                                                                                        >  2:123629/1‑149 (MQ=255)
               ggCGCCCATGTACCCAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGc                                                                                                                                      <  3:115239/148‑1 (MQ=255)
               ggCGCCCATGTACCCAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGc                                                                                                                                      <  2:109721/148‑1 (MQ=255)
                    ccATGTACCCAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGAc                                                                                                                                  >  5:146598/1‑147 (MQ=255)
                           cccAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCt                                                                                                                         <  7:208119/149‑1 (MQ=255)
                            ccAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTg                                                                                                                        >  8:231431/1‑149 (MQ=255)
                            ccAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTg                                                                                                                        >  6:76753/1‑149 (MQ=255)
                            ccAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTg                                                                                                                        >  7:51559/1‑149 (MQ=255)
                            ccAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTg                                                                                                                        >  8:231434/1‑149 (MQ=255)
                             cAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCGAGGCGCCCcttt                                                                                                                                         <  3:1610/131‑3 (MQ=255)
                             cAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCGAGGCGCCCcttt                                                                                                                                         <  3:5503/131‑3 (MQ=255)
                             cAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCGAGGCGCCCcttt                                                                                                                                         >  4:5503/1‑129 (MQ=255)
                             cAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCGAGGCGCCCcttt                                                                                                                                         >  4:1610/1‑129 (MQ=255)
                             cAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTg                                                                                                                        >  6:147366/1‑148 (MQ=255)
                               aGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGcac                                                                                                                     <  5:113731/149‑1 (MQ=255)
                                          tATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTcgc                                                                                                           <  8:51908/148‑1 (MQ=255)
                                               ggCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCa                                                                                                     >  3:170611/1‑149 (MQ=255)
                                                  cATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTg                                                                                                  <  1:123629/149‑1 (MQ=255)
                                                    tGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTgg                                                                                                 >  2:156945/1‑148 (MQ=255)
                                                    tGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGc                                                                                                <  6:90757/149‑1 (MQ=255)
                                                         ttCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAAcgc                                                                                           <  1:258613/149‑1 (MQ=255)
                                                             cGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTcgc                                                                                                           >  2:85739/1‑129 (MQ=255)
                                                             cGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCCTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTcgc                                                                                                           <  1:85739/129‑1 (MQ=255)
                                                              gACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTCGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGaa                                                                                        >  1:70249/1‑148 (MQ=255)
                                                              gACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGaa                                                                                        <  4:170611/147‑1 (MQ=255)
                                                               aCGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTTGCCAGTGCCACGCGAATAc                                                                                     <  6:5412/149‑1 (MQ=38)
                                                                cGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGTGGACAAGCCCGCTGCACAGGTATTCGCGTTGCCAGTGCCACGCGAATAcc                                                                                    <  1:87877/149‑1 (MQ=37)
                                                                   cTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCtgt                                                                                 <  5:147366/149‑1 (MQ=255)
                                                                    tATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCtgtg                                                                                >  1:51757/1‑149 (MQ=255)
                                                                              cTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGa                                                                                                                                   <  2:114588/88‑1 (MQ=255)
                                                                              cTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGa                                                                                                                                   >  1:114588/1‑88 (MQ=255)
                                                                              cTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCtt                                                                      <  2:51757/149‑1 (MQ=255)
                                                                                 aaCATCCTGACCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTgt                                                                    >  5:271970/1‑148 (MQ=255)
                                                                                     tCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTgtcagt                                                                >  3:90956/1‑148 (MQ=255)
                                                                                              gTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTgg                                                                                                 <  8:201778/106‑1 (MQ=255)
                                                                                              gTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTgg                                                                                                 >  7:201778/1‑106 (MQ=255)
                                                                                                 gAGCGCGTTTGAAGCGTGAGGCTGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCTGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCgg                                                   <  3:65759/149‑1 (MQ=255)
                                                                                                 gAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCg                                                    <  8:51559/148‑1 (MQ=255)
                                                                                                 gAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCg                                                    <  6:146598/148‑1 (MQ=255)
                                                                                                       gTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATc                                             <  6:271970/149‑1 (MQ=255)
                                                                                                            aaGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGc                                                                                                                       <  6:124881/70‑1 (MQ=255)
                                                                                                            aaGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGc                                                                                                                       >  5:124881/1‑70 (MQ=255)
                                                                                                            aaGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGTTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAg                                        <  8:101462/149‑1 (MQ=255)
                                                                                                            aaGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGa                                         <  4:90956/148‑1 (MQ=255)
                                                                                                                           tGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAAcgc                                                                                           <  7:76065/84‑1 (MQ=255)
                                                                                                                           tGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAAcgc                                                                                           >  8:76065/1‑84 (MQ=255)
                                                                                                                            gAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCgg                        >  2:223019/1‑149 (MQ=255)
                                                                                                                               gtttGATTTAAAATGTGGCGCCTCGGGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGTCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCAt                     <  2:249284/147‑1 (MQ=38)
                                                                                                                               ggTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCAt                     <  7:231431/149‑1 (MQ=255)
                                                                                                                               ggTTGATGAAAAACGGGGCGCGTCGCGCCCGTCTTGCGGACAAGCCCGCTGCACCGGTATTCGCGGGGCACTGCCACCGCGAATACCTGTGCTGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCAt                     <  7:231434/149‑1 (MQ=255)
                                                                                                                                     tGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATgccgcc               >  8:34492/1‑149 (MQ=255)
                                                                                                                                       tAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGCATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCgg             >  6:20001/1‑149 (MQ=255)
                                                                                                                                       tAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCgg             <  7:34492/149‑1 (MQ=255)
                                                                                                                                       tAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCgg             >  2:11226/1‑149 (MQ=255)
                                                                                                                                                   ccTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCTACCCCGGCATGCCGCCGGGGCCGTCTTAc  <  1:156945/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
ATCGCCGAAGGCTTTGGCGCCCATGTACCCAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTAC  >  NC_002947/2793560‑2793853

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: