Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,162,152 G→C 6.4% intergenic (+22/+18) 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,1520GC 6.4% 395.8 / 15.8 126intergenic (+22/+18)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 G (67/51);  new base C (5/3);  total (72/54)
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.87e-01

Only 100 of 135 total aligned reads displayed.

CAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTCCGGGTTACA  >  NC_000913/4162004‑4162301
                                                                                                                                                    |                                                                                                                                                     
cAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGgc                                                                                                                                                      >  2:532277/1‑150 (MQ=255)
cAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGgc                                                                                                                                                      >  2:368118/1‑150 (MQ=255)
cAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGgc                                                                                                                                                      >  1:131501/1‑150 (MQ=255)
 aGCACGGTAGTTCGTGCGGAATCTGCGGATTTCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATAATGGTGTTCTGGTTTGGCGCGAAGATGCGTAACCGTAAGCAGGTAGTTGCTGAGTAATAGGTTTTATGCCGGATGCTGgcg                                                                                                                                                     <  2:528272/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)
                gCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGtt                                                                                                                                      >  1:961253/1‑150 (MQ=255)
                   gAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTa                                                                                                                                   >  2:411315/1‑150 (MQ=255)
                     atatCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACt                                                                                                                                 >  1:1697865/1‑150 (MQ=255)
                       atCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTc                                                                                                                               >  2:1465837/1‑150 (MQ=255)
                        tCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTcg                                                                                                                              >  1:1568506/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)
                               atCGGCTCACACTTCTTCCCGGCAGTCCGGTCCGTCATTATGGTGGTCTGGAGTGGCGCGAAGATCCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAg                                                                                                                       <  1:370436/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)
                                          cttcttCCCGGCAGTACTCTCCGTCATTATGGGGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTATTCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATAcgccgc                                                                                                            <  1:165545/150‑1 (MQ=255)
                                           ttcttcCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATAcgccgcc                                                                                                           >  2:1409697/1‑149 (MQ=255)
                                            tcttcCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCAc                                                                                                          >  1:808015/1‑150 (MQ=255)
                                             cttcCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACt                                                                                                         >  2:840953/1‑150 (MQ=255)
                                             cttcCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACt                                                                                                         >  1:590233/1‑150 (MQ=255)
                                             cttcCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACt                                                                                                         >  2:876566/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)
                                                              ccGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAg                                                                                        >  2:1511505/1‑150 (MQ=255)
                                                                  cATTATGGTGTTCAGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGa                                                                                    >  2:1049034/1‑150 (MQ=255)
                                                                     tATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACt                                                                                 <  1:653720/150‑1 (MQ=255)
                                                                      aTGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCTGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTg                                                                                <  2:1368806/150‑1 (MQ=255)
                                                                      aTGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTg                                                                                <  2:1613826/150‑1 (MQ=255)
                                                                        ggTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGAt                                                                              >  2:1862236/1‑150 (MQ=255)
                                                                        ggTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGAt                                                                              <  1:532277/150‑1 (MQ=255)
                                                                          tgtTCTGGTTTGGCGCGACTATGCGTAACCGTAAGCATGAAGTTGCTGAGTAATCGGTATTTTGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCa                                                                            <  2:751279/150‑1 (MQ=255)
                                                                           gtTCTGGATTTGCGCGAATATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCaa                                                                           <  1:55611/150‑1 (MQ=255)
                                                                           gtTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCaa                                                                           >  2:1562178/1‑150 (MQ=255)
                                                                              cTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACa                                                                        >  1:248696/1‑150 (MQ=255)
                                                                              cTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACa                                                                        >  1:1289734/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)
                                                                                  aTTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGc                                                                    >  1:27233/1‑150 (MQ=255)
                                                                                   ttGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGcc                                                                   <  1:18906/150‑1 (MQ=255)
                                                                                      gcgcGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGa                                                                >  1:1652662/1‑150 (MQ=255)
                                                                                       cgcgAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGAc                                                               >  2:1119865/1‑150 (MQ=255)
                                                                                       cgcgAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGAc                                                               >  2:1402471/1‑150 (MQ=255)
                                                                                              aTGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGAc                                                        <  2:210335/150‑1 (MQ=255)
                                                                                                   tAACCGAAAGCAGGAAGTTGCTGACTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTtgt                                                   <  2:468547/150‑1 (MQ=255)
                                                                                                    aaCCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTtgtg                                                  <  1:1862236/150‑1 (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:251855/150‑1 (MQ=255)
                                                                                                             cAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCt                                         >  1:243544/1‑150 (MQ=255)
                                                                                                              aGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACATAGCCAGACGTTCGACCTTGTGCGTCTGGCtc                                        <  1:1201536/150‑2 (MQ=255)
                                                                                                              aGGAAGTGGCTGAGTAAGCGGCAGTATGCCGGCTGCTGGCGCAGCCGGCATGGGGTTGACTTCGCGGGCAGTAATCCGCCGCACTCCATCGGGGGCGTGTAGGGGAACTGATCAACCAGAGCCAGACGGGCGACCTTGGGCGTCTGGcgc                                        <  1:1201541/150‑3 (MQ=255)
                                                                                                                 aaGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAAt                                     <  1:1079278/150‑1 (MQ=255)
                                                                                                                 aaGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAAt                                     >  1:848099/1‑150 (MQ=255)
                                                                                                                       cTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTcc                               <  2:4973/150‑1 (MQ=255)
                                                                                                                           gTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGAt                           >  1:518210/1‑150 (MQ=255)
                                                                                                                             aaTCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTc                         >  1:1753141/1‑150 (MQ=255)
                                                                                                                             aaTCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACATTGTGCGTCTGGCTTAATGTTTCCAGATTc                         >  2:1227759/1‑150 (MQ=255)
                                                                                                                              aTCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCt                        >  1:1686718/1‑150 (MQ=255)
                                                                                                                              aTCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGAAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCt                        >  1:1136098/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)
                                                                                                                                     ccatGCCGGATGCGCCAGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATaa                 >  1:680416/3‑150 (MQ=255)
                                                                                                                                      tatGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCTTATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAAc                <  2:1393651/150‑1 (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)
                                                                                                                                         gCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGtt             >  2:882524/1‑150 (MQ=255)
                                                                                                                                         gCCGGATGCGCCAGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGtt             >  1:1468499/1‑150 (MQ=255)
                                                                                                                                          ccGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGttt            <  1:1294702/150‑1 (MQ=255)
                                                                                                                                           cGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCTTATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTc           <  1:1343629/150‑1 (MQ=255)
                                                                                                                                           cGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTc           <  2:471890/150‑1 (MQ=255)
                                                                                                                                            ggATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTcc          >  1:316734/1‑150 (MQ=255)
                                                                                                                                                gCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGTGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTCCGGGt      >  2:220298/1‑150 (MQ=255)
                                                                                                                                                 cTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTCCGGGtt     >  1:1799140/1‑150 (MQ=255)
                                                                                                                                                  tGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTCCGGGTTa    >  2:687569/1‑150 (MQ=255)
                                                                                                                                                   ggCGCATGCGGCATCGGGTTTACTTCGCGGTCAGTAATACGCCGGACTCGGTAGGGTGCGTGTAGGGCAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAAAGTTTCCAGATTCTAGCAGAACGTTTCCGGGTTa  >  2:1139161/1‑149 (MQ=255)
                                                                                                                                                   ggCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTCCGGGTTAc   >  2:795184/1‑150 (MQ=255)
                                                                                                                                                    gcgcATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTCCGGGTTACa  <  1:785960/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTCCGGGTTACA  >  NC_000913/4162004‑4162301

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: