Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 824,572 T→A 5.2% intergenic (‑75/‑58) ybhP ← / → ybhQ endo/exonuclease/phosphatase family protein/inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913824,5720TA 5.2% 305.5 / 4.0 106intergenic (‑75/‑58)ybhP/ybhQendo/exonuclease/phosphatase family protein/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (52/48);  new base A (3/3);  total (55/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

Only 100 of 112 total aligned reads displayed.

GTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCTGGTGGTGGG  >  NC_000913/824423‑824721
                                                                                                                                                     |                                                                                                                                                     
gTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAAt                                                                                                                                                       <  1:1068733/150‑1 (MQ=255)
 ttAAAAGAGGTAAAGCATTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTAGTGTTTGATCGGGCATAATTATCCTGTCTTTGCCTCTTATCTCAGTGAAATAGTGAAGTCTGCGTCACAAAAATGTGCGGTCTTACGGAAtt                                                                                                                                                      <  1:266298/150‑1 (MQ=255)
   aaaCGAGGTAAACCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATCGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAAtttt                                                                                                                                                    <  2:1024481/150‑1 (MQ=255)
    aaCGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTc                                                                                                                                                   <  2:1113815/150‑1 (MQ=255)
    aaCGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTc                                                                                                                                                   >  2:635519/1‑150 (MQ=255)
    aaCGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTc                                                                                                                                                   >  2:1493614/1‑150 (MQ=255)
       gcgGTAAAGCCTTTGTGAATATTGATGGTGATCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCAGCGTCACAAAAAGGTGCGGTCTTACGGAAAATTCCGt                                                                                                                                                <  2:1454897/150‑1 (MQ=255)
       gcgGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGt                                                                                                                                                <  1:185425/150‑1 (MQ=255)
       gcgGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCAGCGTCACAAAAAGGTGCGGTCTTACGGAAAATTCCGt                                                                                                                                                <  2:1448167/150‑1 (MQ=255)
        cgGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGTCGTCACAAAAAGGTGCGGTCTTACGGAAAATTCCGTa                                                                                                                                               <  2:920735/150‑1 (MQ=255)
              aGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTc                                                                                                                                         >  1:1325349/1‑150 (MQ=255)
              aGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCAGTAAAGTTc                                                                                                                                         >  2:341075/1‑150 (MQ=255)
               gCCTTTGTGAATTTTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCg                                                                                                                                        <  2:17668/150‑1 (MQ=255)
               gCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCGTGTCTTTGCCTCTTATCTCATTGATATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTTCGGAATTTTCCGTAAAGTTCg                                                                                                                                        >  2:910543/1‑150 (MQ=255)
               gCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCg                                                                                                                                        <  1:229551/150‑1 (MQ=255)
                ccTTTGTGAATATTGGTGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCgg                                                                                                                                       <  2:1067725/150‑1 (MQ=255)
                ccTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGGCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCgg                                                                                                                                       <  1:1216566/150‑1 (MQ=255)
                ccTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCATCTCTTTGCCTCTTATCTCATTGAAATAGAGAAGTCGTCGTCACAAAAAGGTGCGGTCTTAACGAAATTTCCGTAAAGTTCgg                                                                                                                                       >  2:1696125/1‑150 (MQ=255)
                 cTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGt                                                                                                                                      <  2:853597/150‑1 (MQ=255)
                     gtgAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTActc                                                                                                                                  <  2:572237/150‑1 (MQ=255)
                     gtgAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGACTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTAAGTAAAGTTCTGTActc                                                                                                                                  >  2:765354/1‑150 (MQ=255)
                      tgAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTActct                                                                                                                                 >  1:1559910/1‑150 (MQ=255)
                       gAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTg                                                                                                                                <  1:1007014/150‑1 (MQ=255)
                       gAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCAGTAAAGTTCGGTACTCTg                                                                                                                                >  2:1331903/1‑150 (MQ=255)
                        aaTATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGa                                                                                                                               >  2:1573531/1‑150 (MQ=255)
                          tatTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagt                                                                                                                             <  2:538842/150‑1 (MQ=255)
                          taatGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagt                                                                                                                             <  1:6369/147‑1 (MQ=255)
                              gATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAACGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagtaagt                                                                                                                         <  2:1367711/150‑1 (MQ=255)
                                     aGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATa                                                                                                                  >  1:1584202/1‑150 (MQ=255)
                                        aCCTTGAACTAAAATTCTTGTGTTTGATCGGGCATAATTATCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAATAAGGTGCTGTCTTACGGAATTTTCCGTAATGTTCGGTACTCTGAGTAAGTTGAGATAAAt                                                                                                               <  2:103518/150‑1 (MQ=255)
                                         ccTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTATAGATAAAtt                                                                                                              >  2:124271/1‑150 (MQ=255)
                                         ccTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAtt                                                                                                              >  2:854513/1‑150 (MQ=255)
                                         ccTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAtt                                                                                                              >  2:1608734/1‑150 (MQ=255)
                                           ttGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAttct                                                                                                            <  2:297684/150‑1 (MQ=255)
                                                cGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAgg                                                                                                       >  2:298773/1‑150 (MQ=255)
                                                  aaaaTTGAGGGGTTGCATCGGGCATAATTATCCAGTCTTTGCCACTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGgag                                                                                                     <  2:318566/150‑1 (MQ=255)
                                                        gttgtGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGc                                                                                               >  1:189884/1‑150 (MQ=255)
                                                             gtTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGa                                                                                          >  1:1428167/1‑150 (MQ=255)
                                                                tGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGt                                                                                       <  1:891122/150‑1 (MQ=255)
                                                                 gATCGGTAATAATTTTCCTGTCTTTGCATCTTATCTCATGGAAATAGTGTTGTCTTCGTCACAAAAAGGTGCGGTCTTACGTAATTTTCCGTAAAGTTAGGTACTCTGAGTAATGTGAGATATATTCTTCAGTAGAGAAGACATTAAGTg                                                                                      >  2:305975/1‑150 (MQ=255)
                                                                 gATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTg                                                                                      <  2:1491579/150‑1 (MQ=255)
                                                                   tCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGTAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGc                                                                                    <  2:1369456/150‑1 (MQ=255)
                                                                    cGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGGCGGCGTCACAAAACGGTGCGGGCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGCGAAGCCATGAAGTGGca                                                                                   <  2:343044/150‑1 (MQ=255)
                                                                         aTAATTTTCATTTCTTTGCCTCTTATCTCTTTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCGTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCaacaac                                                                              <  1:908984/150‑1 (MQ=255)
                                                                         aTAAATTTCCTGTCTTTGCCTCTTAACTCATTGAAATAGTCTAGGCGGCGGCACAAACAGGGGCGGTCTTACGGAATTTTCCGTAAAGTGCGGTACTCTGAGTAAGTAGAGATAAAATCTTCAGGAGAGAAGCCATGAAGTGGCaacaac                                                                              <  1:915069/150‑1 (MQ=255)
                                                                          tAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGAGGTCTTACGGAATTTTACGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCTGGAGAGAAGCCATGAAGTGGCaactact                                                                             >  1:262128/1‑149 (MQ=255)
                                                                             ttttCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACgtgt                                                                          >  1:682088/1‑150 (MQ=255)
                                                                               ttCCTGTCTTTGCCTCTTATCTCATTGAAATAGTTTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTc                                                                        <  2:1400796/150‑1 (MQ=255)
                                                                                 ccTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCgt                                                                      <  1:894616/150‑1 (MQ=255)
                                                                                  cTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCgtg                                                                     <  2:682088/150‑1 (MQ=255)
                                                                                  cTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCgtg                                                                     >  1:1732139/1‑150 (MQ=255)
                                                                                    gTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTc                                                                   <  1:29955/150‑1 (MQ=255)
                                                                                        ttGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATATTCCGTAAAGATCGGTACTCTGAGTAAGAAGAGATAAATTCTTCAGGAGAGAATCCATGAAGAGGCAACAACGTGTTCGTGTCGCaa                                                               >  2:530319/1‑150 (MQ=255)
                                                                                         tGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTATGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAAc                                                              >  2:1527512/1‑150 (MQ=255)
                                                                                         tGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCATCAACGTCTTCGTGTCGCAAc                                                              >  2:895015/1‑150 (MQ=255)
                                                                                         tGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAAc                                                              <  1:854513/150‑1 (MQ=255)
                                                                                          gCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGATTTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACg                                                             >  1:155499/1‑150 (MQ=255)
                                                                                          gCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACg                                                             <  1:975810/150‑1 (MQ=255)
                                                                                           ccTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACgg                                                            >  2:460308/1‑150 (MQ=255)
                                                                                            ctctTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACggg                                                           >  1:1564375/1‑150 (MQ=255)
                                                                                              ctTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTGAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTGGTGTCGCAACggggc                                                         >  2:1163200/1‑148 (MQ=255)
                                                                                               ttATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCt                                                        <  2:1413782/150‑1 (MQ=255)
                                                                                               ttATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCt                                                        <  2:378553/150‑1 (MQ=255)
                                                                                               tattatcATTGAAATATTGTAGTCTGCGTAACAAATAGGTGAGTTCATACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCAATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCt                                                        <  2:1293754/145‑1 (MQ=255)
                                                                                                tATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTa                                                       <  2:189884/150‑1 (MQ=255)
                                                                                                   ctcATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGt                                                    >  1:1400413/1‑150 (MQ=255)
                                                                                                   ctcATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGt                                                    >  1:1561024/1‑150 (MQ=255)
                                                                                                     cATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTg                                                  >  1:672045/1‑150 (MQ=255)
                                                                                                         gAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGCTAAATACTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGGCTAAGTTGCTgg                                              >  2:335070/1‑150 (MQ=255)
                                                                                                         gAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTgg                                              >  2:1578127/1‑150 (MQ=255)
                                                                                                           aaTAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTTTTCGTGTCGCAACTGGTCTAAGTTGCTGGCa                                            >  2:121987/1‑150 (MQ=255)
                                                                                                            aTAGTGTAGGCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCTTAAAGTGCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACCTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAg                                           <  1:868408/150‑1 (MQ=255)
                                                                                                             tAGTGTAGTCGGCGTCTCAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGa                                          >  2:466588/1‑150 (MQ=255)
                                                                                                                   aGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCGAAGTTGCTGGCAGATTAtat                                    >  1:48180/1‑148 (MQ=255)
                                                                                                                    gTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGtt                                   <  1:1013524/150‑1 (MQ=255)
                                                                                                                    gTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGtt                                   >  1:331673/1‑150 (MQ=255)
                                                                                                                    gTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGtt                                   <  1:539360/150‑1 (MQ=255)
                                                                                                                      cGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGc                                 >  1:1382809/1‑150 (MQ=255)
                                                                                                                          gTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCTGATTATGTTGCAttt                             >  1:198502/1‑150 (MQ=255)
                                                                                                                              caAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTg                         <  2:867141/150‑1 (MQ=255)
                                                                                                                               aaaaaGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACGCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTgg                        <  2:327591/150‑1 (MQ=255)
                                                                                                                                aaaaGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGTTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGt                       >  2:1125417/1‑150 (MQ=255)
                                                                                                                                aaaaGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGt                       >  2:1371855/1‑150 (MQ=255)
                                                                                                                                   aGGAGCGGTCTAACGGAATTTTCCGGAAAGTTCGTTACTCTGAGTAAGTTGAGATAAATTCTTAAGGAGAGAAGCCATGAAGTGGCAACAACGTGTGCGTGTCGCAACGGGGCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGt                    <  2:244065/150‑1 (MQ=255)
                                                                                                                                    ggTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTg                   >  2:1008181/1‑150 (MQ=255)
                                                                                                                                       gCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGgcg                <  1:362882/150‑1 (MQ=255)
                                                                                                                                           ctttACGGAAAATTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCgctgc            >  1:434344/3‑150 (MQ=255)
                                                                                                                                            cTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCgctgct           >  1:1701258/1‑150 (MQ=255)
                                                                                                                                             ttACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGTTAGTGGCGctgctg          >  1:1291335/1‑150 (MQ=255)
                                                                                                                                               aCGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCACCGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCtggt        <  1:102856/150‑1 (MQ=255)
                                                                                                                                               aCGGAATTTTCAGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCtggt        >  2:1689278/1‑150 (MQ=255)
                                                                                                                                                cGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCtggtg       >  1:1370485/1‑150 (MQ=255)
                                                                                                                                                cGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGCTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCtggtg       >  1:1267376/1‑150 (MQ=255)
                                                                                                                                                cGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTCGTCGTGTCGCAACGGGACTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTAGCGCTGCtggtg       >  1:1284339/1‑150 (MQ=255)
                                                                                                                                                   aaTTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCtggtggtg    <  1:453353/150‑1 (MQ=255)
                                                                                                                                                    aTTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCtggtggtgg   >  1:1676763/1‑150 (MQ=255)
                                                                                                                                                    aTTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCtggtggtgg   <  2:1194410/150‑1 (MQ=255)
                                                                                                                                                     ttttCCGTATAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCTGGTGGTggg  <  2:804339/150‑1 (MQ=255)
                                                                                                                                                     ttttCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCTGGTGGTggg  <  2:1400413/150‑1 (MQ=255)
                                                                                                                                                     ttttCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCTGGTGGTggg  <  2:1561024/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTTAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCTGGTGGTGGG  >  NC_000913/824423‑824721

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: