Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,606,642 C→T 6.9% intergenic (‑241/+27) leuQ ← / ← rsmC tRNA‑Leu/16S rRNA m(2)G1207 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,606,6420CT 6.9% 271.8 / 6.8 87intergenic (‑241/+27)leuQ/rsmCtRNA‑Leu/16S rRNA m(2)G1207 methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (36/45);  new base T (2/4);  total (38/49)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

AATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGG  >  NC_000913/4606494‑4606787
                                                                                                                                                    |                                                                                                                                                 
aaTTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCa                                                                                                                                                  <  2:1132069/150‑1 (MQ=255)
       ggAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGaca                                                                                                                                           <  2:1169179/150‑1 (MQ=255)
         aaaaaTGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATGCGGCGTGAACGCCTTATCCGaccta                                                                                                                                         <  1:1347163/150‑4 (MQ=255)
         aaaaaTGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGacatg                                                                                                                                         >  2:655366/1‑148 (MQ=255)
         aaaaaTGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTCATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACg                                                                                                                                         <  1:1183780/150‑1 (MQ=255)
          aaaaTGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATGCGGCGTGAACGCCTTATCCGacctac                                                                                                                                        <  1:1006466/150‑5 (MQ=37)
          aaaaTGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATGCGGCGTGAACGCCTTATCCGacctac                                                                                                                                        <  1:1116252/150‑5 (MQ=37)
           aaaTGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGt                                                                                                                                       <  2:600968/150‑1 (MQ=255)
           aaaTGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGt                                                                                                                                       <  2:604673/150‑1 (MQ=255)
                   aTCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCg                                                                                                                               <  1:1416682/150‑1 (MQ=255)
                     cAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCgag                                                                                                                             <  1:1114420/150‑1 (MQ=255)
                     cAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCgag                                                                                                                             >  1:325974/1‑150 (MQ=255)
                      aGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCgaga                                                                                                                            <  2:1081604/150‑1 (MQ=255)
                       ggAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGAt                                                                                                                           >  2:1508252/1‑150 (MQ=255)
                       ggAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGAt                                                                                                                           >  1:107921/1‑150 (MQ=255)
                             tAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTaa                                                                                                                     <  1:856144/150‑1 (MQ=255)
                               gCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATTTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAAcc                                                                                                                   <  1:1148830/150‑1 (MQ=255)
                               gCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAAcc                                                                                                                   <  1:85890/150‑1 (MQ=255)
                                   gCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTttc                                                                                                               >  2:634776/1‑150 (MQ=255)
                                    cAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTttct                                                                                                              >  2:256139/1‑150 (MQ=255)
                                       gTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCg                                                                                                           <  1:1542063/150‑1 (MQ=255)
                                       gTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCCGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCg                                                                                                           <  2:5907/150‑1 (MQ=255)
                                           ggTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTg                                                                                                       >  2:1061544/1‑150 (MQ=255)
                                           ggTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTg                                                                                                       >  1:86069/1‑150 (MQ=255)
                                               tgcatgcGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCgg                                                                                                   <  2:530057/150‑1 (MQ=255)
                                                gcatgcGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCggg                                                                                                  <  1:1292263/150‑1 (MQ=255)
                                                 catgcGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGt                                                                                                 >  2:875330/1‑150 (MQ=255)
                                                    gcGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCAt                                                                                              >  2:464563/1‑150 (MQ=255)
                                                      gggggCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATaa                                                                                            <  2:146633/150‑1 (MQ=255)
                                                         ggCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATgg                                                                                         >  1:194102/1‑150 (MQ=255)
                                                           cagatgccagatgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGgcg                                                                                       <  1:634776/150‑1 (MQ=255)
                                                            agatgccagatgcGACGCTGGCGCGGCTTATCTGGCCGACGAAGGGCTAACGTGCAGGGTTTGTAGGGCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGgcgc                                                                                      <  2:1307719/150‑1 (MQ=255)
                                                                 ccagatgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATac                                                                                 <  2:1357394/150‑1 (MQ=255)
                                                                    gatgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACAcc                                                                              <  1:1394349/150‑1 (MQ=255)
                                                                     atgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCt                                                                             <  1:1094456/150‑1 (MQ=255)
                                                                        cGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGa                                                                          <  1:629088/150‑1 (MQ=255)
                                                                         gACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGaa                                                                         <  1:229478/150‑1 (MQ=255)
                                                                             cTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGgcgc                                                                                      <  2:1445915/133‑1 (MQ=255)
                                                                             cTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGgcgc                                                                                      >  1:1445915/1‑133 (MQ=255)
                                                                               ggCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGccc                                                                   >  1:1519845/1‑150 (MQ=255)
                                                                               ggCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGccc                                                                   >  1:206737/1‑150 (MQ=255)
                                                                                    gTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGttt                                                              <  2:325974/150‑1 (MQ=255)
                                                                                     tCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTg                                                             >  1:717156/1‑150 (MQ=255)
                                                                                     acTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCCTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTg                                                             >  1:941625/2‑150 (MQ=255)
                                                                                       ttATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTgcg                                                           >  2:806682/1‑150 (MQ=255)
                                                                                         aTCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGcgcg                                                         <  1:256139/150‑1 (MQ=255)
                                                                                            tGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATc                                                      <  1:464563/150‑1 (MQ=255)
                                                                                            tGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATc                                                      >  2:1123052/1‑150 (MQ=255)
                                                                                             ggCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCa                                                     >  2:947894/1‑150 (MQ=255)
                                                                                                cTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACtt                                                  <  1:1508252/150‑1 (MQ=255)
                                                                                                 tACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTc                                                 >  1:805991/1‑150 (MQ=255)
                                                                                                 tACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTc                                                 >  1:271377/1‑150 (MQ=255)
                                                                                                  aCGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCg                                                <  1:372919/150‑1 (MQ=255)
                                                                                                  aCGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCg                                                <  1:272813/150‑1 (MQ=255)
                                                                                                   cGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGt                                               <  1:1538873/150‑1 (MQ=255)
                                                                                                    gAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCACCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTg                                              <  1:558381/150‑1 (MQ=255)
                                                                                                         gCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGACGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGtt                                                               <  2:167308/128‑1 (MQ=255)
                                                                                                         gCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGACGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGtt                                                               >  1:167308/1‑128 (MQ=255)
                                                                                                                  cAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGtct                                >  2:268435/1‑150 (MQ=255)
                                                                                                                  cAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGtct                                <  1:1154836/150‑1 (MQ=255)
                                                                                                                    ggTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAg                                          <  1:993882/138‑1 (MQ=255)
                                                                                                                    ggTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAg                                          >  2:993882/1‑138 (MQ=255)
                                                                                                                    ggTTTTGTAGGTCGGATAAGGCGTACACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCa                              >  2:801960/1‑150 (MQ=255)
                                                                                                                     gTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCAt                             >  1:885093/1‑150 (MQ=255)
                                                                                                                     gTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCAt                             <  1:382093/150‑1 (MQ=255)
                                                                                                                      ttttGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCagg                            >  1:806686/1‑148 (MQ=255)
                                                                                                                          gTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACATCTTGAATCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGc                        >  2:1587068/1‑150 (MQ=255)
                                                                                                                            aGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCAc                      >  1:490657/1‑150 (MQ=255)
                                                                                                                            aGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCAc                      >  2:295505/1‑150 (MQ=255)
                                                                                                                                cGGATGCGGCGTGAACGCCTTATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCAAGCACGTcc                  >  2:1032909/1‑150 (MQ=255)
                                                                                                                                  gatgcggCGTGAACGCCTTATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCgg                >  2:796247/6‑150 (MQ=255)
                                                                                                                                   aTAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGa               <  1:875330/150‑1 (MQ=255)
                                                                                                                                    tAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGAt              <  2:717156/150‑1 (MQ=255)
                                                                                                                                         cGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACATCTTGAATCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGc         <  1:1587068/150‑1 (MQ=255)
                                                                                                                                         cGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATg                                   <  2:593743/124‑1 (MQ=255)
                                                                                                                                         cGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATg                                   >  1:593743/1‑124 (MQ=255)
                                                                                                                                         cGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGc         <  2:1363663/150‑1 (MQ=255)
                                                                                                                                         cGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGc         >  1:955932/1‑150 (MQ=255)
                                                                                                                                          gTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCa        <  2:643555/150‑1 (MQ=255)
                                                                                                                                           ttCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAg       <  2:1169941/150‑1 (MQ=255)
                                                                                                                                           ttCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAg       <  1:1328897/150‑1 (MQ=255)
                                                                                                                                           ttCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAg       >  1:1339081/1‑150 (MQ=255)
                                                                                                                                            tCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAgg      <  1:1123052/150‑1 (MQ=255)
                                                                                                                                            tCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAgg      >  1:1327681/1‑150 (MQ=255)
                                                                                                                                             cACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGTTAAGGCAGGa     >  2:317248/1‑150 (MQ=255)
                                                                                                                                             cACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGa     >  2:505601/1‑150 (MQ=255)
                                                                                                                                             aaCGCCTTATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCTAGCACGTCCGGATAAGGCAGGa     <  1:1032909/149‑1 (MQ=255)
                                                                                                                                               cgcctcATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCa                       >  1:948578/1‑130 (MQ=255)
                                                                                                                                               cgcctcATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCa                       <  2:948578/130‑1 (MQ=255)
                                                                                                                                               cgccgcATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAg   >  2:1172420/1‑150 (MQ=255)
                                                                                                                                               cgccgcATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAg   >  2:205814/1‑150 (MQ=255)
                                                                                                                                                gccgcATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAgg  >  1:713275/1‑150 (MQ=255)
                                                                                                                                                  cgcATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAg       <  1:63356/143‑1 (MQ=255)
                                                                                                                                                  cgcATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAg       >  2:63356/1‑143 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGG  >  NC_000913/4606494‑4606787

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