Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 729,378 A→T 8.1% intergenic (+38/‑205) 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,3780AT 8.1% 292.6 / 11.8 107intergenic (+38/‑205)ybfA/rhsCDUF2517 family protein/Rhs protein with putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (49/49);  new base T (2/7);  total (51/56)
Fisher's exact test for biased strand distribution p-value = 1.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.34e-01

Only 100 of 113 total aligned reads displayed.

CTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAAT  >  NC_000913/729229‑729527
                                                                                                                                                     |                                                                                                                                                       
cTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCACCAGAGGaa                                                                                                                                                         >  1:1516339/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                                                                                                                                           <  2:1508451/150‑1 (MQ=255)
              cccGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGTTTCCGCTGGCGGtttt                                                                                                                                           <  1:1260901/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)
                                                                                                                                            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)
                                                                                                                                                     |                                                                                                                                                       
CTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAA‑‑ATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAAT  >  NC_000913/729229‑729527

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: