Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,539,020 G→C 6.2% L656V (CTG→GTG)  PP_4872 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,539,0200GC 6.2% 207.7 / 4.7 64L656V (CTG→GTG) PP_4872hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (35/25);  new base C (2/2);  total (37/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

ATAAGTATCGGTCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAG  >  NC_002947/5538877‑5539165
                                                                                                                                                |                                                                                                                                                 
aTAAGTATCGGTCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGc                                                                                                                                              >  7:17262/1‑149 (MQ=255)
     tATCGGTCTGCGGCTGGTGCAGCGGGATTAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGc                                                                                                                                         >  2:143069/1‑149 (MQ=255)
     tATCGGTCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGc                                                                                                                                         >  8:48847/1‑149 (MQ=255)
      aTCGGTCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGc                                                                                                                                         <  2:274598/148‑1 (MQ=255)
       tCGGTCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCt                                                                                                                                        <  1:136412/148‑1 (MQ=255)
           tCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCgcgtcg                                                                                                                                   >  4:298772/1‑149 (MQ=255)
                 gCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCtt                                                                                                                              >  6:107167/1‑148 (MQ=255)
                 gCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCGTCTGCACAAGGCGttctt                                                                                                                                        <  4:318045/138‑4 (MQ=255)
                 gCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCGTCTGCACAAGGCGttctt                                                                                                                                        >  3:318045/1‑135 (MQ=255)
                    ggTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGAc                                                                                                                           >  3:13427/1‑148 (MQ=255)
                        cAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGAtc                                                                                                                       >  8:116860/1‑148 (MQ=255)
                         aGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGTtctc                                                                                                                     >  1:194631/1‑149 (MQ=255)
                          gCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGAtctc                                                                                                                     >  4:327062/1‑148 (MQ=255)
                                   gTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCg                                                                                                           >  8:28384/1‑149 (MQ=255)
                                       gTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGctc                                                                                                        <  2:13461/148‑1 (MQ=255)
                                          cTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGt                                                                                                    <  5:235322/149‑1 (MQ=255)
                                           ttGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTg                                                                                                     >  2:66581/1‑147 (MQ=255)
                                             gcccggcccgATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGTCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTa                                                                                                   <  3:207959/147‑1 (MQ=255)
                                             gcccggcccgATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAg                                                                                                  <  4:188021/148‑1 (MQ=255)
                                                      gATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGctc                                                                                                        >  6:75019/1‑133 (MQ=255)
                                                       aTCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCAcc                                                                                            <  3:180707/144‑1 (MQ=255)
                                                       aTCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCAcc                                                                                            >  4:180707/1‑144 (MQ=255)
                                                             tcggctcggcCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGTTCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCgg                                                                                 <  2:194631/149‑1 (MQ=255)
                                                                 ctcggcCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCgg                                                                                 >  5:54517/1‑145 (MQ=255)
                                                                 ctcggcCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCgg                                                                                 <  6:54517/145‑1 (MQ=255)
                                                                 ctcggcCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTc                                                                             <  8:17262/149‑1 (MQ=255)
                                                                    ggcCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCtt                                                                           <  7:48847/148‑1 (MQ=255)
                                                                     gcCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGc                                                                         <  3:327062/149‑1 (MQ=255)
                                                                      cccATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCACTTGTAGGTAACCTTGCGCAGCGGGTTCTTGcc                                                                        >  6:148766/1‑149 (MQ=255)
                                                                         aTACGCTGCAGCGCTTCACGCAGCATCGGCGAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCt                                                                       >  1:30796/1‑147 (MQ=255)
                                                                         aTACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACTGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCt                                                                       >  1:267324/1‑148 (MQ=255)
                                                                          tACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAg                                                                    >  6:73817/1‑149 (MQ=255)
                                                                          tACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAg                                                                    >  7:62866/1‑149 (MQ=255)
                                                                                       ttCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACAtggccg                                                       >  4:312729/1‑149 (MQ=255)
                                                                                         cACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACAtggccgt                                                      <  4:13427/148‑1 (MQ=255)
                                                                                            gcagcaTCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAg                                                                                                  >  6:178657/1‑101 (MQ=255)
                                                                                            gcagcaTCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACAtggccgtggcc                                                  >  6:162427/1‑149 (MQ=255)
                                                                                               gcaTCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAg                                               <  3:312729/149‑1 (MQ=255)
                                                                                               gcaTCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAg                                               <  5:162427/149‑1 (MQ=255)
                                                                                                  tCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCTCGGCCACGGcc                                            <  2:30796/149‑1 (MQ=255)
                                                                                                     gCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCAtt                                         <  8:62866/149‑1 (MQ=255)
                                                                                                       cAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTgg                                       >  5:127588/1‑149 (MQ=255)
                                                                                                          cccTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGa                                    >  6:341480/1‑149 (MQ=255)
                                                                                                           ccTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGAc                                   <  5:341480/149‑1 (MQ=255)
                                                                                                           ccTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGAc                                   >  1:86921/1‑149 (MQ=255)
                                                                                                            cTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACg                                  >  4:70998/1‑149 (MQ=255)
                                                                                                             ttCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCAttt                                        <  1:224035/142‑2 (MQ=255)
                                                                                                             ttCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCATTCGCTCTTGTAGGTCACCTTGAGAAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTg                                        >  2:224035/1‑142 (MQ=255)
                                                                                                               cGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACgg                                 <  4:304492/147‑1 (MQ=255)
                                                                                                                gggTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGt                               <  5:73817/148‑1 (MQ=255)
                                                                                                                  gTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCa                                           <  1:312945/134‑1 (MQ=255)
                                                                                                                  gTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCa                                           >  2:312945/1‑134 (MQ=255)
                                                                                                                        gTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCt                      <  2:86921/149‑1 (MQ=255)
                                                                                                                             gcAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTAcc                 <  7:28384/149‑1 (MQ=255)
                                                                                                                             cgagcaaC‑GCCGTCTGCACAAGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACAtggccgtg                                                     <  3:268735/108‑1 (MQ=255)
                                                                                                                             cgagcaaC‑GCCGTCTGCACAAGGCCTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACAtggccgtg                                                     >  4:268735/6‑113 (MQ=255)
                                                                                                                              cAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCt                >  5:328950/1‑149 (MQ=255)
                                                                                                                                 gAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGt             >  6:12437/1‑149 (MQ=255)
                                                                                                                                  aaC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTc            >  4:202527/1‑149 (MQ=255)
                                                                                                                                  aaC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTc            >  5:80988/1‑149 (MQ=255)
                                                                                                                                        cTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGgcgc       >  4:189671/1‑149 (MQ=255)
                                                                                                                                        cTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGgcgc       >  2:151500/1‑149 (MQ=255)
                                                                                                                                        cTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTAGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGgcgc       >  7:166132/1‑149 (MQ=255)
                                                                                                                                              cAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAg  <  1:66581/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
ATAAGTATCGGTCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAG  >  NC_002947/5538877‑5539165

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