Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 729,379 A→C 5.3% 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.3% 278.1 / 4.8 82intergenic (+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 (33/44);  new base C (2/3);  total (35/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACAAAAGGGGGAGAAAAGTATGGAACTCTACAGAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCA  >  NC_000913/729114‑729550
                                                                                                                                                                                                                                                                         |                                                                                                                                                                           
tACAAAAGGGGGAGAAAAGTATGGAACTCTACAGAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCgg                                                                                                                                                                    >  2:166371/1‑275 (MQ=255)
    aaaGGGGGAGAAAAGTATGGAACTCTACAGAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGtttt                                                                                                                                                                <  2:241567/275‑1 (MQ=255)
                   tATGGAACTCTACAGAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGAt                                                                                                                                                 >  4:160786/1‑275 (MQ=255)
                         acgctACAGAGAAGTGCCTGCATGTACTATCTTTGTACGCCGTACGTATGCGGTTGCAGCGGGCGTTCGTGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCGACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACAGTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCa                                                                                                                                            <  4:271957/271‑1 (MQ=255)
                                      ttgccTGCATGGCTTATCTTTTTACGCGGTACTGATGCGGTTGCCGCGGGCCTTCTGGCGCTGCCTTTCATACTTTTGTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGCTCAGAGGAATATTTAGa                                                                                                                              <  4:195245/272‑1 (MQ=255)
                                       atCCTGCATGGCTTATCTTTTTACGCCGTACTTATGCGGTTGCATCGGGCGTTCTTGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGa                                                                                                                              <  2:330682/274‑1 (MQ=255)
                                       atCCTGCATGGCTTATCTTTTTACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGCTAAACCGGTGGGGATGGATCCGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAGACGCCAGCGGAAACGGTGGCGGGTTTTTGTATCGGGAAGATGTTCAGAGGAATATTTAGa                                                                                                                              >  2:261916/1‑274 (MQ=255)
                                       atCCTGCATGGCTTATCTTTTTACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGa                                                                                                                              <  4:166371/274‑1 (MQ=255)
                                         ccTGCTTGGATTAGCTTGTTACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTTGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACCGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATAAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACa                                                                                                                            <  2:200886/274‑1 (MQ=255)
                                         ccTGCATGGCTTATCTTTTTACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCCACTGGTGATTTTTATTGATTGCTGAGGTTTGAATACGTTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATACGGAAGCTGTTCAGAGGAATATTTAgcact                                                                                                                           >  4:200886/1‑271 (MQ=255)
                                              aTGGCTTATCTTTTTACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGCGGATGGAACAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTTAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTACCAGGAAGATGCTCAAAAGAATATTTAGACCTCGTta                                                                                                                      >  2:195245/1‑275 (MQ=255)
                                                     aTCTTTTTTCGCCGGACTTTTTCGGTTGCATCGGTCGTTCTTGCGCTGCCTTTCATGCTTTCCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAAt                                                                                                               <  4:261916/275‑1 (MQ=255)
                                                        tttttACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGa                                                                                                            <  2:160786/275‑1 (MQ=255)
                                                                    aCTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATAtt                                                                                                                                  >  4:284892/1‑241 (MQ=255)
                                                                        aTGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGa                                                                                                                              <  2:284892/241‑1 (MQ=255)
                                                                             gTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACTTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCCGAGGAATATTTACAAATCGTTATTCCAATCGAAACATATATTCATGACata                                                                                         >  2:87543/1‑273 (MQ=255)
                                                                                                               ttttCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCa                                                                                                                                                        >  4:370310/1‑176 (MQ=255)
                                                                                                                   cTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGaa                                                                                                                                                    <  2:370311/176‑1 (MQ=255)
                                                                                                                      gAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAAc                                                                                                                                                                            <  7:272685/149‑1 (MQ=255)
                                                                                                                        aaGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGc                                                                                                                                                                          >  7:394595/1‑149 (MQ=255)
                                                                                                                           aCCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTgg                                                                                                                                                                       <  8:394424/149‑1 (MQ=255)
                                                                                                                           aCCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCa                                                                                                                                            >  3:64774/1‑176 (MQ=255)
                                                                                                                            ccGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAATCCGCCAGCGGAACCGCTGGc                                                                                                                                                                      <  8:1293440/149‑1 (MQ=255)
                                                                                                                              gcgcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCgg                                                                                                                                                                    <  7:665681/149‑1 (MQ=255)
                                                                                                                              gcgcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCgg                                                                                                                                                                    >  6:217781/1‑149 (MQ=255)
                                                                                                                               cgcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAgg                                                                                                                                        <  1:64774/176‑1 (MQ=255)
                                                                                                                                    gCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGtttttt                                                                                                                                                              >  8:739980/1‑149 (MQ=255)
                                                                                                                                              cTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGaa                                                                                                                                                    >  6:716384/1‑149 (MQ=255)
                                                                                                                                              cTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGaa                                                                                                                                                    <  5:716384/149‑1 (MQ=255)
                                                                                                                                                aCCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGACTACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTAGTGAAATATATAAAAATATTTTCCTAATTGTTATTATTTGACAGATATACAACTTAAGGCTTTCTGATaaacatgtg                    >  4:9172/1‑269 (MQ=255)
                                                                                                                                                 ccTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCACCTGGCGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGAt                                                                                                                                                 <  7:1376479/149‑1 (MQ=255)
                                                                                                                                                 ccTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGAt                                                                                                                                                 >  7:291477/1‑149 (MQ=255)
                                                                                                                                                      aTCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCa                                                                                                                                            <  8:394595/149‑1 (MQ=255)
                                                                                                                                                      aTCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCa                                                                                                                                            >  5:719412/1‑149 (MQ=255)
                                                                                                                                                        cgcgTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCaga                                                                                                                                          <  5:217781/149‑1 (MQ=255)
                                                                                                                                                         gcgTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCagag                                                                                                                                         <  7:695876/149‑1 (MQ=255)
                                                                                                                                                                   gAAAACCAGCGATAAACCGGTGTGGATGGATCATGCCGAAAAGGCAACTGGTGATTGTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTcagca  <  2:84886/274‑1 (MQ=255)
                                                                                                                                                                   gAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGtt                                                                                                                                                                  <  7:430770/114‑1 (MQ=255)
                                                                                                                                                                   gAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGtt                                                                                                                                                                  >  8:430770/1‑114 (MQ=255)
                                                                                                                                                                          aGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGTTTCCGCTGGCGGt                                                                                                                                                                   >  7:1062768/1‑106 (MQ=39)
                                                                                                                                                                          aGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGTTTCCGCTGGCGGt                                                                                                                                                                   <  8:1062768/106‑1 (MQ=39)
                                                                                                                                                                          aGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGt                                                                                                                        <  5:1082624/149‑1 (MQ=255)
                                                                                                                                                                            cGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTta                                                                                                                      <  7:739980/149‑1 (MQ=255)
                                                                                                                                                                               tAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATAc                                                                                                                   <  5:724550/149‑1 (MQ=255)
                                                                                                                                                                                   ccGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTta                                                                                                                      >  5:32621/1‑142 (MQ=255)
                                                                                                                                                                                   ccGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTta                                                                                                                      <  6:32621/142‑1 (MQ=255)
                                                                                                                                                                                   ccGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAAt                                                                                                               >  4:286046/1‑149 (MQ=255)
                                                                                                                                                                                      gtgtGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATAttt                                                                                                                                 >  2:363149/1‑128 (MQ=255)
                                                                                                                                                                                       tgtgGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGaa                                                                                                           <  2:286046/149‑1 (MQ=255)
                                                                                                                                                                                          ggatggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGAc                                                                                                                             <  4:363149/128‑1 (MQ=255)
                                                                                                                                                                                          ggatggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGaaa                                                                                                          <  5:1464903/147‑1 (MQ=255)
                                                                                                                                                                                          ggatggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGaaa                                                                                                          >  6:1464903/1‑147 (MQ=255)
                                                                                                                                                                                           gatggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACat                                                                                                       >  5:971584/1‑149 (MQ=255)
                                                                                                                                                                                            atggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACata                                                                                                      <  6:198288/149‑1 (MQ=255)
                                                                                                                                                                                             tggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACatat                                                                                                     <  8:1150890/149‑1 (MQ=255)
                                                                                                                                                                                               gatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTa                                                                                                                                <  8:40075/120‑1 (MQ=255)
                                                                                                                                                                                               gatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTa                                                                                                                                >  7:40075/1‑120 (MQ=255)
                                                                                                                                                                                                atCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATAGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATg                                                                                                                                                >  5:967573/1‑103 (MQ=255)
                                                                                                                                                                                                atCAGGCCGAAAAGGCAACTGGTGATTTTTAGTGATTGCTGAGTTTTGAATACATTGAAAACCTCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATg                                                                                                                                                <  6:967573/103‑1 (MQ=255)
                                                                                                                                                                                                 tCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTCGTATCAGGAAGACGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTc                                                                                                 <  8:1070864/149‑1 (MQ=255)
                                                                                                                                                                                                   aGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGaaa                                                                                                          >  6:909872/1‑138 (MQ=255)
                                                                                                                                                                                                   aGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGaaa                                                                                                          <  5:909872/138‑1 (MQ=255)
                                                                                                                                                                                                   aGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCAt                                                                                               <  6:971584/149‑1 (MQ=255)
                                                                                                                                                                                                           aaGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAatata                                                                                       <  7:239195/149‑1 (MQ=255)
                                                                                                                                                                                                            aGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAatatat                                                                                      <  8:291477/149‑1 (MQ=255)
                                                                                                                                                                                                              gCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCaa                                                                                                                >  8:993631/1‑121 (MQ=255)
                                                                                                                                                                                                              gCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCaa                                                                                                                <  7:993631/121‑1 (MQ=255)
                                                                                                                                                                                                                                   ttgCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCt                                                               >  6:1236252/1‑149 (MQ=255)
                                                                                                                                                                                                                                    tgCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATAtt                                                                                                                                  >  2:261784/1‑81 (MQ=255)
                                                                                                                                                                                                                                       tGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATc                                                           >  5:79624/1‑149 (MQ=255)
                                                                                                                                                                                                                                        gAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGa                                                                                                                              <  4:261784/81‑1 (MQ=255)
                                                                                                                                                                                                                                          gTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATc                                                                                                              <  7:754366/95‑1 (MQ=255)
                                                                                                                                                                                                                                          gTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATc                                                                                                              >  8:754366/1‑95 (MQ=255)
                                                                                                                                                                                                                                               gAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGat                                                   >  5:413798/1‑149 (MQ=255)
                                                                                                                                                                                                                                                 aTACATTGAAACCCGCCAGCGGTAACGCTGGCGGTTTTTTGTATCAGGAACATGTTCTGTGGAATATTTAGATATCGTTATACCTATCGAAACTTATATTCATGCAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGatat                                                 <  8:308323/149‑1 (MQ=255)
                                                                                                                                                                                                                                                       aaaaaaCCGCCAGCGTTTCCGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGt                                                                  >  6:573791/3‑126 (MQ=39)
                                                                                                                                                                                                                                                       aaaaaaCCGCCAGCGTTTCCGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGt                                                                  <  5:573791/124‑1 (MQ=39)
                                                                                                                                                                                                                                                            aCCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAgg                                      <  5:1236252/149‑1 (MQ=255)
                                                                                                                                                                                                                                                               gCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGttt              >  4:151130/1‑170 (MQ=255)
                                                                                                                                                                                                                                                                 cAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTc                                 <  6:413798/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                   gCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGcc          <  2:151130/170‑1 (MQ=255)
                                                                                                                                                                                                                                                                      gAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTtata                            >  7:233957/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         |                                                                                                                                                                           
TACAAAAGGGGGAGAAAAGTATGGAACTCTACAGAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCA  >  NC_000913/729114‑729550

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: