Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,407,050 G→T 7.3% R123S (CGC→AGC)  PP_4751 ← amino acid ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,407,0500GT 7.3% 168.9 / 3.8 55R123S (CGC→AGC) PP_4751amino acid ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (28/23);  new base T (2/2);  total (30/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.21e-01

TACGTTGCTGCTGGCCACCAGAGCAGCGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGACCCGCCCCTTGGTCGGCGGCTC  >  NC_002947/5406903‑5407197
                                                                                                                                                    |                                                                                                                                                   
tACGTTGCTGCTGGCCACCAGAGCAGCGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTgcgc                                                                                                                                                    <  4:204908/149‑1 (MQ=255)
  cGTTGCTGCTGGCCACCAGAGCAGCGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCgca                                                                                                                                                   <  8:128132/148‑1 (MQ=255)
   gTTGCTGCTGGCCACCAGAGCAGCGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCgcag                                                                                                                                                  >  3:245707/1‑148 (MQ=255)
   gTTGCTGCTGGCCACCAGAGCAGCGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCgcag                                                                                                                                                  >  2:175325/1‑148 (MQ=255)
        tgctgGCCACCAGAGCAGCGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGt                                                                                                                                            >  4:237573/1‑149 (MQ=255)
          ctgGCCACCAGAGCAGCGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGt                                                                                                                                            >  4:229250/1‑147 (MQ=255)
           tgGCCACCAGAGCAGCGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAg                                                                                                                                         >  4:896/1‑149 (MQ=255)
             gacaccaGAGCAGCGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACTCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAgg                                                                                                                                        <  7:226524/146‑1 (MQ=255)
              ccaccaGAGCAGCGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTg                                                                                                                                      <  3:896/149‑1 (MQ=255)
              ccaccaGAGCAGCGCGCCGGATAATGCTGGGCCTTGCCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGGTGCTGCGCAACGTCAgcct                                                                                                                                      <  8:156982/149‑4 (MQ=255)
               caccaGAGCAGCGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGGTGCTGCGCAACGTCAg                                                                                                                                         <  5:49993/145‑1 (MQ=255)
               caccaGAGCAGCGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGGTGCTGCGCAACGTCAg                                                                                                                                         >  6:49993/1‑145 (MQ=255)
                 ccaGAGCAGCGCGCCGGATAATGCTGGGGCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCgg                                                                                                                                   <  4:304058/149‑1 (MQ=255)
                 ccaGAGCAGCGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCgg                                                                                                                                   <  4:245707/149‑1 (MQ=255)
                 ccaGAGCAGCGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCAAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCgg                                                                                                                                   >  5:243554/1‑149 (MQ=255)
                    gagCAGCGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCAcc                                                                                                                                              >  3:102416/1‑135 (MQ=255)
                    gagCAGCGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCAcc                                                                                                                                              <  4:102416/135‑1 (MQ=255)
                     agcagcGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAAc                                                                                                                               <  2:283423/149‑1 (MQ=255)
                             gcCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGaa                                                                                                                       <  3:237573/149‑1 (MQ=255)
                              ccGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAAc                                                                                                                      <  8:148713/149‑1 (MQ=255)
                                     aTGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTgg                                                                                                                <  2:143344/148‑1 (MQ=255)
                                     aTGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTgg                                                                                                                >  6:128438/1‑148 (MQ=255)
                                        cTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAg                                                                                                                             <  6:57633/132‑1 (MQ=255)
                                        cTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAg                                                                                                                             >  5:57633/1‑132 (MQ=255)
                                        cTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAaca                                                                                                            <  4:144473/149‑1 (MQ=255)
                                         tGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGa                                                                                                                                 >  4:161606/1‑127 (MQ=255)
                                         tGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGa                                                                                                                                 <  3:161606/127‑1 (MQ=255)
                                           ggCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAacac                                                                                                           <  2:35439/147‑1 (MQ=255)
                                              cTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCGGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATAc                                                                                                      <  4:261538/149‑1 (MQ=255)
                                                    agccagccCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATgc                                                                                                >  7:305504/1‑149 (MQ=255)
                                                    agccagccCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCAACGTCAGGTGCGGGAACAGGGTGAACGACTGGAACACCATACCAATgc                                                                                                >  4:268453/1‑149 (MQ=255)
                                                     gccagccCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCATGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATgc                                                                                                >  1:286079/1‑149 (MQ=255)
                                                     gccagccCAACGCGTACGATCAGTAGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGATGCGACGCACCGTCAGGTGCGGGAACAGCTTGAACGACAGGAACACCATAcccaaacg                                                                                               >  2:339435/1‑143 (MQ=255)
                                                             aaCGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGgcgc                                                                                        >  7:58993/1‑148 (MQ=255)
                                                                   ttCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACtt                                                                                   >  3:333009/1‑147 (MQ=255)
                                                                      gAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTcag                                                                                >  5:169609/1‑147 (MQ=255)
                                                                           gTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCt                                                                         <  6:169609/149‑1 (MQ=255)
                                                                              gCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTa                                                                      >  2:261976/1‑149 (MQ=255)
                                                                              gCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTa                                                                      >  5:365634/1‑149 (MQ=255)
                                                                                aaGGAGCGGGCGTCTGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACg                                                                    >  7:315998/1‑149 (MQ=255)
                                                                                               gCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGc                                                      >  2:219020/1‑148 (MQ=255)
                                                                                                 ttcttTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAg                                                   >  4:166263/1‑149 (MQ=255)
                                                                                                    tttGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCtt                                                >  2:10072/1‑149 (MQ=255)
                                                                                                      tGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCt                                                 <  8:305504/146‑1 (MQ=255)
                                                                                                      tGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTc                                               >  1:278002/1‑148 (MQ=255)
                                                                                                        gcgc‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACAATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGAt                                            <  6:243554/149‑1 (MQ=255)
                                                                                                                   ccAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGcca                                  <  6:365634/149‑1 (MQ=255)
                                                                                                                       ggTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCTAGCTTCGATGTGCACGccatcca                              <  1:175325/149‑1 (MQ=255)
                                                                                                                         tGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTTAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATg                            >  6:188401/1‑149 (MQ=255)
                                                                                                                            ggATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCt                          >  6:297133/1‑148 (MQ=255)
                                                                                                                             gttttGCGCCAGGCTGAAGTTGCGCAGCACCGTCATGTGCGGGTACAGGTTGAACGAATGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGa                        <  3:268453/147‑1 (MQ=255)
                                                                                                                             gATTTGCGCCAGTCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGTTTGAACGACTGGAACACCATACGAATGCGCCGGGGCAACTTCAGCAGTTCTGTACGGGGCAACGCATTGCCAGCTTCGATGTGCACGCCATCCATGCTGa                        <  1:339435/149‑1 (MQ=255)
                                                                                                                                   cgcCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGACCCGcc                  >  1:150282/1‑149 (MQ=255)
                                                                                                                                    gcCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGACCCGccc                 <  2:286079/149‑1 (MQ=255)
                                                                                                                                        ggCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGACCCGCCCCTTg             <  5:128438/149‑1 (MQ=255)
                                                                                                                                        ggCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGACCCGCCCCTTg             >  4:189445/1‑149 (MQ=255)
                                                                                                                                         gCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAATGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGACCCGCCCCTTgg            >  1:40671/1‑149 (MQ=255)
                                                                                                                                          cTGACGGTGCTGCGCAACGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGACCCGCCCCTTGGt           >  4:196569/1‑149 (MQ=255)
                                                                                                                                            gACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGACCCGCCCCTTGGTcg         >  2:42795/1‑149 (MQ=255)
                                                                                                                                                   cgcAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGACCCGCCCCTTGGTCGGCGGCTc  <  5:297133/149‑1 (MQ=255)
                                                                                                                                                   cgcAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGACCCGCCCCTTGGTCGGCGGCTc  <  8:315998/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TACGTTGCTGCTGGCCACCAGAGCAGCGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGC‑TGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGACCCGCCCCTTGGTCGGCGGCTC  >  NC_002947/5406903‑5407197

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