Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,843,401 T→G 5.7% I26L (ATC→CTC)  PP_3392 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,843,4010TG 5.7% 208.5 / 2.3 67I26L (ATC→CTC) PP_3392hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (30/33);  new base G (2/2);  total (32/35)
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

CCACTGTCAACGCAAAGCCGGGGTCGGCATGCAAGACAGCACCCTCCTGCACGTTGCTGTGCTGGCCGATACGGATCGCTTCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTGGTAGATCGCCATCTCAAGGCCTCCTCAAGCCGCGTCGTCTGTTTTCAGCACGCCCGCCGCACACAGCTGGGCAATGT  >  NC_002947/3843254‑3843541
                                                                                                                                                   |                                                                                                                                            
ccACTGTCAACGCAAAGCCGGGGTCGGCATGCAAGACAGCACCCTCCTGCACGTTGCTGTGCTGGCCGATACGGATCGCTTCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCAAGGGTAACGTCGCCGAt                                                                                                                                              >  1:51498/1‑148 (MQ=255)
 cACTGTCAACGCAAAGCCGGGGTCGGCATGCAAGACAGCACCCTCCTGCACGTTGCTGTGCTGGCCGATACGGATCGCTTCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGa                                                                                                                                            >  6:84921/1‑149 (MQ=255)
      tCAACGCAAAGCCGGGGTCGGCATGCAAGACAGCACCCTCCTGCACGTTGCTGTGCTGGCCGATACGGATCGCTTCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTg                                                                                                                                       <  2:202989/149‑1 (MQ=255)
        aaCGCAAAGCCGGGGTCGGCATGCAAGACAGCACCCTCCTGCACGTTGCTGTGCTGGCCGATACGGATCGCTTCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATg                                                                                                                                             >  4:148997/1‑141 (MQ=255)
        aaCGCAAAGCCGGGGTCGGCATGCAAGACAGCACCCTCCTGCACGTTGCTGTGCTGGCCGATACGGATCGCTTCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATg                                                                                                                                             <  3:148997/141‑1 (MQ=255)
         aCGCAAAGCCGGGGTCGGCATGCAAGACAGCACCCTCCTGCACGTTGCTGTGCTGGCCGATACGGATCGCTTCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCa                                                                                                                                    >  7:56599/1‑149 (MQ=255)
             aaaGCCGGGGTCGGCATGCAAGACAGCACCCTCCTGCACGTTGCTGTGCTGGCCGATACGGATCGCTTCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCt                                                                                                                                >  1:8301/1‑149 (MQ=255)
             aaaGCCGGGGTCGGCATGCAAGACAGCACCCTCCTGCACGTTGCTGTGCTGGCCGATACGGATCGCTTCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCt                                                                                                                                >  2:185569/1‑149 (MQ=255)
                   ggggTCGGCATGCAAGACAGCACCCTCCTGCACGTTGCTGTGCTGGCCGATACGGATCGCTTCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGAGGATGCCAGCGTCATCGGCGa                                                                                                                          <  3:211204/149‑1 (MQ=255)
                     ggTCGGCATGCAAGACAGCACCCTCCTGCACGTTGCTGTGCTGGCCGATACGGATCGCTTCGTTGTCGCCCCGTAGCACGGACTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGAt                                                                                                                         >  8:317709/1‑148 (MQ=255)
                        cGGCATGCAAGACAGCACCCTCCTGCACGTTGCTGTGCTGGCCGATACGGATCGCTTCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAAc                                                                                                                     <  4:133579/149‑1 (MQ=255)
                          gCATGCAAGACAGCACCCTCCTGCACGTTGCTGTGCTGGCCGATACGGATCGCTTCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGt                                                                                                                   >  8:295497/1‑149 (MQ=255)
                              gCAAGACAGCACCCTCCTGCACGTTGCTGTGCTGGCCGATACGGATCGCTTCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTcc                                                                                                               >  5:64898/1‑149 (MQ=255)
                                aaGACAGCACCCTCCTGCACGTTGCTGTGCTGGCCGATACGGATCGCTTCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCgg                                                                                                             <  1:286075/149‑1 (MQ=255)
                                aaGACAGCACCCTCCTGCACGTTGCTGTGCTGGCCGATACGGATCGCTTCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGCGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCgg                                                                                                             >  7:96560/1‑149 (MQ=255)
                                       cACCCTCCTGCACGTTGCTGTGCTGGCCGATACGGATCGCTTCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAgg                                                                                                      >  2:370434/1‑149 (MQ=255)
                                          cctcctGCACGTTGCTGTGCTGGCCGATACGGATCGCTTCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTg                                                                                                   >  4:401879/1‑149 (MQ=255)
                                          cctcctGCACGTTGCTGTGCTGGCCGATACGGATCGCTTCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTg                                                                                                   >  4:25509/1‑149 (MQ=255)
                                              ctGCACGTTGCTGTGCTGGCCGATACGGATCGCTTCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTggg                                                                                                 <  3:427834/147‑1 (MQ=255)
                                              ctGCACGTTGCTGTGCTGGCCGATACGGATCGCTTCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGt                                                                                                <  2:264880/148‑1 (MQ=255)
                                               tGCACGTTGCTGTGCTGGCCGATACGGATCGCTACGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGGTCCAGGGTAACGTCGCCGAGGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGggtgggtg                                                                                               >  7:278954/1‑148 (MQ=255)
                                                gCACGTTGCTGTGCTGGCCGATACGGATCGCTTCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACAATAGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGc                                                                                              >  1:146300/1‑148 (MQ=255)
                                                  aCGTTGCTGTGCTGGCCGATACGGATCGCTTCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGATTCCGGGTGCAGTTTGGGTGAAAg                                                                                           >  2:427127/1‑149 (MQ=255)
                                                     ttGCTGTGCTGGCCGATACGGATCGCTTCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGtt                                                                                        <  6:266653/149‑1 (MQ=255)
                                                     ttGCTGTGCTGGCCGATACGGATCGCTTCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGtt                                                                                        <  5:354194/149‑1 (MQ=255)
                                                        ctgtgctgGCCGATCCGGATCGCTTCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCTTCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTc                                                                                     <  8:96560/149‑1 (MQ=255)
                                                          gtgctgGCCGATACGGATCGCTTCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCa                                                                                    >  3:172132/1‑148 (MQ=255)
                                                            gctgGCCGATACGGATCGCTTCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCatat                                                                                 <  7:317709/149‑1 (MQ=255)
                                                            gctgGCCGATACGGATCGCTTCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCatat                                                                                 >  8:138228/1‑149 (MQ=255)
                                                             ctgGCCGATACGGATCGCTTCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATAtt                                                                                <  5:84921/149‑1 (MQ=255)
                                                              tgGCCGATACGGATGGGTTCGTTGTCGCCCCGTAGCACGGCGTGCGGCCATAGGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTg                                                                               <  8:278954/149‑1 (MQ=255)
                                                              tgGCCGATACGGATCGCTTCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTg                                                                               >  8:56952/1‑149 (MQ=255)
                                                              tgGCCGATACGGATCGCTTCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTg                                                                               >  8:50139/1‑149 (MQ=255)
                                                                gCCGATACGGATCGCTTCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTg                                                                               <  3:401879/147‑1 (MQ=255)
                                                                gCCGATACGGATCGCTTCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTGGt                                                                             <  3:25509/149‑1 (MQ=255)
                                                                  cGATACGGATCGCTTCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTGGTAg                                                                           <  1:370434/149‑1 (MQ=255)
                                                                              cTTCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTGGTAGATCGCCAtctc                                                                >  6:273582/1‑148 (MQ=255)
                                                                               ttCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTGGTAGATCGCCAtctc                                                                >  2:340210/1‑147 (MQ=255)
                                                                                 ctttGTCGCCCCGTAGCACTGCCTGCGTCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGATGGCATCCTCTGCGATGAACGTTTCCTGGTGCATGTTGGGTGCAAGGTTGTCATATTGGTAGATCGCCATCTCaa                                                              <  3:264332/145‑1 (MQ=255)
                                                                                 cGTTGTCGCCCCGTAGCACGGGCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTGGTAGATCGCCATCTCAAg                                                             <  1:427127/148‑1 (MQ=255)
                                                                                    tGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTGGTAGATCGCCATCTCAAGGcct                                                         >  1:362899/1‑149 (MQ=255)
                                                                                      tCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCTTATTGGTAGATCGCCATCTCAAGGcctcc                                                       <  1:34637/149‑1 (MQ=255)
                                                                                      tCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTGGTAGATCGCCATCTCAAGGcctcc                                                       >  2:328665/1‑149 (MQ=255)
                                                                                        gCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTGGTAGATCGCCATCTCAAGGcctcct                                                      >  5:395321/1‑148 (MQ=255)
                                                                                          cccGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTGGTAGATCGCCATCTCAAGGCCTCCTCa                                                    >  8:375030/1‑148 (MQ=255)
                                                                                           ccGTAGCACGGCCGGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTGGTAGATCGCCATCTCAAGGCCTCCTCAAg                                                  <  2:8301/149‑1 (MQ=255)
                                                                                            cGTAGCACGGCCTGCGGCCATACGCTGGCGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTGGTAGATCGCCATCTCAAGGCCTCCTCAAGc                                                 <  4:287761/149‑1 (MQ=255)
                                                                                             gTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTGGTAGATCGCCATCTCAAGGCCTCCTCAAGcc                                                <  1:328665/149‑1 (MQ=255)
                                                                                                gCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTGGTAGATCGCCATCTCAAGGCCTCCTCAAGCcgc                                              <  6:64898/148‑1 (MQ=255)
                                                                                                  aCGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTGGTAGATCGCCATCTCAAGGcctc                                                        >  5:259008/1‑136 (MQ=255)
                                                                                                  aCGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTGGTAGATCGCCATCTCAAGGcctc                                                        <  6:259008/136‑1 (MQ=255)
                                                                                                            gCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTGGTAGATCGCCATCTCAAGGCCTCCTCAAGCCGCGTCGTCTGTTTTc                                 <  7:295497/149‑1 (MQ=255)
                                                                                                              cATATGCTGTCGCGTTGTTCCAGGGTAACGTCGCCGATGACGCTTGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTGGTTGATCGCCATCTCAAGGCCTGCTGAAGCGGCGTCGTCTGTTTTCAg                               <  2:51498/149‑1 (MQ=255)
                                                                                                               aTACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTGGTAGATCGCCATCTCAAGGCCTCCTCAAGCCGCGTCGTCTGTTTTCAGc                              <  8:56599/149‑1 (MQ=255)
                                                                                                                   gCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTGGTAGATCGCCATCTCAAGGCCTCCTCAAGCCGCGTCGTCTGTTTTCAGCAcgc                          >  5:221468/1‑149 (MQ=255)
                                                                                                                       aCGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTGGTAGATCGCCATCTCAAGGCCTCCTCAAGCCGCGTCGTCTGTTTTCAGCAcgccc                        <  5:273582/147‑1 (MQ=255)
                                                                                                                         gCCTTGTTCAAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTGGTAGATCGCCATCTCAAGGCCTCCTCAAGCCGCGTCGTCTGTTTTCAGCACGCccgccg                    >  6:388633/1‑149 (MQ=255)
                                                                                                                           cTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTGGTAGATCGCCATCTCAAGGCCTCCTCAAGCCGCGTCGTCTGTTTTCAGCACGCCCGCCGca                  >  7:179813/1‑149 (MQ=255)
                                                                                                                                  cAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTGGTAGATCGCCATCTCAAGGCCTCCTCAAGCCGCGTCGTCTGTTTTCAGCACGCCCGCCGCACACAGCt           <  6:221468/149‑1 (MQ=255)
                                                                                                                                    gggTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTGGTAGATCGCCATCTCAAGGCCTCCTCAAGCCGCGTCGTCTGTTTTCAGCACGCCCGCCGCACACAGCTgg         <  4:172132/149‑1 (MQ=255)
                                                                                                                                     ggTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTGGTAGATCGCCATCTCAAGGCCTCCTCAAGCCGCGTCGTCTGTTTTCAGCACGCCCGCCGCACACAGCTggg        <  7:138228/149‑1 (MQ=255)
                                                                                                                                       gttcatCGCCGAGGATGCCAGCGTCATCGGCGACGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTGGTAGATCGCCATCTCAAGGCCTCCTCAAGCCGCGTCGTCTGTTTTCAGCACGCCCGCCGCACACAGCTgg         <  3:252638/141‑1 (MQ=38)
                                                                                                                                       gttcatCGCCGAGGATGCCAGCGTCATCGGCGACGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTGGTAGATCGCCATCTCAAGGCCTCCTCAAGCCGCGTCGTCTGTTTTCAGCACGCCCGCCGCACACAGCTgg         >  4:252638/6‑146 (MQ=38)
                                                                                                                                         aCGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTGGTAGATCGCCATCTCAAGGCCTCCTCAAGCCGCGTCGTCTGTTTTCAGCAcgcccgc                      <  7:274351/131‑1 (MQ=255)
                                                                                                                                          cGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTGGTAGATCGCCATCTCAAGGCCTCCTCAAGCCGCGTCGTCTGTTTTCAGCACGCCCGCCGCACACAGCTGGGCAAtg   <  2:146300/149‑1 (MQ=255)
                                                                                                                                           gTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTGGTAGATCGCCATCTCAAGGCCTCCTCAAACCGCGTCGTCTGTTTTCAGCACGCCCGCCGCACACAGCTGGGCAAtgt  >  5:346732/1‑149 (MQ=255)
                                                                                                                                            tCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTGGTAGATCGCCATCTCAAGGCCTCCTCAAGCCGCGTCGTCTGTTTTCAGCACGCCCGCCGCACACAGCTGGGCAAtgt  <  7:55962/148‑1 (MQ=255)
                                                                                                                                            tCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTGGTAGATCGCCATCTCAAGGCCTCCTCAAGCCGCGTCGTCTGTTTTCAGCACGCCCGCCGCACACAGCTGGGCAAtgt  >  8:55962/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
CCACTGTCAACGCAAAGCCGGGGTCGGCATGCAAGACAGCACCCTCCTGCACGTTGCTGTGCTGGCCGATACGGATCGCTTCGTTGTCGCCCCGTAGCACGGCCTGCGGCCATACGCTGACGCCTTGTTCCAGGGTAACGTCGCCGATGACGCTGGCATCCTCGGCGATGAACGTTTCCGGGTGCAGGTTGGGTGCAAGGTTGTCATATTGGTAGATCGCCATCTCAAGGCCTCCTCAAGCCGCGTCGTCTGTTTTCAGCACGCCCGCCGCACACAGCTGGGCAATGT  >  NC_002947/3843254‑3843541

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