Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 806,327 T→C 10.7% K335E (AAG→GAG)  PP_0693 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947806,3270TC10.7% 148.7 / 7.4 56K335E (AAG→GAG) PP_0693hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (25/25);  new base C (3/3);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGACCAGCTCTTCCTT  >  NC_002947/806180‑806474
                                                                                                                                                   |                                                                                                                                                           
ccAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCtt                                                                                                                                                             <  1:139278/148‑1 (MQ=255)
ccAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCtt                                                                                                                                                             >  3:40996/1‑148 (MQ=255)
  aGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAg                                                                                                                                                          <  1:127045/149‑1 (MQ=255)
  aGGTCGCGTACGCGCTCGTCGACGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAg                                                                                                                                                          <  2:242558/149‑1 (MQ=255)
      cgcgTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAc‑‑‑‑c                                                                                                                                               <  5:153863/149‑1 (MQ=255)
       gcgTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAc‑‑‑‑c                                                                                                                                               <  4:40996/148‑1 (MQ=255)
       gcgTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAc‑‑‑‑c                                                                                                                                               >  7:221325/1‑148 (MQ=255)
         gTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCAGCTCCAGTACCTTGCCGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGa                                                                                                                                             >  8:90576/1‑148 (MQ=255)
            cgcgCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTCAAG‑‑‑CAC‑‑‑‑Cgctttg                                                                                                                                         <  4:49326/149‑6 (MQ=255)
               gcTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCg                                                                                                                                     >  5:268609/1‑149 (MQ=255)
                   gtcgAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGcc                                                                                                                                   <  3:103815/147‑1 (MQ=255)
                      gAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGc                                                                                                                               <  6:268609/148‑1 (MQ=255)
                          tCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCCCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCtt                                                                                                                                                             >  1:175203/1‑122 (MQ=255)
                          tCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCCCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCtt                                                                                                                                                             <  2:175203/122‑1 (MQ=255)
                           cTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAgcgc                                                                                                                         >  2:234577/1‑149 (MQ=255)
                              cGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCt                                                                                                                        <  1:234577/147‑1 (MQ=255)
                                           ccGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAAtt                                                                                                          >  4:32252/1‑148 (MQ=255)
                                             gTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTcc                                                                                                       >  6:69336/1‑148 (MQ=255)
                                                  tGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTg                                                                                                                       >  6:267136/1‑128 (MQ=255)
                                                  tGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTg                                                                                                                       <  5:267136/128‑1 (MQ=255)
                                                    ggCAAGCGCGCCGGCCAGCGCGTCGTTGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCa                                                                                                <  4:159522/149‑1 (MQ=255)
                                                        aGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGc                                                                                                                               >  4:252083/1‑114 (MQ=255)
                                                        aGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGgcg                                                                                            <  6:51464/149‑1 (MQ=255)
                                                           gcgcCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCAt                                                                                         <  1:5985/149‑1 (MQ=255)
                                                                  ccAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTCGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTcc                                                                                  >  3:211085/1‑149 (MQ=255)
                                                                     gCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTgg                                                                               <  8:221325/149‑1 (MQ=255)
                                                                          tcgtGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGGTAAACGGGAGGATGAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTcc                                                                                  >  1:86821/1‑149 (MQ=255)
                                                                          tcgtGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTgcc                                                                          <  1:31468/149‑1 (MQ=255)
                                                                             tGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTgccgc                                                                        >  4:67684/1‑148 (MQ=255)
                                                                             tGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTgccgc                                                                        >  4:26310/1‑148 (MQ=255)
                                                                               ggCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCa                                                                      <  3:32252/148‑1 (MQ=255)
                                                                                    cTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTgg                                                                >  5:42660/1‑149 (MQ=255)
                                                                                     tttGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGt                                                               >  7:25608/1‑149 (MQ=255)
                                                                                      ttGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTc                                                              >  1:134998/1‑149 (MQ=255)
                                                                                                gcagTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAgg                                                    <  3:67684/149‑1 (MQ=255)
                                                                                                      gctgctCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTcc                                              >  1:34052/1‑149 (MQ=255)
                                                                                                      gctgctCCAGTACCTTGACGTCATCGACCAGCGCAAAACGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTAGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTcc                                              >  1:47374/1‑149 (MQ=255)
                                                                                                         gctCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGc                                           >  2:3960/1‑149 (MQ=255)
                                                                                                          ctCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCt                                             <  6:195766/146‑1 (MQ=255)
                                                                                                          ctCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCt                                             >  5:195766/1‑146 (MQ=255)
                                                                                                                  ccTTGACGGCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACgg                                  <  7:116739/149‑1 (MQ=255)
                                                                                                                   cTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACgg                                  >  2:148075/1‑148 (MQ=255)
                                                                                                                      gACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAg                               <  1:148075/148‑1 (MQ=255)
                                                                                                                            aTCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCg                         >  4:230765/1‑148 (MQ=255)
                                                                                                                             tCGACCAGCGCAAAGCGGTGCTCGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGa                       <  4:211085/149‑1 (MQ=255)
                                                                                                                                 ccAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTgg                   <  1:3960/149‑1 (MQ=255)
                                                                                                                                       cAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTc                                               >  8:81489/1‑115 (MQ=255)
                                                                                                                                       cAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTc                                               <  7:81489/115‑1 (MQ=255)
                                                                                                                                        aaaGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGAcc             >  2:251865/1‑148 (MQ=255)
                                                                                                                                        aaaGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGAcc             <  1:251865/148‑1 (MQ=255)
                                                                                                                                        aaaGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGAc              >  1:64876/1‑147 (MQ=255)
                                                                                                                                         tgtcGGTGCTCAAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGACCAg           >  2:234879/4‑149 (MQ=255)
                                                                                                                                          aGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGACCAGc          <  2:134998/149‑1 (MQ=255)
                                                                                                                                           tcGGTGCTCAAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGACCAGct         >  5:102781/2‑149 (MQ=255)
                                                                                                                                              gTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGACCAGCTCtt      >  5:171612/1‑149 (MQ=255)
                                                                                                                                              gTGCTCAAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGACCAGCTCtt      <  6:102781/149‑1 (MQ=255)
                                                                                                                                               tGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGACCAGCTCTTc     <  8:25608/149‑1 (MQ=255)
                                                                                                                                                 cTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCGTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGACCAGCTCTTCCt   >  5:194173/1‑149 (MQ=255)
                                                                                                                                                 cTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGACCAGCTCTTcc    <  3:230765/148‑1 (MQ=255)
                                                                                                                                                 cTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGACCAGCTCTTcc    <  2:34052/148‑1 (MQ=255)
                                                                                                                                                  ttGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGACCAGCTCTTCCtt  >  7:247248/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                           
CCAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGACCAGCTCTTCCTT  >  NC_002947/806180‑806474

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