Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,162,150 T→G 6.1% intergenic (+20/+20) yijD → / ← trmA DUF1422 family inner membrane protein/tRNA m(5)U54 methyltransferase, SAM‑dependent; tmRNA m(5)U341 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,162,1500TG 6.1% 391.3 / 12.5 131intergenic (+20/+20)yijD/trmADUF1422 family inner membrane protein/tRNA m(5)U54 methyltransferase, SAM‑dependent; tmRNA m(5)U341 methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (70/52);  new base G (5/3);  total (76/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.12e-01

Only 100 of 132 total aligned reads displayed.

TACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTCCGGGTTA  >  NC_000913/4162002‑4162299
                                                                                                                                                    |                                                                                                                                                     
tACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTg                                                                                                                                                      >  2:845371/1‑150 (MQ=255)
  cAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGgc                                                                                                                                                    >  2:368118/1‑150 (MQ=255)
  cAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGgc                                                                                                                                                    >  1:131501/1‑150 (MQ=255)
  cAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGgc                                                                                                                                                    >  2:1294702/1‑150 (MQ=255)
   aGCACGGTAGTTCGTGCGGAATCTGCGGATTTCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATAATGGTGTTCTGGTTTGGCGCGAAGATGCGTAACCGTAAGCAGGTAGTTGCTGAGTAATAGGTTTTATGCCGGATGCTGgcg                                                                                                                                                   <  2:528272/150‑1 (MQ=255)
     cACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCa                                                                                                                                                 <  1:1199279/150‑1 (MQ=255)
     cACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCa                                                                                                                                                 >  1:1636984/1‑150 (MQ=255)
       cGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATc                                                                                                                                               >  1:1393651/1‑150 (MQ=255)
        ggTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATcc                                                                                                                                              <  1:608343/150‑1 (MQ=255)
         gTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGAACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCg                                                                                                                                             <  2:1762227/150‑1 (MQ=255)
           aGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTGCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGc                                                                                                                                           <  1:1156081/150‑1 (MQ=255)
             ttCGTGCGGAATATCCGGATCTCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCGCCAGCATCCGGCAt                                                                                                                                         <  1:650506/150‑1 (MQ=255)
             ttCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTGCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATAATGCCGGATGCGCCAGCATCCGGCAt                                                                                                                                         <  1:1106363/150‑1 (MQ=255)
              tCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATg                                                                                                                                        >  2:1558723/1‑150 (MQ=255)
                gTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATggg                                                                                                                                      <  1:111305/150‑1 (MQ=255)
                  gCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGTTGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCTTCCGGCATgcttt                                                                                                                                    >  1:168556/1‑146 (MQ=255)
                        tataCGGATATCGGCTCTAACTTCTTCACGGCAGTACTCTACGACATTATGGTGTTCTGGATTGGCGAGAAGAAGCGTAACCGTAAGCAGGAAGTTGCTGAGTAAACGGTATAAAGCCGGATGCTGGCGCATCCGGCATGGGTTTTACtt                                                                                                                              >  2:443509/1‑150 (MQ=255)
                         atCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTc                                                                                                                             >  2:1465837/1‑150 (MQ=255)
                          tCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTcgc                                                                                                                           >  2:622564/1‑150 (MQ=255)
                          tCCGGATATCGGCTCTACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTcgc                                                                                                                           >  2:1150916/1‑150 (MQ=255)
                           ccGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTcgc                                                                                                                           <  1:1091125/150‑1 (MQ=255)
                            cGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTcgcg                                                                                                                          >  2:1343629/1‑150 (MQ=255)
                            cGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGACATAATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTcgcg                                                                                                                          >  1:1469014/1‑150 (MQ=255)
                             ggATCTCGGCTCTTACTTCTTACCGGCAGTACTCTCCGTCATTAGGGTGTTCTGGAGTGGCGCGAAGATGCGTTACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCgg                                                                                                                         <  2:931016/150‑1 (MQ=255)
                                 atCGGCTCTAACTTCTTCCCGGCAGAACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAg                                                                                                                     <  1:97273/150‑1 (MQ=255)
                                   cGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTAATTCGCGGTCAGTa                                                                                                                   >  2:1428103/1‑150 (MQ=255)
                                          aaCTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATAcgcc                                                                                                            <  2:76741/150‑1 (MQ=255)
                                           aCTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATAcgccg                                                                                                           >  2:654925/1‑150 (MQ=255)
                                            cttcttCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATAcgccgc                                                                                                          >  2:408867/1‑150 (MQ=255)
                                             ttcttcCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATAcgccgcc                                                                                                         >  2:1409697/1‑149 (MQ=255)
                                              tcttcCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCAc                                                                                                        >  1:808015/1‑150 (MQ=255)
                                               cttcCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACt                                                                                                       >  2:876566/1‑150 (MQ=255)
                                               cttcCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACt                                                                                                       >  2:840953/1‑150 (MQ=255)
                                               cttcCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACt                                                                                                       >  1:590233/1‑150 (MQ=255)
                                                  cccGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCa                                                                                                    >  1:1141427/1‑150 (MQ=255)
                                                           aCTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCg                                                                                           <  1:368118/150‑1 (MQ=255)
                                                             tctcCGTCATTATGGTGTTCTGTATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTCTGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCgtg                                                                                         <  1:1504597/150‑1 (MQ=255)
                                                              ctcCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCgtgt                                                                                        >  1:561945/1‑150 (MQ=255)
                                                                ccGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAg                                                                                      >  2:1511505/1‑150 (MQ=255)
                                                                ccGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAg                                                                                      >  1:1246015/1‑150 (MQ=255)
                                                                 cGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAgg                                                                                     >  2:1005006/1‑150 (MQ=255)
                                                                       tATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACt                                                                               <  1:653720/150‑1 (MQ=255)
                                                                        aTGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCTGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTg                                                                              <  2:1368806/150‑1 (MQ=255)
                                                                        aTGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTg                                                                              >  1:285023/1‑150 (MQ=255)
                                                                        aTGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTg                                                                              <  2:1613826/150‑1 (MQ=255)
                                                                          ggTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGAt                                                                            >  2:1862236/1‑150 (MQ=255)
                                                                          ggTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGAt                                                                            >  2:1398720/1‑150 (MQ=255)
                                                                          ggTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGAt                                                                            <  1:532277/150‑1 (MQ=255)
                                                                            tgtTCTGGTTTGGCGCGACTATGCGTAACCGTAAGCATGAAGTTGCTGAGTAATCGGTATTTTGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCa                                                                          <  2:751279/150‑1 (MQ=255)
                                                                             gtTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCaa                                                                         >  1:1652227/1‑150 (MQ=255)
                                                                             gtTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCaa                                                                         >  2:1562178/1‑150 (MQ=255)
                                                                                cTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACa                                                                      >  1:248696/1‑150 (MQ=255)
                                                                                  ggATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACaga                                                                    <  2:1462050/150‑1 (MQ=255)
                                                                                   gATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTTCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACagag                                                                   <  2:1424014/150‑1 (MQ=255)
                                                                                   gATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACagag                                                                   >  2:1406692/1‑150 (MQ=255)
                                                                                        gcgcGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGa                                                              >  1:1652662/1‑150 (MQ=255)
                                                                                         cgcgAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGAc                                                             >  2:1402471/1‑150 (MQ=255)
                                                                                         cgcgAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGAc                                                             >  2:1119865/1‑150 (MQ=255)
                                                                                            gAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGtt                                                          >  1:225033/1‑150 (MQ=255)
                                                                                                aTGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGAc                                                      <  2:210335/150‑1 (MQ=255)
                                                                                                 tGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGAcc                                                     >  1:1550446/1‑150 (MQ=255)
                                                                                                     tAACCGAAAGCAGGAAGTTGCTGACTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTtgt                                                 <  2:468547/150‑1 (MQ=255)
                                                                                                      aaCCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTtgtg                                                <  2:615623/150‑1 (MQ=255)
                                                                                                      aaCCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTtgtg                                                <  2:393638/150‑1 (MQ=255)
                                                                                                      aaCCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTtgtg                                                <  1:1862236/150‑1 (MQ=255)
                                                                                                       aCCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGc                                               >  1:4973/1‑150 (MQ=255)
                                                                                                        ccGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCg                                              >  2:1616212/1‑150 (MQ=255)
                                                                                                         cGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGCCGTTCGACCTTGTGCGt                                             >  2:1805542/1‑150 (MQ=255)
                                                                                                          gTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTc                                            >  2:1185075/1‑150 (MQ=255)
                                                                                                           tacgCAGGAAGTTGCTGAGTAATCGGGATTATGCCGGATGCTGGCGCATCCGGCATGGGATTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCt                                           <  1:1211800/147‑1 (MQ=255)
                                                                                                            aaGCAGGAAGTTGCTGAGTAAGCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTg                                          <  2:906461/150‑1 (MQ=255)
                                                                                                              gCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGc                                        <  2:248696/150‑1 (MQ=255)
                                                                                                              gCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGc                                        >  2:707095/1‑150 (MQ=255)
                                                                                                               cAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCt                                       >  1:243544/1‑150 (MQ=255)
                                                                                                                aGGAAGTGGCTGAGTAAGCGGCAGTATGCCGGCTGCTGGCGCAGCCGGCATGGGGTTGACTTCGCGGGCAGTAATCCGCCGCACTCCATCGGGGGCGTGTAGGGGAACTGATCAACCAGAGCCAGACGGGCGACCTTGGGCGTCTGGcgc                                      <  1:1201541/150‑3 (MQ=255)
                                                                                                                   aaGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAAt                                   <  1:1079278/150‑1 (MQ=255)
                                                                                                                   aaGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAAt                                   >  1:848099/1‑150 (MQ=255)
                                                                                                                   aaGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAAt                                   <  1:1330039/150‑1 (MQ=255)
                                                                                                                    aGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATg                                  <  1:1465837/150‑1 (MQ=255)
                                                                                                                         cTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTcc                             <  2:4973/150‑1 (MQ=255)
                                                                                                                          tGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCATGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCa                            >  1:888503/1‑150 (MQ=255)
                                                                                                                             gTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGAt                         >  1:518210/1‑150 (MQ=255)
                                                                                                                               aaTCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTc                       >  1:1753141/1‑150 (MQ=255)
                                                                                                                                aTCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCt                      >  1:1686718/1‑150 (MQ=255)
                                                                                                                                  cGGTCTTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTg                    <  2:431655/150‑1 (MQ=255)
                                                                                                                                  cGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCTGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCTAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTg                    <  2:1469014/150‑1 (MQ=255)
                                                                                                                                   ggTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGc                   <  1:1370668/150‑1 (MQ=255)
                                                                                                                                    gTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCa                  <  1:1760113/150‑1 (MQ=255)
                                                                                                                                       ccatGCCGGATGCGCCAGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATaa               >  1:680416/3‑150 (MQ=255)
                                                                                                                                        tatGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAAc              <  2:1777783/150‑1 (MQ=255)
                                                                                                                                        catGCCGGATGCGCCAGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAAc              >  1:1683041/2‑150 (MQ=255)
                                                                                                                                         atGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACg             <  2:1636984/150‑1 (MQ=255)
                                                                                                                                         atGCCGGATGCGCCAGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGCTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACg             >  2:525715/1‑150 (MQ=255)
                                                                                                                                          tGCCGGATGCGCCAGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGt            >  1:98606/1‑150 (MQ=255)
                                                                                                                                           gCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGTACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGtt           >  2:888744/1‑150 (MQ=255)
                                                                                                                                            ccGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGttt          <  1:1294702/150‑1 (MQ=255)
                                                                                                                                              ggATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTcc        >  1:316734/1‑150 (MQ=255)
                                                                                                                                                  gCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGTGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTCCGGGt    >  2:220298/1‑150 (MQ=255)
                                                                                                                                                  gCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTCCGGGt    <  2:1141427/150‑1 (MQ=255)
                                                                                                                                                   cTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTCCGGGtt   >  1:1799140/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTCCGGGTTA  >  NC_000913/4162002‑4162299

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