Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,517,057 C→A 10.0% intergenic (+55/+42) patD → / ← yncL gamma‑aminobutyraldehyde dehydrogenase/uncharacterized protein YncL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,517,0570CA10.0% 347.8 / 21.7 116intergenic (+55/+42)patD/yncLgamma‑aminobutyraldehyde dehydrogenase/uncharacterized protein YncL
Reads supporting (aligned to +/- strand):  ref base C (58/46);  new base A (5/7);  total (63/53)
Fisher's exact test for biased strand distribution p-value = 3.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.52e-01

Only 100 of 117 total aligned reads displayed.

CGGTGGGCAGAAACTTTCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGAAACATTCATTTTTGCCACCTT  >  NC_000913/1516909‑1517206
                                                                                                                                                    |                                                                                                                                                     
cGGTGGGCAGAAACTTTCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACt                                                                                                                                                      <  2:1440207/150‑1 (MQ=255)
 ggTGGGCAGAAACTTTCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTg                                                                                                                                                     <  2:566009/150‑1 (MQ=255)
   tGGGCAGAAACTTTCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGAt                                                                                                                                                   >  2:1308148/1‑150 (MQ=255)
   tGGGCAGAAACTTTCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGAt                                                                                                                                                   <  2:543514/150‑1 (MQ=255)
    gggCAGAAACTTTCTGGTTACGGGAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGAtt                                                                                                                                                  <  1:1318014/150‑1 (MQ=255)
    gggCAGAAACTTTCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGAtt                                                                                                                                                  >  2:1114255/1‑150 (MQ=255)
         gAAACTTTCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTctct                                                                                                                                             >  1:576706/1‑150 (MQ=255)
         gAAACTTTCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTctct                                                                                                                                             >  2:582275/1‑150 (MQ=255)
          aaaCTTTCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAATGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTg                                                                                                                                            <  2:127265/150‑1 (MQ=255)
          aaaCTTTCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTg                                                                                                                                            <  1:719658/150‑1 (MQ=255)
          aaaCTTTCTGGTTACGGCAAGGATATGGCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTg                                                                                                                                            <  2:1057789/150‑1 (MQ=255)
             cTTTCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGTTTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACa                                                                                                                                         >  2:1528016/1‑150 (MQ=255)
             cTTTCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGAAATCAGTGCTCTGACa                                                                                                                                         >  2:1051623/1‑150 (MQ=255)
              tttCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGAAATCAGTGCTCTGACAg                                                                                                                                        <  1:949334/150‑1 (MQ=255)
              tttCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGAAATCAGTGCTCTGACAg                                                                                                                                        <  1:10787/150‑1 (MQ=255)
                tCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTa                                                                                                                                      <  2:727649/150‑1 (MQ=255)
                 cTGGTTACGCCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGAAATCAGTGCTCTGACAGTAc                                                                                                                                     <  1:1234116/150‑1 (MQ=255)
                  tGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGAAATCAGTGCTCTGACAGTACg                                                                                                                                    <  2:815256/150‑1 (MQ=255)
                    gTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGAAATCAGTGCTCTGACAGTACgca                                                                                                                                  <  2:840903/150‑3 (MQ=39)
                     ttACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGtgc                                                                                                                                 <  2:1039820/150‑1 (MQ=255)
                       aCGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGtgctg                                                                                                                               >  1:183301/1‑150 (MQ=255)
                        cGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGAAATCAGTGCTCTGACAGTACgcaagat                                                                                                                              <  1:1051623/150‑7 (MQ=38)
                                atatGGCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCaaaa                                                                                                                      <  2:1312713/150‑1 (MQ=255)
                                    gTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGAAATCAGTGCTCTGACAGt                                                                                                                                       <  2:284382/129‑1 (MQ=39)
                                    gTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGAAATCAGTGCTCTGACAGt                                                                                                                                       >  1:284382/1‑129 (MQ=39)
                                       aCTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCTTACTGTCAGAGCAATGATTTCTATGACAGTACGTGCTGCGTGCAAAATCACGAc                                                                                                               >  2:470574/1‑150 (MQ=255)
                                             tGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCaaa                                                                                                         >  1:771674/1‑150 (MQ=255)
                                             tGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCaaa                                                                                                         >  1:1347734/1‑150 (MQ=255)
                                             tGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCaaa                                                                                                         >  2:1409609/1‑150 (MQ=255)
                                               ggCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATa                                                                                                       >  2:56010/1‑150 (MQ=255)
                                               ggCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATa                                                                                                       >  1:485222/1‑150 (MQ=255)
                                               ggCTGGAGGATTACACCGTAGTCCCCCACGTCATGTTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTAAGAGATGATCATGCGTACTGTCACCGCACTGATTTCTCTTACAGTACGTGCTGCGTGCAAAATCACTACATCAAATa                                                                                                       >  2:528947/1‑150 (MQ=255)
                                                gCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATaa                                                                                                      <  2:157956/150‑1 (MQ=255)
                                                   ggaggaTTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAAcc                                                                                                   >  2:69977/1‑150 (MQ=255)
                                                   ggaggaTTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAAcc                                                                                                   >  2:1158067/1‑150 (MQ=255)
                                                   ggaggaTTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAAcc                                                                                                   >  2:1147068/1‑150 (MQ=255)
                                                   ggaggaTTACACCGTCGTCAGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATAAAATAAAcc                                                                                                   >  2:633587/1‑150 (MQ=255)
                                                    gaggaTTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCa                                                                                                  >  2:456744/1‑150 (MQ=255)
                                                     aggaTTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAg                                                                                                 <  1:678678/150‑1 (MQ=255)
                                                     aggaTTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAg                                                                                                 <  2:576706/150‑1 (MQ=255)
                                                      ggaTTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGc                                                                                                >  2:1221483/1‑150 (MQ=255)
                                                       gaTTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGcc                                                                                               >  2:715768/1‑150 (MQ=255)
                                                            cacCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAAtc                                                                                          >  2:1028753/1‑150 (MQ=255)
                                                             acCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAAtct                                                                                         >  1:884490/1‑150 (MQ=255)
                                                               cgtcgtCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCa                                                                                       <  1:437895/150‑1 (MQ=255)
                                                                  cgtcCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATag                                                                                    <  2:669715/150‑1 (MQ=255)
                                                                  cgtcCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATag                                                                                    >  1:749264/1‑150 (MQ=255)
                                                                    tcCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATagag                                                                                  <  2:1215237/150‑1 (MQ=255)
                                                                     ccGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAgaga                                                                                 >  1:32827/1‑150 (MQ=255)
                                                                      cGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGAt                                                                                >  1:1123554/1‑150 (MQ=255)
                                                                      cGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTATCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGAt                                                                                >  1:1098868/1‑150 (MQ=255)
                                                                           cGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAagag                                                                           <  2:48408/150‑1 (MQ=255)
                                                                             tCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAagagca                                                                         <  2:474838/150‑3 (MQ=255)
                                                                             tCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTa                                                                         <  1:355600/150‑1 (MQ=255)
                                                                                 ggTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTaaaca                                                                     <  2:1347734/150‑1 (MQ=255)
                                                                                  gTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTaaacaa                                                                    >  2:197827/1‑150 (MQ=255)
                                                                                   ttAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTaaacaaa                                                                   <  1:69977/150‑1 (MQ=255)
                                                                                   ttAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTaaacaaa                                                                   <  2:183301/150‑1 (MQ=255)
                                                                                    gtaataTTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTaaacaaac                                                                  >  2:861991/6‑150 (MQ=255)
                                                                                     aaaCATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAAcc                                                                 >  1:229994/1‑150 (MQ=255)
                                                                                      aaCATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCa                                                                <  1:582275/150‑1 (MQ=255)
                                                                                          ttAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAagag                                                            >  2:1359459/1‑150 (MQ=255)
                                                                                          ttAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAAcag                                                            >  1:1009020/1‑150 (MQ=255)
                                                                                            aaGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAAcagca                                                          <  1:715768/150‑1 (MQ=255)
                                                                                            aaGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAAcagca                                                          <  1:456744/150‑1 (MQ=255)
                                                                                             aGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTaaac                                                                      >  2:1388828/1‑137 (MQ=255)
                                                                                             aGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTaaac                                                                      <  1:1388828/137‑1 (MQ=255)
                                                                                              ggATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACagcagc                                                        >  1:992460/1‑150 (MQ=255)
                                                                                              ggATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACagcagc                                                        >  1:1294782/1‑150 (MQ=255)
                                                                                              ggATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACagcagc                                                        >  1:394599/1‑150 (MQ=255)
                                                                                               gATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAgcagca                                                       >  1:111885/1‑150 (MQ=255)
                                                                                                       tGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAAtt                                               <  2:657940/150‑1 (MQ=255)
                                                                                                        gCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAAttt                                              <  1:188238/150‑1 (MQ=255)
                                                                                                        gCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAAttt                                              >  1:632414/1‑150 (MQ=255)
                                                                                                        gCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGACCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAAttt                                              <  2:1264789/150‑1 (MQ=255)
                                                                                                         cAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTa                                             >  2:505711/1‑150 (MQ=255)
                                                                                                           gATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATc                                           <  2:1169987/150‑1 (MQ=255)
                                                                                                           gATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATc                                           >  1:1255139/1‑150 (MQ=255)
                                                                                                                  agagTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGAGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATAAGaacta                                    >  2:779616/1‑150 (MQ=255)
                                                                                                                   gagTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGtactac                                   <  1:1503722/150‑1 (MQ=255)
                                                                                                                        tGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTc                              <  1:1038527/150‑1 (MQ=255)
                                                                                                                          agaTGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTa                            >  1:966443/1‑150 (MQ=255)
                                                                                                                           gatgatCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTAc                           <  2:1102170/150‑1 (MQ=255)
                                                                                                                             tgatCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTg                         >  2:182351/1‑150 (MQ=255)
                                                                                                                             tgatCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTg                         <  1:616272/150‑1 (MQ=255)
                                                                                                                                tCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGaa                      <  1:1528016/150‑1 (MQ=255)
                                                                                                                                  ttGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGAAAc                    <  2:1273014/150‑1 (MQ=255)
                                                                                                                                  ttGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGAAAc                    >  2:194948/1‑150 (MQ=255)
                                                                                                                                   tGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCCGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCCGCAAAGCAATTTATCAGTACTACAGTTCTACTGGAAAca                   >  1:840195/1‑150 (MQ=255)
                                                                                                                                      gTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGAAAcattc                <  1:1221483/150‑1 (MQ=255)
                                                                                                                                       tACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGAAAcattca               >  1:1290229/1‑150 (MQ=255)
                                                                                                                                        aCTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGAAAcattcat              >  2:914675/1‑150 (MQ=255)
                                                                                                                                         cTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGAAAcattcatt             >  1:207124/1‑150 (MQ=255)
                                                                                                                                         cTGTCAGAGAAATCAGTGCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGAAAcattcatt             >  1:1239260/1‑150 (MQ=255)
                                                                                                                                           gTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGAAACATTCAtttt           >  1:1333931/1‑150 (MQ=255)
                                                                                                                                             cAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGAAACATTCATTTTTg         >  1:981400/1‑150 (MQ=255)
                                                                                                                                              agagCACTGATTTCTCTGACAGTACGTGCTTCGTGCAAAATCACTAAATCAAATAAACCAGCCAACTCTCATAGAGATAAGATTAACCAAACCAACAGCAGCAAAGAAATTTATCATTACTACAGTTCTACTGGAAACATTCATTTTTGc        <  1:516876/150‑1 (MQ=255)
                                                                                                                                                   aCTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTc                              >  2:1079406/1‑123 (MQ=255)
                                                                                                                                                   aCTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTc                              <  1:1079406/123‑1 (MQ=255)
                                                                                                                                                    cTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGAAACATTCATTTTTGCCACCtt  >  1:797935/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGGTGGGCAGAAACTTTCTGGTTACGGCAAGGATATGTCACTTTATGGGCTGGAGGATTACACCGTCGTCCGCCACGTCATGGTTAAACATTAAGGATAATATTGCAGATCGTAAGAGTATGAGATGATCTTGCGTACTGTCAGAGCACTGATTTCTCTGACAGTACGTGCTGCGTGCAAAATCACTACATCAAATAAACCAGCCAAATCTCATAGAGATAAGAGTAAACAAACCAACAGCAGCAAAGAAATTTATCAGTACTACAGTTCTACTGGAAACATTCATTTTTGCCACCTT  >  NC_000913/1516909‑1517206

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