Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,075,747 G→C 15.1% A272G (GCC→GGC)  PP_5328 ← phosphate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,075,7470GC15.1% 183.4 / 24.6 73A272G (GCC→GGC) PP_5328phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (27/34);  new base C (6/5);  total (33/40)
Fisher's exact test for biased strand distribution p-value = 7.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-01

CAGCCCCCCCTTGGAGTCACCGATGATCAGCGAAATGCCGCCAACCAGCTGGGCCGTGGCAGTAATTTCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCGGTCATTTGCGGCAGTTCGATGCGGTTCTGCTCGGTGGTGACCTCTTCGGTCATCATGTTTTCCGTGCGGGTCAG  >  NC_002947/6075605‑6075886
                                                                                                                                                |                                                                                                                                               
cAGCCCCCCCTTGGAGTCACCGATGATCAGCGAAATGCCGCCAACCAGCTGGGCCGTGGCAGTAATTTCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGc                                                                                                                                          >  1:273466/1‑149 (MQ=255)
   cccccccTTCGAGTCACCGATGATCAGCGAAATGCCGCCACCCAGCTGGCCCGTGGCAGTAATTTCGGTGGCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACTGCCGt                                                                                                                                       <  6:220926/149‑1 (MQ=255)
    cgccccTTGGAGTCACCGATGATCAGCGAAATGCCGCCAACCAGCTGGGCCGTGGCAGTAATTTCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGtt                                                                                                                                      <  3:259171/147‑1 (MQ=255)
    ccccccTTGTAGTCACCGATGATCATCGAAATTCCGCCAACCAGCTGGGCCGTGGCAGTAATTTCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGCCACGGCCgac                                                                                                                                      <  6:21239/149‑3 (MQ=255)
    ccccccTTGGAGTCACCGATGATCAGCGAAATGCCGCCAACCAGCTGGGCCGTGGCAGTAATTTCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGtt                                                                                                                                      <  2:140045/149‑1 (MQ=255)
     cccccTTGGAGTCACCGATGATCAGCGAAATGCCGCCAACCAGCTGGGCCGTGGCAGTAATTTCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTAGGCCGTGGCACGGCCGtt                                                                                                                                      >  3:332738/1‑148 (MQ=255)
      ccccTTGGAGTCACCGATGATCAGCGAAATGCCGCCAACCAGCTGGGCCGTGGCAGTAATTTCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGa                                                                                                                                    <  5:18572/149‑1 (MQ=255)
      ccccGTGGAGGCACCGTTGATCAGCGAAATGCCGCCAACCCGCTGGGCCGTGGCAGTAATTTCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGCCACGGCCgacgt                                                                                                                                    <  1:376381/149‑5 (MQ=255)
        ccTTGGAGTCACCGATGATCAGCGAAATGCCGCCAACCAGCTGGGCCGTGGCAGTAATTTCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATc                                                                                                                                  <  3:172151/149‑1 (MQ=255)
        ccTTGGAGTCACCGATGATCAGCGAAATGCCGCCAACCAGCTGGGCCGTGGCAGTAATTTCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATc                                                                                                                                  <  1:64135/149‑1 (MQ=255)
        ccTTGGAGTCACCGATGATCAGCGAAATGCCGCCAACCAGCTGGGCCGTGGCAGTAATTTCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATc                                                                                                                                  <  3:166449/149‑1 (MQ=255)
         cTTGGAGTCACCGATGATCAGCGAAATGCCGCCAACCAGCTGGGCCGTGGCAGTAATTTCGGCGTCGGCGCTATCAGACAACTTGTAT‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCGGTGCCACGGCCgacgtgtt                                                                                                                                 <  4:133507/149‑8 (MQ=255)
            ggAGTCACCGATGATCAGCGAAATGCCGCCAACCAGCTGGGCCGTGGCAGTAATTTCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACg                                                                                                                               <  2:252969/148‑1 (MQ=255)
             gAGTCACCGATGATCAGCGAAATGCCGCCAACCAGCTGGGCCGTGGCAGTAATTTCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTa                                                                                                                             <  3:173632/149‑1 (MQ=255)
                     gatgatCAGCGAAATGCCGCCAACCAGCTGGGCCGTGGCAGTAATTTCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGc                                                                                                                         >  7:335668/1‑145 (MQ=255)
                     gatgatCAGCGAAATGCCGCCAACCAGCTGGGCCGTGGCAGTAATTTCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGc                                                                                                                         <  8:335668/145‑1 (MQ=255)
                     gatgatCAGCGAAATGCCGCCAACCAGCTGGGCCGTGGCAGTAATTTCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGt                                                                                                                     >  8:65733/1‑149 (MQ=255)
                       tgatCAGCGAAATGCCGCCAACCAGCTGGGCCGTGGCAGTAATTTCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGTCCGTTGAGGCTCTTGTCGCGCAGGCTTAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGt                                                                                                                     <  7:250716/147‑1 (MQ=255)
                         atCAGCGAAATGCCGCCAACCAGCTGGGCCGTGGCAGTAATTTCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑g                                                                                                              <  3:10563/149‑1 (MQ=255)
                         atCAGCGAAATGCCGCCAACCAGCTGGGCCGTGGCAGTAATTTCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑g                                                                                                              <  1:333171/149‑1 (MQ=255)
                                     ccgccAACCAGCTGGGCCGTGGCAGTAATTTCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGAt                                                                                                  >  7:324748/1‑149 (MQ=255)
                                      cgccAACCAGCTGGGCCGTGGCAGTAATTTCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGt                                                                                                          >  3:166183/1‑140 (MQ=255)
                                      cgccAACCAGCTGGGCCGTGGCAGTAATTTCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGt                                                                                                          <  4:166183/140‑1 (MQ=255)
                                       gccAACCAGCTGGGCCGTGGCAGTAATTTCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATg                                                                                                 >  5:154381/1‑148 (MQ=255)
                                         cAACCAGCTGGGCCGTGGCAGTAATTTCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAg                                                                                              <  5:168784/149‑1 (MQ=255)
                                           aCCAGCTGGGCCGTGGCAGTAATTTCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGa                                                                                             <  1:269836/149‑1 (MQ=255)
                                                   ggCCGTGGCAGTAATTTCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCa                                                                                    <  8:344166/149‑1 (MQ=255)
                                                     ccGTGGCAGTAATTTCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGAGCAGGCGGAACACGTCGGCCGTGTCACGGCCGTTGATCACGTACAGCCAGGGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCAcc                                                                                  <  6:289215/149‑1 (MQ=255)
                                                          gCAGTAATTTCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGtctc                                                                             <  2:28063/149‑1 (MQ=255)
                                                           cAGTAATTTCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCg                                                                            <  4:125621/149‑1 (MQ=255)
                                                            aGTAATTTCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGt                                                                                                      <  7:241434/122‑1 (MQ=255)
                                                            aGTAATTTCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGt                                                                                                      >  8:241434/1‑122 (MQ=255)
                                                            aGTAATTTCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCgg                                                                           >  5:237192/1‑149 (MQ=255)
                                                            aGTAATTTCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCgg                                                                           >  3:179734/1‑149 (MQ=255)
                                                            aGTAATTTCGGAGTCGGCGCTATCAGACAACTAGTAA‑‑CGGCCGTTGAGGCTCATGTCGCGCAGGCTGAACACGTACGCCGTGGCACGGCCGTTGATCACCTACAGCCAGTGCT‑‑‑CGCGTGGGTCGATGAACATATTCTTCACCGtcaagg                                                                           >  7:288758/1‑145 (MQ=255)
                                                             gTAATTTCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCCCGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCGGt                                                                          >  8:2588/1‑149 (MQ=255)
                                                              tAATTTCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGata                                                                                           >  6:235047/1‑131 (MQ=255)
                                                                 tttCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCGGTCAtt                                                                      <  2:351027/149‑1 (MQ=255)
                                                                   tCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCGGTCAttt                                                                     <  7:65733/148‑1 (MQ=255)
                                                                     ggTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCGGTCATTTGc                                                                   <  6:237192/148‑1 (MQ=255)
                                                                         tCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCGGTCATTTGCGGCAg                                                              <  2:273466/149‑1 (MQ=255)
                                                                                    aGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCGGTCATTTGCGGCAGTTCGATGCgg                                                    <  8:324748/148‑1 (MQ=255)
                                                                                    aGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCGGTCATTTGCGGCAGTTCGATGCGGt                                                   >  3:90203/1‑149 (MQ=255)
                                                                                    aGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCGGTCATTTGCGGCAGTTCGATGCGGt                                                   >  1:192226/1‑149 (MQ=255)
                                                                                     gACAACTTGTAACCCGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGGGGGTCGATGAAGATATTCTTCACCGTCTCGGTCATTTGCGGCAGTTCGATGCgg                                                    <  1:129385/149‑1 (MQ=255)
                                                                                     gACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCGGTCATTTGCGGCAGTTCGATGCGGtt                                                  >  5:180326/1‑149 (MQ=255)
                                                                                     gACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCGGTCATTTGCGGCAGTTCGATGCGGtt                                                  <  2:192226/149‑1 (MQ=255)
                                                                                     gACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCGGTCATTTGCGGCAGTTCGATGCGGt                                                   >  1:137767/1‑148 (MQ=255)
                                                                                        aaCTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCGGTCATTTGCGGCAGTTCGATGCgg                                                    <  3:339707/144‑1 (MQ=255)
                                                                                        aaCTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCGGTCATTTGCGGCAGTTCGATGCgg                                                    >  4:339707/1‑144 (MQ=255)
                                                                                           ttGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCGGTCATTTGCGGCAGTTCGATGCGGTTCTGCTc                                            >  6:208726/1‑149 (MQ=255)
                                                                                             gTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCGGTCATTTGCGGCAGTTCGATGCGGTTCTGCTCg                                           <  6:154381/148‑1 (MQ=255)
                                                                                              tAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCGGTCATTTGCGGCAGTTCGATGCGGTTCTGCTCgg                                          <  4:90203/148‑1 (MQ=255)
                                                                                                    ggCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCGGTCATTTGCGGCAGTTCGATGCGGTTCTGCTCGgtggtg                                     <  4:179734/149‑1 (MQ=255)
                                                                                                             ggCTCTTGTCGCGCAGGCTGAACTCGTCGTCCGTGGCTCGGACGTGGATCACGTACAGCCAGTGCT‑‑‑GGCGTGTGTCGCTGAAGATATTCTTGACCGTCTGGGTCATTTGCGGCAGTTCGATGCTGGTCTGCTCGGTGGTGACCTCTTCg                            <  8:288758/149‑1 (MQ=255)
                                                                                                             ggCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCTGAGGGAGTGTGTCTATGAAGATTATCTTCACAGTCTCGGTCTTTAGCGGGAGATCGATGCGTTTCATCTCTGTTGTGCGctct                               >  1:298855/1‑149 (MQ=255)
                                                                                                              gCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCGGTCATTTGCGGCAGTTCGATGCGGTTCTGCTCGGTGGTGACCTCTTCgg                           >  1:229902/1‑149 (MQ=255)
                                                                                                                tctTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCGGTCATTTGCGGCAGTTCGATGCGGTTCTGCTCGGTGGTGACCTCTTCGGt                          >  4:352682/1‑148 (MQ=255)
                                                                                                                    gTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCGGTCATTTGCGGCAGTTCGATGCGGTTCTGCTCGGTGGTGACCTCTTCGGtcatca                     >  6:155904/1‑149 (MQ=255)
                                                                                                                    gTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCGGTCATTTGCGGCAGTTCGATGCGGTTCTGCTCGGTGGTGACCTCTTCGGtcatca                     >  7:22933/1‑149 (MQ=255)
                                                                                                                        cgcAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCGGTCATTTGCGGCAGTTCGATGCGGTTCTGCTCGGTGGTGACCTCTTCGGTCATCATGtt                 <  7:2588/149‑1 (MQ=255)
                                                                                                                               tGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCGGTCATTTGCGGCAGTTCGATGCGGTTCTGCTCGGTGGTGACCTCTTCGGTCATCATGTTTTCCGTg          >  8:194246/1‑149 (MQ=255)
                                                                                                                                  acacTTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCGGTCATTTGCTGCAGTTCGATGCGGTTCTGCTCGGTGGTGACCTCTTCGGTCATCATGTTTTCCGTGCgt       <  4:332738/149‑2 (MQ=255)
                                                                                                                                   cacGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCGGTCATTTGCGGCAGTTCGATGCGGTTCTGCTCGGTGGTGACCTCTTCGGTCATCATGTTTTCCGTGCggg      >  6:76932/1‑149 (MQ=255)
                                                                                                                                    tcaacGGCCGTGCCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCGGTCATTTGCGGCAGTTCGATGCGGTTCTGCTCGGTGGTGACCTCTTCGGtcatc                      <  2:40592/128‑1 (MQ=255)
                                                                                                                                    tcaacGGCCGTGCCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCGGTCATTTGCGGCAGTTCGATGCGGTTCTGCTCGGTGGTGACCTCTTCGGtcatc                      >  1:40592/5‑132 (MQ=255)
                                                                                                                                    tcaacGGCCGTGCCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCGGTCATTTGCGGCAGTTCGATGCGGTTCTGCTCGGTGGTGACCTCTTCGGTCATCATGTTTTCCGTGCGGGt     >  8:36237/5‑149 (MQ=255)
                                                                                                                                     caacGGCCGTGCCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCGGTCATTTGCGGCAGTTCGATGCGGTTCTGCTCGGTGGTGACCTCTTCGGTCATCATGTTTTCCGTGCGGGt     >  4:220670/4‑148 (MQ=255)
                                                                                                                                      gTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCGGTCATTTGCGGCAGTTCGATGCGGTTCTGCTCGGTGGTGACCTCTTCGGTCATCATGTTTTCCGTGCGGGt     >  3:174144/1‑147 (MQ=255)
                                                                                                                                       acGGCCGTGCCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCGGTCATTTGCGGCAGTTCGATGCGGTTCTGCTCGGTGGTGACCTCTTCGGTCATCATGTTTTCCGTGCGGGTCAg  >  7:51332/2‑149 (MQ=255)
                                                                                                                                       acGGCCGTGCCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCGGTCATTTGCGGCAGGTCGATGCGGTTCTGCTCGGTGGTGACCTCTTCGGTCATCATGTTTTGCGTGCGGGTCAg  >  5:343025/2‑149 (MQ=255)
                                                                                                                                           ccGTGCCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCGGTCATTTGCGGCAGTTCGATGCGGTTCTGCTCGGTGGTGACCTCTTCGGTCATCATGTTTTCCGTGCgg       <  5:236655/140‑1 (MQ=255)
                                                                                                                                           ccGTGCCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCGGTCATTTGCGGCAGTTCGATGCGGTTCTGCTCGGTGGTGACCTCTTCGGTCATCATGTTTTCCGTGCgg       >  6:236655/1‑140 (MQ=255)
                                                                                                                                                |                                                                                                                                               
CAGCCCCCCCTTGGAGTCACCGATGATCAGCGAAATGCCGCCAACCAGCTGGGCCGTGGCAGTAATTTCGGTGTCGGCGCTATCAGACAACTTGTAA‑‑CGGCCGTTGAGGCTCTTGTCGCGCAGGCTGAACACGTCGGCCGTGGCACGGCCGTTGATCACGTACAGCCAGTGCT‑‑‑GGCGTGGGTCGATGAAGATATTCTTCACCGTCTCGGTCATTTGCGGCAGTTCGATGCGGTTCTGCTCGGTGGTGACCTCTTCGGTCATCATGTTTTCCGTGCGGGTCAG  >  NC_002947/6075605‑6075886

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: