Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 729,379 A→C 9.8% intergenic (+39/‑204) 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,3790AC 9.8% 264.9 / 14.4 107intergenic (+39/‑204)ybfA/rhsCDUF2517 family protein/Rhs protein with putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (47/48);  new base C (4/7);  total (51/56)
Fisher's exact test for biased strand distribution p-value = 5.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 113 total aligned reads displayed.

TGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATC  >  NC_000913/729230‑729528
                                                                                                                                                     |                                                                                                                                                       
tGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGaaa                                                                                                                                                         >  1:157847/1‑150 (MQ=255)
tGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGaaa                                                                                                                                                         <  1:1488806/150‑1 (MQ=255)
tGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGaaa                                                                                                                                                         <  2:1004500/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)
          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)
                                                                                                                             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)
                                                                                                                                           gCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCttt            >  2:4757/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)
                                                                                                                                                    aaCGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAAt   >  2:631738/1‑150 (MQ=255)
                                                                                                                                                     aCGCTGGAGTTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATc  <  1:193430/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
TGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATC  >  NC_000913/729230‑729528

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: