Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 729,379 A→C 5.2% 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 5.2% 404.9 / 7.2 140intergenic (+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 (69/62);  new base C (5/3);  total (75/65)
Fisher's exact test for biased strand distribution p-value = 7.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 142 total aligned reads displayed.

GGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAAT  >  NC_000913/729231‑729527
                                                                                                                                                    |                                                                                                                                                      
ggAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAAc                                                                                                                                                       >  2:1361092/1‑150 (MQ=255)
     gACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGATTACCAGCGATTAACCGGTGTGGATGGTTCAGGCCGAAAAGGCAACTGGTGATTTTTTTTGATTGGTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTgg                                                                                                                                                  <  1:348704/150‑1 (MQ=255)
        ctcgcCCGCTTTTACAGCTACCGGCATCGCGTCTGGGCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCgg                                                                                                                                               <  1:1393774/148‑1 (MQ=255)
        cgcgCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAACCCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGTTTCCGCTGGCgg                                                                                                                                               <  1:1311567/150‑1 (MQ=255)
         gcgcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGt                                                                                                                                              >  2:1892803/1‑150 (MQ=255)
          cgcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGtt                                                                                                                                             <  1:1929410/150‑1 (MQ=255)
           gcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGGGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGttt                                                                                                                                            <  2:1548486/150‑1 (MQ=255)
            cccGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGTTTCCGCTGGCGGtttt                                                                                                                                           <  2:737179/150‑1 (MQ=255)
            cccGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGTTTCCGCTGGCGGtttt                                                                                                                                           <  2:730075/150‑1 (MQ=255)
               gCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTg                                                                                                                                      <  2:1197408/150‑1 (MQ=255)
               gCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTg                                                                                                                                      <  1:919148/150‑1 (MQ=255)
                cTTTTACAGCTACCTGCCTCGCGTCTGGTCGAATACCAGCGATTAGCCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGt                                                                                                                                     >  2:302179/1‑150 (MQ=255)
                      cAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAgg                                                                                                                               >  1:537855/1‑150 (MQ=255)
                      cAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAgg                                                                                                                               >  1:1341248/1‑150 (MQ=255)
                      cAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAgg                                                                                                                               >  1:242905/1‑150 (MQ=255)
                      cAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGGTTT‑‑TTGTATCAgg                                                                                                                               >  1:1799437/1‑150 (MQ=255)
                      cAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATGTTTATTGATTGCGGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAgg                                                                                                                               >  1:317363/1‑150 (MQ=255)
                        gCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCATCGGAAACGCTTGCGGTTTT‑‑TTGTATCAGGaa                                                                                                                             >  2:377517/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                                                                                                                       >  2:1097682/1‑150 (MQ=255)
                              tGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGtt                                                                                                                       >  1:987616/1‑150 (MQ=255)
                               gCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTc                                                                                                                      <  2:1229134/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)
                                      gTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGaa                                                                                                               >  2:1047004/1‑150 (MQ=255)
                                         tGGTCGAAAACCAGCGTTAAACGGGTGTTGATGGTTCAGGCCGAAAAGGCAACTGGTGATTTTAATTGATTGTTGAGTTTTGAATACATTGAAAACCTCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAatat                                                                                                            <  2:406808/150‑1 (MQ=255)
                                           gTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATAttt                                                                                                          >  1:1400438/1‑150 (MQ=255)
                                           gTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATAttt                                                                                                          >  2:68159/1‑150 (MQ=255)
                                           gTCGAAAACCAGCGATAAACCGGTGTGAATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATAttt                                                                                                          >  2:1142119/1‑150 (MQ=255)
                                            tCGAAAACCAGCGATAAACCGGTGTGGATGGATCCGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCTCCAGCGGAAACGCTTGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTa                                                                                                         >  2:806403/1‑150 (MQ=255)
                                            tCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTa                                                                                                         >  2:1298057/1‑150 (MQ=255)
                                             cGAAAACCAGCTATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAg                                                                                                        <  2:525624/150‑1 (MQ=255)
                                              gAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGa                                                                                                       >  2:1012488/1‑150 (MQ=255)
                                              gAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGa                                                                                                       >  2:1375656/1‑150 (MQ=255)
                                                  aCCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATc                                                                                                   >  1:1839449/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)
                                                               cGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTAATTGATTGCTGAGTTTTGAAAACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCg                                                                                      <  2:1016812/150‑1 (MQ=255)
                                                                ggTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGa                                                                                     >  1:1295583/1‑150 (MQ=255)
                                                                  tgtgGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGaaa                                                                                   <  2:560390/150‑1 (MQ=255)
                                                                    tggatggaTCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACa                                                                                 >  1:350953/1‑150 (MQ=255)
                                                                     ggatggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACat                                                                                >  1:1359006/1‑150 (MQ=255)
                                                                     ggatggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACat                                                                                >  1:1014297/1‑150 (MQ=255)
                                                                      gatggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACata                                                                               >  1:950118/1‑150 (MQ=255)
                                                                       atggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACatat                                                                              >  2:1177423/1‑150 (MQ=255)
                                                                            tCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCa                                                                         <  1:1358694/150‑1 (MQ=255)
                                                                            tCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCa                                                                         <  1:1024439/150‑1 (MQ=255)
                                                                             cAGGCAGAAAAGGCAACTGGTGATTTTTATTGATTGATGAGTAGTGAATACATTGAAACCCCCCAGCGGAAACGCGGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCAt                                                                        <  2:1070812/150‑1 (MQ=255)
                                                                              aGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCCCCCGCGGAAACGCTGGCGGGTTT‑‑TTGTATCAGGAAGAATTTCAGAGGAATATTTATACATCGATATACCAATCGAAACatttaataatt                                                                       >  1:1964532/1‑144 (MQ=255)
                                                                              aGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATAAATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATg                                                                       <  1:377517/150‑1 (MQ=255)
                                                                                gCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGaa                                                                     >  1:210887/1‑150 (MQ=255)
                                                                                         gCAACTGTTGATTTTTATTTATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAAtatatata                                                            <  2:881776/150‑1 (MQ=255)
                                                                                           aaCTGGTTATTTTTATTTATTGCTGATTTTTGAATACATTTAAAACCGCCAGCGGAAACTCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATaaa                                                          <  1:1801689/150‑1 (MQ=255)
                                                                                           aaCTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATaaa                                                          <  2:1113263/150‑1 (MQ=255)
                                                                                            aCTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAat                                                         <  2:1560159/150‑1 (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)
                                                                                                   aTTTTTATTGATTGCTGAGTTTTTAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTcc                                                  <  1:1004700/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)
                                                                                                             attgCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCt                                        <  2:1839449/150‑1 (MQ=255)
                                                                                                                gCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCCGGAAGATGTTCAGAGGAATATTTAGACATCGTGATACGAATCGATACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTAt                                     >  2:326558/1‑150 (MQ=255)
                                                                                                                 cTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATc                                    >  1:1108448/1‑150 (MQ=255)
                                                                                                                 cTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATc                                    <  2:1400438/150‑1 (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)
                                                                                                                            aTACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATAtc                         >  2:551710/1‑150 (MQ=255)
                                                                                                                              aCATTTAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATAtctc                       <  1:1814666/150‑1 (MQ=255)
                                                                                                                               cATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTGTTATCTGACAGATAtctaa                      >  1:155684/1‑148 (MQ=255)
                                                                                                                               cATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCa                      <  1:787251/150‑1 (MQ=255)
                                                                                                                                aTTGAAAACCGCCCGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCGCAc                     >  1:1393294/1‑150 (MQ=255)
                                                                                                                                 caaaaaaCCGCCAGCGTTTCCGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACt                    >  2:1948338/4‑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‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATTAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAg                <  2:1290634/150‑1 (MQ=255)
                                                                                                                                     aaaCCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAg                >  1:1636714/1‑150 (MQ=255)
                                                                                                                                     aaaCCGCAAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAg                <  1:835218/150‑1 (MQ=255)
                                                                                                                                        ccGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCt             >  1:608285/1‑150 (MQ=255)
                                                                                                                                        ccGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCt             >  2:452762/1‑150 (MQ=255)
                                                                                                                                        ccGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCt             >  2:1443542/1‑150 (MQ=255)
                                                                                                                                            cAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCt         >  2:1053913/1‑150 (MQ=255)
                                                                                                                                             aGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCtt        >  1:1707310/1‑150 (MQ=255)
                                                                                                                                             aGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCtt        <  2:963511/150‑1 (MQ=255)
                                                                                                                                              gCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTta       <  2:1187835/150‑1 (MQ=255)
                                                                                                                                               cGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTtat      <  1:460405/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
GGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTT‑‑TTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAAT  >  NC_000913/729231‑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: