Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 729,377 A→T 8.1% intergenic (+37/‑206) ybfA → / → rhsC DUF2517 family protein/Rhs protein with putative toxin domain; putative neighboring cell growth inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913729,3770AT 8.1% 295.3 / 11.8 107intergenic (+37/‑206)ybfA/rhsCDUF2517 family protein/Rhs protein with putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (48/49);  new base T (2/7);  total (51/56)
Fisher's exact test for biased strand distribution p-value = 1.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.02e-01

Only 100 of 113 total aligned reads displayed.

TCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAA  >  NC_000913/729228‑729526
                                                                                                                                                     |                                                                                                                                                       
tCTGGAAAGACCGCGACCGCTTTTACAGCTACCTGCATTGCGTATGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACAGCCAGCGGa                                                                                                                                                         <  2:331076/150‑1 (MQ=255)
 cTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCACCAGAGGaa                                                                                                                                                        >  1:1516339/1‑150 (MQ=255)
  tGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGaaa                                                                                                                                                       <  2:1004500/150‑1 (MQ=255)
  tGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGaaa                                                                                                                                                       <  1:1488806/150‑1 (MQ=255)
    gAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACg                                                                                                                                                     >  1:932371/1‑150 (MQ=255)
        gACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTGTTGATTGCTGAGTTTTGAATACATTGAAAACCGCCATCCGAAACGCTgg                                                                                                                                                 >  2:369754/1‑150 (MQ=255)
        gACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTgg                                                                                                                                                 >  1:775778/1‑150 (MQ=255)
          ccGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCg                                                                                                                                               <  2:1228295/150‑1 (MQ=255)
           cgcgCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGTTTCCGCTGGCgg                                                                                                                                              <  1:624796/150‑1 (MQ=255)
           cgcgCCCGCTTTTACAGCTACCGGCATCGCGGCTGGGCGAAAACCACCGATAAACCGGTGGGGATGGATCAGGCCGAAAAGGCAACTGCTGATTTTTATTGAATGCTGAGTTTTGAATACATTGAAAACCGCCAGCGTTTCCGCTGGCgg                                                                                                                                              <  2:483267/150‑1 (MQ=255)
            gcgcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAGCCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGt                                                                                                                                             >  2:75167/1‑150 (MQ=255)
            gcgcCCGCTTTTACAGCTACCTGCATCGCGTCTGGGCGAAAACCAGCGATAAACCGGTGTGGAGGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGt                                                                                                                                             <  1:1145721/150‑1 (MQ=255)
              gcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGTTTCCGCTGGCGGttt                                                                                                                                           <  1:1115459/150‑1 (MQ=255)
               cccGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGTTTCCGCTGGCGGtttt                                                                                                                                          <  1:1260901/150‑1 (MQ=255)
               cccGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGTTTCCGCTGGCGGtttt                                                                                                                                          <  2:1508451/150‑1 (MQ=255)
                ccGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGttttt                                                                                                                                         <  1:1311868/150‑1 (MQ=255)
                  gCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCCGCGGAACCGCTGGCGGGTTTTTg                                                                                                                                       >  1:1523373/1‑150 (MQ=255)
                    ttttACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGtt                                                                                                                                     >  1:1516706/1‑149 (MQ=255)
                    tgttACAGCGACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGTTTCCGCTGGCGGttttcaatg                                                                                                                                     <  2:84050/148‑6 (MQ=255)
                       tACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCa                                                                                                                                  >  2:1692547/1‑150 (MQ=255)
                        aCAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGAGCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAg                                                                                                                                 <  1:1447472/150‑1 (MQ=255)
                        aCAGCTACCTGCATCGCGTCCGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCCCCAGCGGAAACGCTGGCGGTTTTTTGTATCAg                                                                                                                                 >  2:1272742/1‑150 (MQ=255)
                                cTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCTACTGGAGATTTTTATTGATTGCTGAGTTTTGATTACATTGAAAAACGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGAAGt                                                                                                                         >  1:1117497/1‑150 (MQ=255)
                                 tGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAAGTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTGTTTTGTATCAGGAAGATGtt                                                                                                                        >  1:262529/1‑150 (MQ=255)
                                  gcttCGCGTGAGGTCGAAAACCTGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGTTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTc                                                                                                                       <  2:1090530/147‑1 (MQ=255)
                                   cATCTCGTATGTTCTAAAACCAGCTATAAACGGGTGTGTTTGTATAAGGCCTAAAAGGCAACTGGGGATTTTTATTGATTGCTTAGTTTTGAATACATTGAAAACCGCAAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCa                                                                                                                      <  2:1517246/150‑1 (MQ=255)
                                      cgcgTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTTAATACATTGAAAAACGCCAGCGGAAACGCTGGCGGTTTTTTGTATCATGAAGATGTTCacag                                                                                                                   >  2:1507282/1‑150 (MQ=255)
                                      cgcgTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCagag                                                                                                                   <  2:915697/150‑1 (MQ=255)
                                      cgcgTCTGGTCGAAAACCAGCGATAAACCGGGGTGGATGGATCAGGCCTAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCagag                                                                                                                   <  1:1702699/150‑1 (MQ=255)
                                              gTCGAAAACCAGCGATAAACCGGTGTGGATGGAGCAGGCCGAAAGGCAACTGGTGATTATTATTGATTGCTGAGTTTTCAATACATTGCAAACCGCCAGCGGAACCGCTGGCGGTTTTTTGTATCAGGAAGATGTACAGAGGAATATTTa                                                                                                          >  1:231379/1‑150 (MQ=255)
                                                cGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAg                                                                                                         <  2:416891/150‑1 (MQ=255)
                                                cGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAAACGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGCAAGATGTTAAGAGGAATATTTAg                                                                                                         >  2:497528/1‑150 (MQ=255)
                                                cGAAAACCAGCGATAAACCGGGGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGCGTTTTGAAAACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAg                                                                                                         <  2:1585195/150‑1 (MQ=255)
                                                 gAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGa                                                                                                        <  2:157847/150‑1 (MQ=255)
                                                 gAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGa                                                                                                        >  2:105054/1‑150 (MQ=255)
                                                  aaaaCCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCTAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGAc                                                                                                       <  2:1187550/150‑1 (MQ=255)
                                                  aaaaCCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGCAACGCTGGCGGTTTTATGTAACACGAAGATGTTCAGAGGAATATTTAGAc                                                                                                       >  1:90258/1‑150 (MQ=255)
                                                     aCCAGCTATAAACCTGTGTGGATGGATCAGGCCGAAAAGTCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATc                                                                                                    <  1:1046311/150‑1 (MQ=255)
                                                     aCCAGCGATAAAGCGGTGTGGATGGAGCAGGCCGAATAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAAAACATTGAAAACCGCCAGCGGAAACGCGGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATc                                                                                                    <  2:856352/150‑1 (MQ=255)
                                                     aCCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATc                                                                                                    >  2:1040654/1‑150 (MQ=255)
                                                       cAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGt                                                                                                  >  1:1242538/1‑150 (MQ=255)
                                                       cAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGACAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGt                                                                                                  >  1:349370/1‑150 (MQ=255)
                                                          cGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCACCGGAAACGCTGGCGGTTTTTTGTATCACGAAGAAGTTAAGATGAATATTTAGACATCGTtaa                                                                                               >  2:1442813/1‑149 (MQ=255)
                                                            aTAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATAc                                                                                             <  1:681167/150‑1 (MQ=255)
                                                             tAAACCGGTGTGGATGTATCAGGCAGAAAAGTCAACTTGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACAGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATAcc                                                                                            <  2:272286/150‑1 (MQ=255)
                                                              aaaCCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCa                                                                                           >  2:416298/1‑150 (MQ=255)
                                                                     tgtgGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGaaa                                                                                    <  1:1160553/150‑1 (MQ=255)
                                                                     ggtgGATGGATCATTCCGAAAAGGAAACAGGTTATTTTAATTGTTTGCTTATTTTTGAATTCTTTGAAAACCGCCAGCGTAAACGGTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGaaa                                                                                    <  1:1296094/149‑1 (MQ=255)
                                                                      gtgGATGGATCAGGCCGAACAGGCAACTGGTGATGTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCTGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAAc                                                                                   <  1:263075/150‑1 (MQ=255)
                                                                      gtgGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAAc                                                                                   <  2:1206121/150‑1 (MQ=255)
                                                                        ggatggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACat                                                                                 >  1:335768/1‑150 (MQ=255)
                                                                           tggatCAGGCCGAAAAGGCAACTGGGGATTTTTATTGATTGCTGAGTTTTTATTACTTTGATAACCGCAAGCGGAAACGCGGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACatata                                                                              <  2:1555119/150‑1 (MQ=255)
                                                                             gatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATAtt                                                                            >  2:565538/1‑150 (MQ=255)
                                                                                 aGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATg                                                                        <  1:1385224/150‑1 (MQ=255)
                                                                                 aGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCGGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATg                                                                        <  2:851085/150‑1 (MQ=255)
                                                                                    ccGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGaa‑‑a                                                                   >  2:1379282/1‑150 (MQ=255)
                                                                                     cGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑at                                                                  <  1:135041/150‑1 (MQ=255)
                                                                                     cGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTAAGAGGAATATTTAGACATAGTTATACCAATCGAAACATATAATAATGaa‑‑aa                                                                  >  2:1668781/1‑149 (MQ=255)
                                                                                      gAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATTTTCATGAA‑‑ata                                                                 >  2:1143609/1‑150 (MQ=255)
                                                                                       aaaaGGCAACTGGTGATTTTTATTGATTGCTGTGTTTTGAATTCATTGAAAACCGCCATCGGAATCTCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCTTGAA‑‑atat                                                                <  1:1507282/150‑1 (MQ=255)
                                                                                        aaaGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑atata                                                               <  2:173547/150‑1 (MQ=255)
                                                                                           ggCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑atatatat                                                            >  2:237896/1‑150 (MQ=255)
                                                                                           ggCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑atatatat                                                            >  1:1653650/1‑150 (MQ=255)
                                                                                                 tGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAatat                                                      <  2:932371/150‑1 (MQ=255)
                                                                                                  ggTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATAtt                                                     >  1:915149/1‑150 (MQ=255)
                                                                                                  ggTGATTTTTATTGATTGATGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGAGGTTTTTTGTATAAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATAtt                                                     >  1:31819/1‑150 (MQ=255)
                                                                                                    tGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATAtttt                                                   <  2:1693028/150‑1 (MQ=255)
                                                                                                     taTTTTTATTGATTGCTGAGTTTTGAATACATTAAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTc                                                  <  1:1272102/149‑1 (MQ=255)
                                                                                                     gATTTTTCTTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTc                                                  <  1:406134/150‑1 (MQ=255)
                                                                                                     gATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCATGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTc                                                  <  1:486066/150‑1 (MQ=255)
                                                                                                     gATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTc                                                  >  1:662453/1‑150 (MQ=255)
                                                                                                       tttttATTGATTGCTGCGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACCTCGTTATACCACTCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTCCt                                                <  2:265746/150‑1 (MQ=255)
                                                                                                       tttttATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCTGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATTTATATATAtttcatg                                                <  1:99913/150‑5 (MQ=255)
                                                                                                       tttttATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTCCt                                                <  1:416298/150‑1 (MQ=255)
                                                                                                        ttttATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTCCTa                                               >  1:311367/1‑150 (MQ=255)
                                                                                                               gattgCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTATACATCGTTATAACAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTACTAATTGTTc                                        >  2:361905/1‑150 (MQ=255)
                                                                                                                attgCTGAGTTTTTAATACATTGAAAACCGCCAGCTGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTCCTAATTGTTCt                                       <  1:1528188/150‑1 (MQ=255)
                                                                                                                attgCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTCCTAATTGTTCt                                       <  2:473331/150‑1 (MQ=255)
                                                                                                                        gTTTTGAATACATTGAAAACCGCCAGCGGAAACGCGGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑AAATATATAAATATTTTCCTAATTGTTCTTATATGaa                               >  1:1429746/1‑149 (MQ=255)
                                                                                                                              aaTACATTGAAAACCGCCAGCGGAAACGCGGGCGGGTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGatat                         >  1:163792/1‑150 (MQ=255)
                                                                                                                               aTACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATAtc                        <  1:1040654/150‑1 (MQ=255)
                                                                                                                                  cATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCa                     >  2:511709/1‑150 (MQ=255)
                                                                                                                                  cATTGAAAACCCCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTAAGAGGAATAATTAGACATCGTTATACAAAAAGAAACATAAATTCATGAA‑‑AAATATATAAATATTTACCTAATTGTTCTTATATTACAGATATATCa                     >  2:1705615/1‑150 (MQ=255)
                                                                                                                                     tGAAAACCGCCAGCTGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACtt                  <  1:717508/150‑1 (MQ=255)
                                                                                                                                      gAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTa                 >  2:1133508/1‑150 (MQ=255)
                                                                                                                                       aaaaCCGCCAGCGTTTCCGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTaa                >  2:550544/1‑150 (MQ=255)
                                                                                                                                        aaaCCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGGTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAg               >  1:1721066/1‑150 (MQ=255)
                                                                                                                                         aaCCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGTAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAgg              <  1:1112111/150‑1 (MQ=255)
                                                                                                                                          aCCGCCAGCGTTTCCGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACCTCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGc             >  2:1097266/1‑150 (MQ=255)
                                                                                                                                            cGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCtt           >  2:813864/1‑150 (MQ=255)
                                                                                                                                             gCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCttt          >  2:540094/1‑150 (MQ=255)
                                                                                                                                             gCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCttt          >  2:615402/1‑150 (MQ=255)
                                                                                                                                              ccAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTc         <  1:75167/150‑1 (MQ=255)
                                                                                                                                              ccAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTc         <  2:1523373/150‑1 (MQ=255)
                                                                                                                                               cAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGAAATATAACTTAAGGCttt          >  1:1051095/1‑150 (MQ=255)
                                                                                                                                               cAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCt        <  2:775778/150‑1 (MQ=255)
                                                                                                                                               cAGCGGAAACGCGGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTTTACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCt        <  2:231379/150‑1 (MQ=255)
                                                                                                                                               cAGAGGAAACGCTGGCGGTTTTTTGTATCAGTAAGATGTTCAGAGGAATATTTAGACATCTTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCt        <  1:497528/150‑1 (MQ=255)
                                                                                                                                                aGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCtt       <  2:1516706/150‑1 (MQ=255)
                                                                                                                                                    gAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATaa   >  1:1360732/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
TCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAA  >  NC_000913/729228‑729526

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: