Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,089,198 T→A 6.1% intergenic (‑70/+13) yeeZ ← / ← yoeB putative epimerase/toxin of the YoeB‑YefM toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,089,1980TA 6.1% 369.0 / 8.6 127intergenic (‑70/+13)yeeZ/yoeBputative epimerase/toxin of the YoeB‑YefM toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base T (66/53);  new base A (6/2);  total (72/55)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 139 total aligned reads displayed.

CTTGCCAGCCTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCTGGCTTCCCCTTACCT  >  NC_000913/2089050‑2089346
                                                                                                                                                    |                                                                                                                                                    
cTTGCCAGCCTCGCGCTGAAAGTGACATCGCCCGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTatg                                                                                                                                                     <  2:462734/150‑1 (MQ=255)
cTTGCCAGCCTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTatg                                                                                                                                                     <  1:576396/150‑1 (MQ=255)
cTTGCCAGCCTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTatg                                                                                                                                                     <  2:778580/150‑1 (MQ=255)
cTTGCCAGCCTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTatg                                                                                                                                                     <  2:1572029/150‑1 (MQ=255)
   gCCAGCCTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGCCCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTatgaat                                                                                                                                                  >  1:21127/1‑150 (MQ=255)
   gCCAGCCTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTatgaat                                                                                                                                                  >  1:635937/1‑150 (MQ=255)
   gCCAGCCTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTatgaat                                                                                                                                                  >  2:655361/1‑150 (MQ=255)
    ccAGCCTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTAATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTatgaatg                                                                                                                                                 <  1:146408/150‑1 (MQ=255)
      aGCCTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAtgaatgaa                                                                                                                                               >  1:186699/1‑150 (MQ=255)
         cTCGCGCTGAAAGTGACAGCGCCAGCGGCATGCCTAACCACCGTAACCCGATAATTGAGAACTTTTTCATCCGTTCTCTCCTGATACACCTTATATCTATAAGGCTACGCTATCGTATCAAAACTGACAATTCTTTCATAgaatgaattg                                                                                                                                            <  2:347106/150‑3 (MQ=255)
           cgcgCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCAGACATTTGCGACGTATTTCATCCGATATCTCCTGATCCACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAAGCTGt                                                                                                                                          >  1:93476/1‑150 (MQ=255)
              gcTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCa                                                                                                                                       <  1:12769/150‑1 (MQ=255)
                   aaGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAACACTGACAATTCATTCTATGAATGAATCTGTTCaattat                                                                                                                                  >  1:1642181/1‑150 (MQ=255)
                    aGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTGTTTCATCCGTTATCTCCTGATACACCTTTGATCTATAAGGTTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCaataatt                                                                                                                                 >  2:94575/1‑149 (MQ=255)
                     gggACATCGCCAGCGGCATGCCTAACCACCCTAACCCGTCAATTGCGACGTTTTTCATCCGTTATCTCCTGATTCACCTTAGATATATTAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATTATGa                                                                                                                                <  2:918279/148‑1 (MQ=255)
                     gTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGa                                                                                                                                >  2:804523/1‑150 (MQ=255)
                     gTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGa                                                                                                                                <  1:417900/150‑1 (MQ=255)
                       gACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATa                                                                                                                              <  2:1188489/150‑1 (MQ=255)
                       gACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGATAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATa                                                                                                                              >  1:514380/1‑150 (MQ=255)
                          aTCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAATGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACg                                                                                                                           >  2:386770/1‑150 (MQ=255)
                          aTCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACg                                                                                                                           <  2:814048/150‑1 (MQ=255)
                          aTCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACATTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTAATTCTATGAATGAATCTGTTCAATAATGATTACg                                                                                                                           >  2:1505205/1‑150 (MQ=255)
                            cGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGAc                                                                                                                         >  2:115510/1‑150 (MQ=255)
                             gCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACa                                                                                                                        >  1:1387549/1‑150 (MQ=255)
                             gCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACa                                                                                                                        <  2:1270030/150‑1 (MQ=255)
                             gCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGAGCTAGAAGGCTACGCGAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACa                                                                                                                        >  1:1008153/1‑150 (MQ=255)
                              ccAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGCCCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACAt                                                                                                                       >  2:1567395/1‑150 (MQ=255)
                                aGCGGCATGCCTAACCACCCTAACCAGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACTtgc                                                                                                                     <  1:33589/150‑1 (MQ=255)
                                 gCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACAtgct                                                                                                                    >  2:195413/1‑150 (MQ=255)
                                  cGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACAtgctg                                                                                                                   <  1:1005451/150‑1 (MQ=255)
                                    gCATGCCTACCCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCa                                                                                                                 >  2:142026/1‑150 (MQ=255)
                                    gCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCa                                                                                                                 >  2:237252/1‑150 (MQ=255)
                                             accaccCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGt                                                                                                        <  2:1544164/150‑1 (MQ=255)
                                              ccaccCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGAAAACGACATGCTGCAATGAGCAGTg                                                                                                       <  2:614901/150‑1 (MQ=255)
                                               caccCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGa                                                                                                      >  1:786021/1‑150 (MQ=255)
                                                 cccTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCAGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAAt                                                                                                    <  2:471056/150‑1 (MQ=255)
                                                     aaCCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAATCTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAAtcgtc                                                                                                >  2:25371/1‑150 (MQ=255)
                                                       cccGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCgg                                                                                              <  1:1388029/150‑1 (MQ=255)
                                                       cccGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCgg                                                                                              <  1:716535/150‑1 (MQ=255)
                                                       cccGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCgg                                                                                              >  2:587457/1‑150 (MQ=255)
                                                       cccGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCgg                                                                                              >  1:1097596/1‑150 (MQ=255)
                                                        ccGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGAAGTGAATCGTCGGt                                                                                             >  2:719237/1‑150 (MQ=255)
                                                        ccGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCATAGATCTATAAGGCTACGCTAGCGTATCAAAACTCACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTTCAATGAGCAGTGAATCGTCGGt                                                                                             >  1:300668/1‑150 (MQ=255)
                                                               ttGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACcgcg                                                                                      >  2:897173/1‑150 (MQ=255)
                                                                tGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGt                                                                                     <  2:1638867/150‑1 (MQ=255)
                                                                 gCGACCTTTTTAATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTAAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCTta                                                                                    >  1:216608/1‑150 (MQ=255)
                                                                     ccTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATAcc                                                                                <  1:22944/150‑1 (MQ=255)
                                                                     ccTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATAcc                                                                                <  2:442373/150‑1 (MQ=255)
                                                                      cTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCa                                                                               >  2:1653752/1‑150 (MQ=255)
                                                                        ttttCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGa                                                                             >  1:163861/1‑150 (MQ=255)
                                                                              tCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGc                                                                       <  1:836155/150‑1 (MQ=255)
                                                                               ccGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGct                                                                      >  2:1320302/1‑150 (MQ=255)
                                                                               acGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGct                                                                      <  1:1431161/149‑1 (MQ=255)
                                                                                     tctcCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATCACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTg                                                                <  2:154555/150‑1 (MQ=255)
                                                                                      ctcCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGt                                                               >  2:1538783/1‑150 (MQ=255)
                                                                                          tGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATg                                                           <  1:878084/150‑1 (MQ=255)
                                                                                          tGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTAACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGTTGCTCATCTGTAATg                                                           >  1:158518/1‑150 (MQ=255)
                                                                                            aTACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCg                                                         <  1:587457/150‑1 (MQ=255)
                                                                                                acCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCgg                                                     <  1:1648354/150‑1 (MQ=255)
                                                                                                 ccTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCggg                                                    <  2:977415/150‑1 (MQ=255)
                                                                                                    tAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCTTATACCAGACGGTGCTCCTCTGTAATGAGTCGGGAcc                                                 <  1:1001253/150‑1 (MQ=255)
                                                                                                      gATCTATAAGGATACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAg                                               >  2:421164/1‑150 (MQ=255)
                                                                                                       aTCTATAAGGCTACGCTAGAGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGa                                              <  2:805225/150‑1 (MQ=255)
                                                                                                         cTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGCCGGTGCTCCTCTGTAATGCGTCGGGACCAGaaa                                            >  2:1041244/1‑150 (MQ=255)
                                                                                                         cTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGaaa                                            >  2:1124245/1‑150 (MQ=255)
                                                                                                              aGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGa                                       >  2:1037531/1‑150 (MQ=255)
                                                                                                               ggCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGAc                                      >  1:6761/1‑150 (MQ=255)
                                                                                                               ggCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGAc                                      <  2:1177570/150‑1 (MQ=255)
                                                                                                               ggCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGAc                                      >  1:604317/1‑150 (MQ=255)
                                                                                                               ggCTACGCTAGCGTATCAAAACTGACAATACATTCTATGAATGAATCTGTTCAATAATGAAAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGAc                                      >  1:598218/1‑150 (MQ=255)
                                                                                                                    cGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAAtt                                 >  2:1326799/1‑150 (MQ=255)
                                                                                                                      cTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATTATTATAACGCCATGCTGCAATGAGCAGTGAATCGTCGGTAACCGAGTATACCAGACGGTGCTCCTCTGTATTGCGTCGGGACCACAAACCTGGCAAAGTAt                               >  2:444595/1‑150 (MQ=255)
                                                                                                                        aGCTTATCAAAACTGACAATTCTTTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCTGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGt                             <  1:930393/150‑1 (MQ=255)
                                                                                                                             aTCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAg                        >  2:64653/1‑150 (MQ=255)
                                                                                                                             aTCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAg                        <  1:386770/150‑1 (MQ=255)
                                                                                                                             aTCAAAACTGACAATTCATTCTAAGAAAGAATCTGTTCAATAATGATAACGACATGCTCCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAg                        <  1:393077/150‑1 (MQ=255)
                                                                                                                               cAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAggg                      <  2:403843/150‑1 (MQ=255)
                                                                                                                               cAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAggg                      >  2:410185/1‑150 (MQ=255)
                                                                                                                                aaaaCTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAgggg                     >  1:432760/1‑150 (MQ=255)
                                                                                                                                aaaaCTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAgggg                     >  2:1486798/1‑150 (MQ=255)
                                                                                                                                   aCTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACTACATGCTGCAATGAGCAGTGAATCTTCGGTAACATCGTATACCAGACGGTGATCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTc                  >  1:78782/1‑150 (MQ=255)
                                                                                                                                   aCTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGtta                  >  1:446656/1‑149 (MQ=255)
                                                                                                                                    cTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCt                 <  1:237252/150‑1 (MQ=255)
                                                                                                                                     aacagattcattcATAGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCTg                >  1:1199303/6‑150 (MQ=255)
                                                                                                                                      gACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCTgg               <  2:601212/150‑1 (MQ=255)
                                                                                                                                       cagattcattcATAGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGTTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGTAACCTGACAAATTTTGTTTCAGGGGTTCTGGc              >  2:680006/4‑150 (MQ=255)
                                                                                                                                       cagattcattcATAGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGTTTCTGGc              >  2:1681045/4‑150 (MQ=255)
                                                                                                                                       cagattcattcATAGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCTGGc              >  1:922740/4‑150 (MQ=255)
                                                                                                                                       cagattcattcATAGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCTGGc              >  1:1509792/4‑150 (MQ=255)
                                                                                                                                         aaTTCAATCTATGAATGAATCTGTTCAAAAATGATAACGACATGCGGCAATGAGCAGTGAATCGTCGGTAACCGCGAATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATAATGTTTCAGGGGTTCTGGCtt            <  2:230130/150‑1 (MQ=255)
                                                                                                                                           ttcattcTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCTGGCTTcc          <  2:186699/150‑1 (MQ=255)
                                                                                                                                           ttcattcTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGACGGGACCAGAAACCTGACTAATTATGTTTCAGGGGTTCTGGCGTcc          >  1:169259/1‑150 (MQ=255)
                                                                                                                                             cattcTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCTGGCTTcccc        <  1:1037531/150‑1 (MQ=255)
                                                                                                                                             cattcTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAAGCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCTGGCTTcccc        <  2:482969/150‑1 (MQ=255)
                                                                                                                                               ttcTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCTGGCTTCCCCtt      <  2:1248412/150‑1 (MQ=255)
                                                                                                                                                tcTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCTGGCTTCCCCTTa     >  1:105253/1‑150 (MQ=255)
                                                                                                                                                tcTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTATGGCTTCCCCTTa     >  1:1395542/1‑150 (MQ=255)
                                                                                                                                                 cTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCTGGCTTCCCCTTAc    >  2:1163265/1‑150 (MQ=255)
                                                                                                                                                 cTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCTGGCTTCCCCTTAc    >  2:1266578/1‑150 (MQ=255)
                                                                                                                                                   atgaatgaATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCTGGCTTCCCCTTTCCt  >  1:692015/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTTGCCAGCCTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTATGAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCTGGCTTCCCCTTACCT  >  NC_000913/2089050‑2089346

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: