Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 77,503 A→T 9.1% V14E (GTG→GAG)  rsmB ← ribosomal RNA small subunit methyltransferase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294777,5030AT 9.1% 199.9 / 9.1 66V14E (GTG→GAG) rsmBribosomal RNA small subunit methyltransferase B
Reads supporting (aligned to +/- strand):  ref base A (36/24);  new base T (3/3);  total (39/27)
Fisher's exact test for biased strand distribution p-value = 6.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.39e-01

CAGCAGCTGTGCGGCCAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTCAGGCTCAGGGCTT  >  NC_002947/77358‑77652
                                                                                                                                                  |                                                                                                                                                     
cagcagCTGTGCGGCCAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGc                                                                                                                                                   >  1:373664/1‑150 (MQ=255)
 agcagcTGTGCGGCCAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCg                                                                                                                                                  >  2:1206978/1‑150 (MQ=255)
 agcagcTGTGCGGCCAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCg                                                                                                                                                  >  2:1027008/1‑150 (MQ=255)
     gcTGTGCGGCCAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGTGCTCAGCGGcaa                                                                                                                                              <  1:1309490/150‑3 (MQ=255)
         tgCGGCCAGCAGGTCGACGCGTGGCAGACAGCGTGTGGTACGGAACGCGAAGGGCTGGGTGAAGCCGAGGTGGCGTTCATCTTCTTTATCTAGTTGCGCCGGCCGCGAGCTGTTCAGAGAAGCCTTGCCGCTGAGCACAGCGGCCAgcgc                                                                                                                                          >  2:521405/1‑150 (MQ=255)
             gCCAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACgg                                                                                                                                      <  1:676620/150‑1 (MQ=255)
              ccAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGTGCTCAGCGGCAAg                                                                                                                                             >  2:539662/1‑142 (MQ=255)
              ccAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGTGCTCAGCGGCAAg                                                                                                                                             <  1:539662/142‑1 (MQ=255)
                 gcagGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGc                                                                                                                                                >  1:1177118/1‑136 (MQ=255)
                 gcagGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGc                                                                                                                                                <  2:1177118/136‑1 (MQ=255)
                 gcagGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGggcgg                                                                                                                                  >  2:139504/1‑150 (MQ=255)
                               ggCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGTCGACGCCACGGTCATCTACCTTGTCCAGTTGAACCTGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTggg                                                                                                                    <  1:1385987/150‑1 (MQ=255)
                                  tGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTc                                                                                                                 >  2:150515/1‑150 (MQ=255)
                                   cctggCGTTCGGTACCGCACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGTGCTCAGCGGcaa                                                                                                                                              >  1:1405772/5‑118 (MQ=255)
                                   cctggCGTTCGGTACCGCACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGTGCTCAGCGGcaa                                                                                                                                              <  2:1405772/116‑3 (MQ=255)
                                     cAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCAtt                                                                                                              <  1:1376344/150‑1 (MQ=255)
                                     cAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCAtt                                                                                                              >  1:183580/1‑150 (MQ=255)
                                     cAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCAtt                                                                                                              >  2:684517/1‑150 (MQ=255)
                                      aGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTg                                                                                                             <  2:268208/150‑1 (MQ=255)
                                       gCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTgg                                                                                                            <  1:150515/150‑1 (MQ=255)
                                             ggTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGc                                                                                                      >  2:778926/1‑150 (MQ=255)
                                               tACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGt                                                                                                                       <  2:812542/131‑1 (MQ=255)
                                               tACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTACAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGt                                                                                                                       >  1:812542/1‑131 (MQ=255)
                                                    aaCGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGc                                                                                               >  2:1077653/1‑150 (MQ=255)
                                                     aCGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGcc                                                                                              <  1:213348/150‑1 (MQ=255)
                                                          aGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCcggcgg                                                                                         <  1:684517/150‑1 (MQ=255)
                                                             tCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGa                                                                                      >  2:390006/1‑150 (MQ=255)
                                                              ccTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGaa                                                                                     <  1:584928/150‑1 (MQ=255)
                                                               cTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAGGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGTCCCAGCACCTTGCCGGCGGCGAAc                                                                                    >  1:1203069/1‑150 (MQ=255)
                                                                  ggTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTtc                                                                                 >  1:171732/1‑150 (MQ=255)
                                                                        gCCGCGGTCGCGTTCATCTACCTTTTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGCACTTCTCgcgg                                                                           >  1:57277/1‑150 (MQ=255)
                                                                          cgcgGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCtt                                                                                                 <  1:608382/126‑1 (MQ=255)
                                                                          cgcgGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCtt                                                                                                 >  2:608382/1‑126 (MQ=255)
                                                                           gcgGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGcggcgg                                                                        >  1:815826/1‑150 (MQ=255)
                                                                              gTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTg                                                                     >  2:1352984/1‑150 (MQ=255)
                                                                                 gcgTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTg                                                                  <  2:171732/150‑1 (MQ=255)
                                                                                   gTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAAc                                                                                    >  2:1089327/1‑130 (MQ=255)
                                                                                   gTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAAc                                                                                    <  1:1089327/130‑1 (MQ=255)
                                                                                    tgtaTCTACCTTGTGCAGTTGCTCCTGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAAc                                                               >  2:1090560/4‑150 (MQ=255)
                                                                                        tCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCTTCGGCTGTTAAGCAGGt                                                           >  2:793003/1‑150 (MQ=255)
                                                                                        tCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGt                                                           <  2:815826/150‑1 (MQ=255)
                                                                                         cTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTc                                                          >  2:82842/1‑150 (MQ=255)
                                                                                           aCCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCAc                                                        >  2:537076/1‑150 (MQ=255)
                                                                                            ccTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACt                                                       >  1:986301/1‑150 (MQ=255)
                                                                                            ccTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACt                                                       <  1:96624/150‑1 (MQ=255)
                                                                                              ttGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGa                                                     >  1:1090990/1‑150 (MQ=255)
                                                                                               tGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGaa                                                    <  2:1090990/150‑1 (MQ=255)
                                                                                               tGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGaa                                                    <  1:139504/150‑1 (MQ=255)
                                                                                                gTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGACGTGTTCGCCGGCAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAg                                                   <  1:1027008/150‑1 (MQ=255)
                                                                                                      ttGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAAt                                             <  1:1206978/150‑1 (MQ=255)
                                                                                                              gCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGACGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTgcc                                     >  1:110703/1‑150 (MQ=255)
                                                                                                                  cGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGCACAGGTCACTGAAGCTCAATGCCTTGCCGCCa                                 <  1:741672/150‑1 (MQ=255)
                                                                                                                        gTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTgcc                                     >  2:492243/1‑140 (MQ=255)
                                                                                                                        gTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTgcc                                     <  1:492243/140‑1 (MQ=255)
                                                                                                                         ttCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGtt                          >  1:947976/1‑150 (MQ=255)
                                                                                                                         ttCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACCTCTTTCCTGTGCCGGCTAACTGGTGAGGGAAGCTCAATGCCTTGCCGCCAGGCAGtt                          <  1:603991/150‑1 (MQ=255)
                                                                                                                           cAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGc                        <  1:476299/150‑1 (MQ=255)
                                                                                                                             gCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAg                      >  2:438920/1‑150 (MQ=255)
                                                                                                                              cGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAgg                     >  1:504652/1‑150 (MQ=255)
                                                                                                                                aaGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCg                   >  1:136105/1‑150 (MQ=255)
                                                                                                                                    cTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTc               >  1:744995/1‑150 (MQ=255)
                                                                                                                                     ttGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTCa              >  2:901623/1‑150 (MQ=255)
                                                                                                                                      tggCCGCTGTGCTCAGCGGCCAGCGCACGGGCGGCGGCGAGGAGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTCa              >  2:392049/3‑150 (MQ=255)
                                                                                                                                         ccGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAgg                               >  2:1154629/1‑130 (MQ=255)
                                                                                                                                         ccGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGCTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAgg                               <  1:1154629/130‑1 (MQ=255)
                                                                                                                                          cGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTCAGGCt          >  1:893811/1‑150 (MQ=255)
                                                                                                                                                  acaGCGGCCAGCGCACGGGCGGCGGCGCGCGGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTCAGGCTCAGGGCtt  <  1:1071008/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CAGCAGCTGTGCGGCCAGCAGGTCGAGGCGTGGCTGCCAGCGTGCGGTACCGAACGCCAGGTCCTGGGTCAGGCCGCGGTCGCGTTCATCTACCTTGTCCAGTTGCGCCGGCAGCGAGCTGTTCAGCGAAGCCTTGCCGCTGAGCACAGCGGCCAGCGCACGGGCGGCGGCGAGGCGTGGGTTCATTGGCCCAGCACCTTGCCGGCGGCGAACTTCTCGCGGCGGCTGTTGAACAGGTCACTGAAGCTCAATGCCTTGCCGCCAGGCAGTTGCAGGCGGGTCAGGCTCAGGGCTT  >  NC_002947/77358‑77652

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: