Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 729,378 A→T 6.4% 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 6.4% 398.4 / 13.8 140intergenic (+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 (69/62);  new base T (6/3);  total (75/65)
Fisher's exact test for biased strand distribution p-value = 5.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.41e-01

Only 100 of 143 total aligned reads displayed.

CTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAAT  >  NC_000913/729229‑729527
                                                                                                                                                     |                                                                                                                                                       
  ggAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAAc                                                                                                                                                       >  1:885122/1‑150 (MQ=255)
       gACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGATTACCAGCGATTAACCGGTGTGGATGGTTCAGGCCGAAAAGGCAACTGGTGATTTTTTTTGATTGGTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTgg                                                                                                                                                  <  1:348704/150‑1 (MQ=255)
         ccGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCg                                                                                                                                                <  1:477171/150‑1 (MQ=255)
          ctcgcCCGCTTTTACAGCTACCGGCATCGCGTCTGGGCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCgg                                                                                                                                               <  1:1393774/148‑1 (MQ=255)
           gcgcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGt                                                                                                                                              >  2:1892803/1‑150 (MQ=255)
            cgcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGtt                                                                                                                                             <  1:1929410/150‑1 (MQ=255)
            cgcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGtt                                                                                                                                             >  1:1771418/1‑150 (MQ=255)
              cccGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGTTTCCGCTGGCGGtttt                                                                                                                                           <  2:737179/150‑1 (MQ=255)
              cccGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGTTTCCGCTGGCGGtttt                                                                                                                                           <  2:730075/150‑1 (MQ=255)
                 gCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTg                                                                                                                                      <  2:1197408/150‑1 (MQ=255)
                        cAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAgg                                                                                                                               >  1:537855/1‑150 (MQ=255)
                        cAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAgg                                                                                                                               >  1:1341248/1‑150 (MQ=255)
                        cAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGGTTT‑‑TTGTATCAgg                                                                                                                               >  1:1799437/1‑150 (MQ=255)
                        cAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATGTTTATTGATTGCGGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAgg                                                                                                                               >  1:317363/1‑150 (MQ=255)
                        cAGCTACCTGCATCGCGGCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGAGCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAgg                                                                                                                               <  1:1312852/150‑1 (MQ=255)
                         aGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGa                                                                                                                              <  2:1064723/150‑1 (MQ=255)
                          gCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGaa                                                                                                                             >  1:963511/1‑150 (MQ=255)
                          gCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGaa                                                                                                                             >  2:1104568/1‑150 (MQ=255)
                             aCCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGAt                                                                                                                          >  1:1187835/1‑150 (MQ=255)
                              ccTGCATCGCGTCTGTTCGAAAACCAGCGATAAACCGGTGTGTATGTATCAGGCGTAAAAGTCAACTGGGGATTTTAATTGATTGATGATTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATg                                                                                                                         <  2:403320/150‑1 (MQ=255)
                              ccTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATg                                                                                                                         <  2:1417953/150‑1 (MQ=255)
                               cagCATCGCGGCTTGTCGAAACCCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGt                                                                                                                        <  2:82685/148‑1 (MQ=255)
                               cTGCATCGCGTCTGTTCGAAATCCAGCGATAAACCTGTGTGGATGGATCAGGCCGATAAGGCAACTGGTGATTTTTATTGATTGCTGTGTTTTGAATACATTGAAAACCGCCAGCGGCAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGt                                                                                                                        <  1:312717/150‑1 (MQ=255)
                                tGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGtt                                                                                                                       >  1:987616/1‑150 (MQ=255)
                                tGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGtt                                                                                                                       >  2:1097682/1‑150 (MQ=255)
                                 gCATCGCGTCTGGGCGAAAACCAGCGATAAACCGGTGTGGATGTATCAGGCCGAAAAGGCAACTTGGTATTTTTATTGATTGCTGTGTTTTTAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTc                                                                                                                      <  1:1482489/150‑1 (MQ=255)
                                   aTCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATGTTTATAGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCag                                                                                                                    <  1:505966/150‑1 (MQ=255)
                                      gcgTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAgg                                                                                                                 >  2:254012/1‑150 (MQ=255)
                                       cgTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCCCCAGCGGAAACGCTGGCGGTTTT‑‑TTTTATCAGGAAGATGTTCAGAGGa                                                                                                                >  1:1270770/1‑150 (MQ=255)
                                           tGGTCGAAAACCAGCGTTAAACGGGTGTTGATGGTTCAGGCCGAAAAGGCAACTGGTGATTTTAATTGATTGTTGAGTTTTGAATACATTGAAAACCTCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAatat                                                                                                            <  2:406808/150‑1 (MQ=255)
                                            ggTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATAtt                                                                                                           >  1:670557/1‑150 (MQ=255)
                                             gTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATAttt                                                                                                          >  2:68159/1‑150 (MQ=255)
                                              tCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTa                                                                                                         >  2:1298057/1‑150 (MQ=255)
                                                gAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGa                                                                                                       >  2:1375656/1‑150 (MQ=255)
                                                    aCCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATc                                                                                                   >  1:1839449/1‑150 (MQ=255)
                                                      cAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCGGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGt                                                                                                 <  2:1025035/150‑1 (MQ=255)
                                                      cAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGt                                                                                                 >  2:1442135/1‑150 (MQ=255)
                                                       aGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGtt                                                                                                >  2:1667807/1‑150 (MQ=255)
                                                       aGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGtt                                                                                                >  2:430424/1‑150 (MQ=255)
                                                       aGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGtt                                                                                                >  2:2014482/1‑150 (MQ=255)
                                                        gCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTta                                                                                               <  2:1708339/150‑1 (MQ=255)
                                                            tAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGAATTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATAcc                                                                                           >  2:840898/1‑150 (MQ=255)
                                                               aCCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAAt                                                                                        <  1:1813281/150‑1 (MQ=255)
                                                                 cGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTAATTGATTGCTGAGTTTTGAAAACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCg                                                                                      <  2:1016812/150‑1 (MQ=255)
                                                                  ggTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGa                                                                                     >  1:1295583/1‑150 (MQ=255)
                                                                    tgtgGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGaaa                                                                                   <  2:560390/150‑1 (MQ=255)
                                                                       ggatggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACat                                                                                >  2:468387/1‑150 (MQ=255)
                                                                       ggatggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACat                                                                                >  1:1014297/1‑150 (MQ=255)
                                                                       ggatggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACat                                                                                >  1:1359006/1‑150 (MQ=255)
                                                                        gatggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACata                                                                               >  1:950118/1‑150 (MQ=255)
                                                                         atggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACatat                                                                              >  2:1177423/1‑150 (MQ=255)
                                                                             atCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTc                                                                          <  2:120573/150‑1 (MQ=255)
                                                                              tCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCa                                                                         <  1:1358694/150‑1 (MQ=255)
                                                                              tCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCa                                                                         <  1:1024439/150‑1 (MQ=255)
                                                                                aGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATAAATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATg                                                                       <  1:377517/150‑1 (MQ=255)
                                                                                  gCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGaa                                                                     >  1:210887/1‑150 (MQ=255)
                                                                                   ccGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTGTT‑‑TTGTATCAGGAAGATGTTCAGAGGATTATTTAGACATCGTTATACCAATCGAAACATATATTCATGaaa                                                                    >  2:337234/1‑150 (MQ=255)
                                                                                           gCAACTGTTGATTTTTATTTATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAAtatatata                                                            <  2:881776/150‑1 (MQ=255)
                                                                                             aaCTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATaaa                                                          <  2:1113263/150‑1 (MQ=255)
                                                                                              aCTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAat                                                         <  2:1560159/150‑1 (MQ=255)
                                                                                               cTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAata                                                        >  2:787251/1‑150 (MQ=255)
                                                                                                    taTTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTc                                                   <  1:1603519/149‑1 (MQ=255)
                                                                                                    gATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCCGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGTATATTTAGACATCGTTTTACCAATCGAAACATATATTCATGAAATATATATAAATAtttta                                                   >  2:1465012/1‑149 (MQ=255)
                                                                                                     tttttttttGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTcc                                                  <  2:1270770/149‑1 (MQ=255)
                                                                                                     aTTTTTATTTATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTcc                                                  <  2:1515390/150‑1 (MQ=255)
                                                                                                     aTTTTTATTGATTGCTGAGTTTTTAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTcc                                                  <  1:1004700/150‑1 (MQ=255)
                                                                                                      tttttATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATTAAATATATATAAATATTTTCCt                                                 >  2:1184927/1‑150 (MQ=255)
                                                                                                           attgattgCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTg                                            <  2:1385324/150‑1 (MQ=255)
                                                                                                              tattgCTGAGTTTTGAATACATTGAAAACCGCCATCGGAAACGCTGGCTGTTTT‑‑TTGTATCAGGAAGATGTTCAGATGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTc                                         <  2:1234431/149‑1 (MQ=255)
                                                                                                              tattgATGAGTTTTGAATACATTGAAAACCGGCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCGTAATTGTTc                                         <  2:1433401/149‑1 (MQ=255)
                                                                                                              gattgCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTc                                         <  2:1580579/150‑1 (MQ=255)
                                                                                                               atagCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCt                                        <  2:1279411/147‑1 (MQ=255)
                                                                                                                  gCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCCGGAAGATGTTCAGAGGAATATTTAGACATCGTGATACGAATCGATACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTAt                                     >  2:326558/1‑150 (MQ=255)
                                                                                                                   cTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATc                                    >  1:1108448/1‑150 (MQ=255)
                                                                                                                    tGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCt                                   <  2:885122/150‑1 (MQ=255)
                                                                                                                       gTTTTGAATACATTGAAAACCGCCAGAGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGAc                                <  2:944963/150‑1 (MQ=255)
                                                                                                                         ttttAATACATTTAAAACCGCCAGCGGAAACTCTTGCGGTTTT‑‑TTGTATCATGAAGATGTTCATAGTAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAg                              <  1:1523493/150‑1 (MQ=255)
                                                                                                                         tttGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAg                              >  2:266883/1‑150 (MQ=255)
                                                                                                                          ttGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATAAATAAATATTTTCCTAATTGTTCTTATCTGACAGa                             <  2:1960143/150‑1 (MQ=255)
                                                                                                                            taaTACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGata                           <  1:251043/149‑1 (MQ=255)
                                                                                                                             aaTACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGatat                          >  2:1788976/1‑150 (MQ=255)
                                                                                                                              aTACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATAtc                         >  2:551710/1‑150 (MQ=255)
                                                                                                                               tACATTGAAAACCGCCAGCGTAAACGCTGGCGGTTTT‑‑TTTTATCAGGAAGATTTTCAGAGGAATATTTATAAAGCGGTATAGCAATCGAAACATATATTCATGAAATTTATATAAATATTTTCCTAATTGTTATTTTCTGACAGatatat                        >  2:284336/1‑148 (MQ=255)
                                                                                                                                aCATTTAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATAtctc                       <  1:1814666/150‑1 (MQ=255)
                                                                                                                                 cATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTGTTATCTGACAGATAtctaa                      >  1:155684/1‑148 (MQ=255)
                                                                                                                                  aTTGAAAACCGCCCGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCGCAc                     >  1:1393294/1‑150 (MQ=255)
                                                                                                                                  aTTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCAc                     <  2:13128/150‑1 (MQ=255)
                                                                                                                                    aaaaaaCCGCCAGCGTTTCCGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATATGACAGATATCTCACtt                   >  1:2010682/3‑150 (MQ=255)
                                                                                                                                     gAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGATACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGGTTTCTCTCttt                  >  2:762231/1‑149 (MQ=255)
                                                                                                                                     gAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTa                  >  1:312996/1‑150 (MQ=255)
                                                                                                                                       aaaCCGCCAGCGTTTCCGCTGGCGGTTTTCAATGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATATCAGTTa                  >  1:1943192/1‑150 (MQ=38)
                                                                                                                                       aaaCCGCCAGCGTTTCCGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTAACTTAAg                >  1:1327288/1‑150 (MQ=255)
                                                                                                                                       aaaCCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAg                >  1:1636714/1‑150 (MQ=255)
                                                                                                                                       aaaCCGCAAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAg                <  1:835218/150‑1 (MQ=255)
                                                                                                                                         aCCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAgga              >  1:371937/1‑149 (MQ=255)
                                                                                                                                          ccGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCt             >  2:452762/1‑150 (MQ=255)
                                                                                                                                          ccGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCt             >  1:608285/1‑150 (MQ=255)
                                                                                                                                              cAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCt         >  2:1053913/1‑150 (MQ=255)
                                                                                                                                               aGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCtt        <  2:963511/150‑1 (MQ=255)
                                                                                                                                                    aaaCGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATaaa   <  2:1899013/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
CTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAAT  >  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: