Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,920,205 (T)7→8 intergenic (+62/‑18) yebT → / → rsmF MCE domain protein/16S rRNA m(5)C1407 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,920,1981.T97.3% 230.6 / ‑1.1 73intergenic (+55/‑25)yebT/rsmFMCE domain protein/16S rRNA m(5)C1407 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base T (37/34);  total (38/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 = 7.34e-01

GGAAGGCAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCG  >  NC_000913/1920068‑1920330
                                                                                                                                   |                                                                                                                                    
ggAAGGCAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGtttttttt                                                                                                                              >  1:1202598/1‑138 (MQ=255)
      cAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGt                                                                                                                         <  2:1202598/139‑1 (MQ=255)
            cTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACa                                                                                                                   >  1:1855321/1‑139 (MQ=255)
                     gCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTgg                                                                                                          >  1:64793/1‑139 (MQ=255)
                       aGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAg                                                                                                                          >  1:700301/1‑121 (MQ=255)
                       aGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAg                                                                                                                          <  2:700301/121‑1 (MQ=255)
                              tGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAacaca                                                                                                 >  1:1956142/1‑139 (MQ=255)
                                                 gAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGc                                                                              <  2:64793/139‑1 (MQ=255)
                                                 gAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGc                                                                              <  2:1956142/139‑1 (MQ=255)
                                                       ggAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAg                                                                                                                          <  1:1417048/89‑1 (MQ=255)
                                                       ggAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAg                                                                                                                          >  2:1417048/1‑89 (MQ=255)
                                                       ggAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGcc                                                                                                        >  2:441942/1‑107 (MQ=255)
                                                       ggAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGcc                                                                                                        <  1:441942/107‑1 (MQ=255)
                                                           cTGCGCTTCCCAAATAATGCCCAATGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGaca                                                                    <  1:1390423/139‑1 (MQ=255)
                                                            tGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTc                                                                         <  2:199144/133‑1 (MQ=255)
                                                            tGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTc                                                                         >  1:199144/1‑133 (MQ=255)
                                                              cgcTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTg                                                                                                           <  2:344810/97‑1 (MQ=255)
                                                              cgcTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTg                                                                                                           >  1:1672235/1‑97 (MQ=255)
                                                              cgcTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTg                                                                                                           >  1:344810/1‑97 (MQ=255)
                                                              cgcTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTg                                                                                                           <  2:1672235/97‑1 (MQ=255)
                                                              cgcTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGttttttttgcc                                                                                                                           <  1:1101659/80‑5 (MQ=255)
                                                              cgcTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGttttttttgcc                                                                                                                           >  2:1101659/1‑76 (MQ=255)
                                                                 ttCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATg                                                              >  1:1068437/1‑139 (MQ=255)
                                                                      aaaTAATGTCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCg                                                                                   <  1:368926/113‑1 (MQ=255)
                                                                      aaaTAATGTCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCg                                                                                   >  2:368926/1‑113 (MQ=255)
                                                                             gCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCt                                                                                                               >  2:1971811/1‑78 (MQ=255)
                                                                             gCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCt                                                                                                               <  1:1971811/78‑1 (MQ=255)
                                                                             gCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAacacac                                                                                                <  2:1747047/93‑1 (MQ=255)
                                                                             gCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAacacac                                                                                                >  1:1747047/1‑93 (MQ=255)
                                                                             gCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGAt                                                   <  1:543684/138‑1 (MQ=255)
                                                                             gCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGAt                                                   >  2:543684/1‑138 (MQ=255)
                                                                              cccACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGAc                                                                                >  2:1520096/1‑108 (MQ=255)
                                                                              cccACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGAc                                                                                <  1:1520096/108‑1 (MQ=255)
                                                                              cccACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGttttttt                                                                                                                               <  2:969240/60‑1 (MQ=255)
                                                                              cccACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGttttttt                                                                                                                               >  1:969240/1‑60 (MQ=255)
                                                                                 aCTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACg                                                                               >  2:956830/1‑106 (MQ=255)
                                                                                 aCTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACg                                                                               <  1:956830/106‑1 (MQ=255)
                                                                                    gCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGa                                           <  2:1855321/139‑1 (MQ=255)
                                                                                        cGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGct                                       <  2:1068437/139‑1 (MQ=255)
                                                                                           cGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGcct                                                                                                       <  2:453049/72‑2 (MQ=255)
                                                                                           cGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGcct                                                                                                       >  1:453049/1‑71 (MQ=255)
                                                                                             tGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAac                                                                                                    >  1:1446314/1‑73 (MQ=255)
                                                                                             tGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAac                                                                                                    <  2:1446314/73‑1 (MQ=255)
                                                                                             tGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGctc                                      >  1:505319/1‑135 (MQ=255)
                                                                                             tGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGctc                                      <  2:505319/135‑1 (MQ=255)
                                                                                             tGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCAt                                  >  1:1280692/1‑139 (MQ=255)
                                                                                             tGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCAt                                  >  1:225590/1‑139 (MQ=255)
                                                                                             tGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCAt                                  >  2:939035/1‑139 (MQ=255)
                                                                                              gCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCAtt                                 <  1:939035/139‑1 (MQ=255)
                                                                                               ccTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGttttttttgc                                                                                                                            >  2:307209/1‑43 (MQ=255)
                                                                                               ccTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGttttttttgc                                                                                                                            <  1:307209/46‑4 (MQ=255)
                                                                                               ccTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTGCCAGTGGTACATGCTCGTg                                                                                                           <  1:1212536/63‑1 (MQ=255)
                                                                                               ccTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTGCCAGTGGTACATGCTCGTg                                                                                                           >  2:1212536/1‑63 (MQ=255)
                                                                                                cTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAg                                                                                                                          <  2:1453953/48‑1 (MQ=39)
                                                                                                cTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAg                                                                                                                          >  1:1453953/1‑48 (MQ=255)
                                                                                                    gcCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGa                                                                                 >  2:41839/1‑85 (MQ=255)
                                                                                                    gcCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGa                                                                                 <  1:41839/85‑1 (MQ=255)
                                                                                                        gAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACa                                                                                                                   >  1:238708/1‑47 (MQ=255)
                                                                                                        gAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACa                                                                                                                   <  2:238708/47‑1 (MQ=255)
                                                                                                         aGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTg                                                                                                           >  2:1633800/1‑54 (MQ=255)
                                                                                                         aGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTg                                                                                                           <  1:1633800/54‑1 (MQ=255)
                                                                                                         aGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTa                                                                                          >  2:472607/1‑71 (MQ=255)
                                                                                                         aGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTa                                                                                          <  1:472607/71‑1 (MQ=255)
                                                                                                         aGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCg                                                     <  2:777671/108‑1 (MQ=255)
                                                                                                         aGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCg                                                     >  1:777671/1‑108 (MQ=255)
                                                                                                                     aaaCTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTTACACAAATgcg                                                            >  1:1012143/1‑89 (MQ=255)
                                                                                                                     aaaCTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATgcg                                                            <  2:1012143/89‑1 (MQ=255)
                                                                                                                        cTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGa                                           >  2:250647/1‑103 (MQ=255)
                                                                                                                        cTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGa                                           <  1:250647/103‑1 (MQ=255)
                                                                                                                           cgcgcCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCg                                                                                   <  1:1520235/60‑1 (MQ=255)
                                                                                                                           cgcgcCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCg                                                                                   >  2:1520235/1‑60 (MQ=255)
                                                                                                                           cgcgcCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGcc    >  2:350423/1‑139 (MQ=255)
                                                                                                                           cgcgcCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGcc    >  1:612392/1‑139 (MQ=255)
                                                                                                                           cgcgcCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTAGATGATTTTCTTGCCGCCTGTCAGCGcc    >  2:1234055/1‑139 (MQ=255)
                                                                                                                            gcgcCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTtgat                             >  2:887619/1‑113 (MQ=255)
                                                                                                                            gcgcCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTtgat                             <  1:887619/113‑1 (MQ=255)
                                                                                                                             cgcCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCg                                            >  1:1747938/1‑97 (MQ=255)
                                                                                                                             cgcCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCg                                            <  2:1747938/97‑1 (MQ=255)
                                                                                                                             cgcCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCg  >  2:431500/1‑139 (MQ=255)
                                                                                                                             cgcCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCg  <  1:350423/139‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                    
GGAAGGCAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCG  >  NC_000913/1920068‑1920330

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: