Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,439,588 G→C 9.5% intergenic (‑79/+28) yhdN ← / ← rplQ DUF1992 family protein/50S ribosomal subunit protein L17

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,439,5880GC 9.5% 409.3 / 27.2 135intergenic (‑79/+28)yhdN/rplQDUF1992 family protein/50S ribosomal subunit protein L17
Reads supporting (aligned to +/- strand):  ref base G (57/65);  new base C (9/4);  total (66/69)
Fisher's exact test for biased strand distribution p-value = 1.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.08e-01

Only 100 of 143 total aligned reads displayed.

AGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCG  >  NC_000913/3439439‑3439737
                                                                                                                                                        |                                                                                                                                                     
aGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACccgccg                                                                                                                                                       >  2:1044750/1‑150 (MQ=255)
aGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACccgccg                                                                                                                                                       >  2:390844/1‑150 (MQ=255)
 ggTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCgg                                                                                                                                                      >  2:940448/1‑150 (MQ=255)
  gTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCggg                                                                                                                                                     >  2:423749/1‑150 (MQ=255)
  gTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑‑AAA‑AAACCCGCCgggg                                                                                                                                                    <  2:1704666/150‑1 (MQ=255)
   ttATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCgggg                                                                                                                                                    >  2:1368268/1‑150 (MQ=255)
    tATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGc                                                                                                                                                   <  1:717548/150‑1 (MQ=255)
    tATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGc                                                                                                                                                   <  1:1276312/150‑1 (MQ=255)
    tATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGc                                                                                                                                                   <  1:1578759/150‑1 (MQ=255)
    tATCAAACTCACCTTTCGCTAGCGCTGCTGCTATATGCCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGc                                                                                                                                                   <  2:1055196/150‑1 (MQ=255)
     aTCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCg                                                                                                                                                  >  2:362293/1‑150 (MQ=255)
       cAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCggg                                                                                                                                                <  2:1220481/150‑1 (MQ=255)
       cAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCggg                                                                                                                                                <  2:896025/150‑1 (MQ=255)
       cAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCAGGGGCggg                                                                                                                                                >  1:1502678/1‑150 (MQ=255)
       aaaaTTCAGCTTGCGCTTGCGCTTCTGCTATATCGCGCTCGGCCCAGTGGTAAATAAACAACATAAATCCTCCCGGGCTTGACCAAAACGTGTACAGTTATTTTAGTAAAGTGGGTATATTAGGGTTAT‑‑AAAA‑AAACCCGCCGGGGCggg                                                                                                                                                <  2:469433/149‑1 (MQ=255)
         aaCTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGtt                                                                                                                                              >  1:1015477/1‑150 (MQ=255)
          aCTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGttt                                                                                                                                             <  1:512545/150‑1 (MQ=255)
             cACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGtttttt                                                                                                                                          <  1:1359823/150‑1 (MQ=255)
             cACCTTTCGCTAGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAA‑‑AAAA‑AAACCCGCCGGGGCGGGtttttt                                                                                                                                          <  1:48707/150‑1 (MQ=255)
              aCCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCCCGGCGGGttttttt                                                                                                                                         <  1:816101/150‑1 (MQ=255)
                cattCGCGTGCGCTTCTGCTATCTGGCGCTCTGCCCACTGGGCAAGTAACCACATAAATACTCCCGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAA‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTAc                                                                                                                                       <  1:1408071/148‑1 (MQ=255)
                cTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAACTACTCCAGGGGCTGAACAAAACGATTACAGCTATTGTAGATAAGTGTGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCTGGTTTTTTTAc                                                                                                                                       <  1:1401063/150‑1 (MQ=255)
                         gcgcTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCa                                                                                                                              <  1:1490042/150‑1 (MQ=255)
                               cTGCTATATGGCGCTCAGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATCAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTAc                                                                                                                        <  1:604066/150‑1 (MQ=255)
                               cTGATATATGGCGCTCTGCCCACTGGTCAAGGAACCACATAAATACTCCAGGGGATGAACAAAAAGAGGACAGCTATTGTAGATAAGTGGGGATCCGACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGGTGTTGGTCGTTGCTTCAGATGAc                                                                                                                        >  1:131105/1‑150 (MQ=255)
                                tGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGTGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTAct                                                                                                                       <  2:195768/150‑1 (MQ=255)
                                tGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACAAAAATAATCCAGGGGATTAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTAct                                                                                                                       <  2:210282/150‑1 (MQ=255)
                                 gCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTActc                                                                                                                      >  1:541966/1‑150 (MQ=255)
                                 gCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTActc                                                                                                                      >  1:1177504/1‑150 (MQ=255)
                                 gCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTAAGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGATTTTTTACGTTGCTTCAGATGActc                                                                                                                      >  1:1331738/1‑150 (MQ=255)
                                     tatGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTgca                                                                                                                  >  1:1242394/1‑150 (MQ=255)
                                       tGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTgcagc                                                                                                                >  2:769926/1‑150 (MQ=255)
                                        ggCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGTTGAACAAAAAGAGTACAGCTATTGTAGATAAGGGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTgcagca                                                                                                               <  1:875364/150‑1 (MQ=255)
                                           gcTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTGCAGAGTAATCTgaagcaacg                                                                                                            <  1:1044750/150‑4 (MQ=255)
                                                 cccACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGGTTTTTTACGTTGCGTCAGATTACTCTGCAGCAGCTTCTGCt                                                                                                      >  2:1120553/1‑150 (MQ=255)
                                                      tGGTCAAGTAACCACATAAATACTCCCGGGGATGAACAAAAAGAGTCCAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTtctc                                                                                                 <  2:892716/150‑1 (MQ=255)
                                                      tGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTtctc                                                                                                 <  1:736797/150‑1 (MQ=255)
                                                       ggTCAAGTAACCACATAAATATTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTctct                                                                                                >  1:1102210/1‑150 (MQ=255)
                                                       ggTCAAGTAACAACATAAATACTCCAGGGGATTAACAAAAATAGTACAGCTATTGTAGATAATTGTGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTctct                                                                                                <  2:270508/150‑1 (MQ=255)
                                                        gTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTg                                                                                               <  2:1527527/150‑1 (MQ=255)
                                                          cAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGaa                                                                                             <  1:643742/150‑1 (MQ=255)
                                                            aGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACg                                                                                           >  1:1144758/1‑150 (MQ=255)
                                                                aCCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCa                                                                                       <  2:1271916/150‑1 (MQ=255)
                                                                   acaTAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAAcc                                                                                    <  2:1618675/150‑1 (MQ=255)
                                                                   acaTAAATACTCAAGGGGATTAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCTGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAAcc                                                                                    <  2:1316810/150‑1 (MQ=255)
                                                                    caTAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCa                                                                                   >  2:451941/1‑150 (MQ=255)
                                                                      tAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCCGc                                                                                 >  2:159130/1‑150 (MQ=255)
                                                                             tCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGt                                                                          <  1:1249933/150‑1 (MQ=255)
                                                                              ccAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTa                                                                         >  2:731420/1‑150 (MQ=255)
                                                                                aGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTAGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAg                                                                       <  2:1144384/150‑1 (MQ=255)
                                                                                 ggggATGAACATAAAGAGTACAGTTTTTGTAGTTTAGTTGGGTTACTACGGGTTT‑‑AAAATAACCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTAATCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAg                                                                       <  1:1530935/150‑1 (MQ=255)
                                                                                       gAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCgg                                                                >  2:1383140/1‑150 (MQ=255)
                                                                                       gAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCgg                                                                >  2:1348536/1‑150 (MQ=255)
                                                                                         aCAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCgg                                                                >  1:1555063/1‑150 (MQ=255)
                                                                                         aCAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGgcg                                                              >  2:461681/1‑150 (MQ=255)
                                                                                          cAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGgcgc                                                             >  1:1176575/1‑150 (MQ=255)
                                                                                           aaaaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGcgcg                                                            >  1:724320/1‑150 (MQ=255)
                                                                                             aaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGtt                                                          >  2:952705/1‑150 (MQ=255)
                                                                                              aaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AATCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTg                                                         <  1:626052/150‑1 (MQ=255)
                                                                                                    aCAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCt                                                   >  2:565704/1‑150 (MQ=255)
                                                                                                     cAGCTTTTGCAGATAAGTGGGGTTACTACGGGTAT‑‑AAAA‑TAACCGGGCGGGGCGGGTTTTTTTACGTTGGTTCAGATTTCTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCATCGTGCTCGATGTAAGCCATCGGCGCGTTGTCGCGTg                                                  <  2:1603756/150‑1 (MQ=255)
                                                                                                     cAGCTATTGTAGTTAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTg                                                  <  1:220535/150‑1 (MQ=255)
                                                                                                      aGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGc                                                 >  1:1153104/1‑150 (MQ=255)
                                                                                                      aGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGc                                                 <  2:1233819/150‑1 (MQ=255)
                                                                                                      aGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGc                                                 <  1:362293/150‑1 (MQ=255)
                                                                                                       gCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCa                                                >  2:1390779/1‑150 (MQ=255)
                                                                                                         tATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACg                                              >  1:1038416/1‑150 (MQ=255)
                                                                                                            tGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAAGCAGATCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGaa                                           <  2:848555/150‑1 (MQ=255)
                                                                                                             gTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTGTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAg                                          >  1:337004/1‑150 (MQ=255)
                                                                                                               aGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTACTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGcc                                        >  1:1673530/1‑150 (MQ=255)
                                                                                                                gATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGccc                                       >  2:607939/1‑149 (MQ=255)
                                                                                                                 aTAAGTGGGGATACTACGGGTAT‑‑AAAA‑ACACCCGCCGGGGCGGGTTTTTTTACGTTGCGTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCcac                                      <  1:295060/150‑1 (MQ=255)
                                                                                                                 aTAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGcccc                                      >  2:860788/1‑148 (MQ=255)
                                                                                                                    aGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACt                                   >  2:404418/1‑150 (MQ=255)
                                                                                                                     gTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACtt                                  >  2:1382628/1‑150 (MQ=255)
                                                                                                                      gggggATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTc                                 <  2:1456444/149‑1 (MQ=255)
                                                                                                                        gggATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAg                               <  2:296480/150‑1 (MQ=255)
                                                                                                                         ggATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACCCTTCAGa                              >  2:1329916/1‑150 (MQ=255)
                                                                                                                          gATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGaa                             <  1:1383140/150‑1 (MQ=255)
                                                                                                                          gATACTACGGGTAA‑‑AAAA‑ATACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGaa                             <  2:725795/150‑1 (MQ=255)
                                                                                                                            tactacGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATa                           <  2:1102210/150‑1 (MQ=255)
                                                                                                                                ccGGCTAT‑‑ACAA‑AAACCCGCCGGGGCGGGTTCTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCTCCTGCACGGAAGCCACACTTCAGAATACGAg                       <  2:1134604/149‑1 (MQ=255)
                                                                                                                                acGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTTCTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCCTCGGCGCGTTGTCGCCTGCACGGAAGCCGCACTTCAGAATACGAg                       >  1:94311/1‑150 (MQ=255)
                                                                                                                                  gggTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAgtg                     <  1:584597/150‑1 (MQ=255)
                                                                                                                                   aacgt‑‑aAAA‑AAACCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCGTCTGCTTTCTCTGAACGCTCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGGCTGCACGGAAGCCACACTTCTGAATACGAgtgt                    >  1:1467142/5‑150 (MQ=255)
                                                                                                                                    acgt‑‑aAAA‑AAACCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTTCACGGAAGCCACACTTCAGAATACGAGTGTa                   >  1:39934/4‑150 (MQ=255)
                                                                                                                                    acgt‑‑aAAA‑AAACCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTa                   >  1:1714578/4‑150 (MQ=255)
                                                                                                                                       t‑‑aAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTCGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAacc                >  2:1352117/1‑150 (MQ=255)
                                                                                                                                          aaaa‑aaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGACCGATCATCCAGCTCGATGTAAGCCATCGGCGCGTTGTCGACTGCAAGGAAGCCACCCTGCAGAATAAGAGTGTCaccc               >  1:1667051/1‑149 (MQ=255)
                                                                                                                                          aaaa‑aaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAacca               >  2:1169534/1‑150 (MQ=255)
                                                                                                                                          aaaa‑aaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAacca               >  2:527130/1‑150 (MQ=255)
                                                                                                                                           aaa‑aaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCCTCGGCGCGTTGTCGCCTGCACGGAAGCCCCACTGCAGAATACGAGTGTAacccc              >  2:1578403/1‑148 (MQ=255)
                                                                                                                                           aaa‑aaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAacctc              >  1:11222/1‑148 (MQ=255)
                                                                                                                                           aaa‑aaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAaccac              >  2:789497/1‑150 (MQ=255)
                                                                                                                                           aaa‑aaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAaccac              >  2:812270/1‑150 (MQ=255)
                                                                                                                                            aa‑aaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAGCGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGGAcccacc             >  1:265487/1‑150 (MQ=255)
                                                                                                                                               aaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCgg           <  1:99479/150‑1 (MQ=255)
                                                                                                                                                  cccGCCGTGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCTATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCAc        <  1:291463/150‑1 (MQ=255)
                                                                                                                                                  cccGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCAc        >  1:434348/1‑150 (MQ=255)
                                                                                                                                                        ggggCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCACCCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGACGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTcg  >  2:1698595/1‑150 (MQ=255)
                                                                                                                                                        |                                                                                                                                                     
AGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAT‑‑AAAA‑AAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCG  >  NC_000913/3439439‑3439737

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