Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 824,573 T→A 6.9% intergenic (‑76/‑57) ybhP ← / → ybhQ endonuclease/exonuclease/phosphatase domain‑containing protein YbhP/putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913824,5730TA 6.9% 318.8 / 8.8 100intergenic (‑76/‑57)ybhP/ybhQendonuclease/exonuclease/phosphatase domain‑containing protein YbhP/putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (45/48);  new base A (3/4);  total (48/52)
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

Only 100 of 103 total aligned reads displayed.

TAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCTGGTGGTGGGC  >  NC_000913/824425‑824722
                                                                                                                                                    |                                                                                                                                                     
tAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAAttt                                                                                                                                                      >  2:55257/1‑150 (MQ=255)
  aaCGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTc                                                                                                                                                    <  1:911833/150‑1 (MQ=255)
    cgcgTTAAAGTCTTTGTTAATATTGATGGTGAGCACCTTGAAAGAATATTGTTGTGTTTGATAGGGCATAATTTTCCTGTCTTTGCCTATTATCTCTTTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTAATGAATTTTCCg                                                                                                                                                  <  2:416171/150‑1 (MQ=255)
    cgcgGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCg                                                                                                                                                  <  1:1222756/150‑1 (MQ=255)
       ggTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTaa                                                                                                                                               <  2:308492/150‑1 (MQ=255)
        gTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAAAATTCCGTaag                                                                                                                                              <  1:1444198/150‑2 (MQ=255)
         tAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAg                                                                                                                                             <  2:1426749/150‑1 (MQ=255)
         tAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAg                                                                                                                                             <  1:25381/150‑1 (MQ=255)
         tAAAGCCTTTGTGAATATTGATCGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGACTTTGCCTCTTATCTCATTGAAATAGTGAAGTCGGCGTCAAAAAAAGGTGCGGTCTTACGGAAATATCCGTAAAg                                                                                                                                             >  2:1628285/1‑150 (MQ=255)
          aaaGCCTTTTTGAATATTGATGGTGTTAAACTTGAAAGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTATTTTCCTCTTATATCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCTGTCTTACGGAAAATTCCGTaagac                                                                                                                                            <  2:278206/150‑4 (MQ=255)
          aaaGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGt                                                                                                                                            <  1:510448/150‑1 (MQ=255)
          aaaGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAAAATTCCGTaagac                                                                                                                                            <  1:1617069/150‑4 (MQ=255)
           aaGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAAAATTCCGTaagacc                                                                                                                                           <  1:330463/150‑5 (MQ=255)
             gCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCg                                                                                                                                         >  2:619393/1‑150 (MQ=255)
               cTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGt                                                                                                                                       <  2:79750/150‑1 (MQ=255)
                 ttGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTAc                                                                                                                                     <  2:50714/150‑1 (MQ=255)
                 ttGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTAc                                                                                                                                     <  2:1614389/150‑1 (MQ=255)
                 ttGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGCTCTTACGGAATTTTCCGTAAAGTTCGGTAc                                                                                                                                     <  2:392856/150‑1 (MQ=255)
                     gAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTTCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTg                                                                                                                                 <  2:341437/150‑1 (MQ=255)
                        tatTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagt                                                                                                                              <  1:1065990/150‑1 (MQ=255)
                         atTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCTCAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagta                                                                                                                             >  2:1423110/1‑150 (MQ=255)
                         atTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagta                                                                                                                             >  2:1292828/1‑150 (MQ=255)
                         atTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagta                                                                                                                             >  2:288417/1‑150 (MQ=255)
                            gATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagtaagt                                                                                                                          >  1:1498617/1‑150 (MQ=255)
                            gATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagtaagt                                                                                                                          <  1:231778/150‑1 (MQ=255)
                                      aCCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAt                                                                                                                <  2:1261910/150‑1 (MQ=255)
                                          tGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAttctt                                                                                                            >  1:191505/1‑150 (MQ=255)
                                             aCGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAg                                                                                                         >  2:330481/1‑150 (MQ=255)
                                             aCGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAg                                                                                                         <  1:1249393/150‑1 (MQ=255)
                                              cGAAAATTGTTGTGTTTGATCGGGCATAATATTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAgg                                                                                                        >  1:521593/1‑150 (MQ=255)
                                               gAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAAATTTCCGTAAAGTTCGGTACTCTGAGTAAGAAGAGATAAATTCTTCAGga                                                                                                       >  2:1061088/1‑150 (MQ=255)
                                                 aaaTTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGgaga                                                                                                     <  2:692355/150‑1 (MQ=255)
                                                     tgttgtGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAg                                                                                                 <  2:62438/150‑1 (MQ=255)
                                                       ttgtTTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGcc                                                                                               <  1:1139240/150‑1 (MQ=255)
                                                         gtgtTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCAt                                                                                             <  2:572322/150‑1 (MQ=255)
                                                          tgtTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATg                                                                                            >  1:741951/1‑150 (MQ=255)
                                                           gtTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATATCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGa                                                                                           >  1:1384995/1‑150 (MQ=255)
                                                              tcaTCGGCCATAATTTTCCTGTCTTTGCCACTTATCTCAATGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGCTCTTACGAAATTTTACGTAAAGTTCGGTACTCAGAGTAAGTAGAGATAAATTCTTCAAGAGAGAACCCATGAAGt                                                                                        >  2:8298/3‑150 (MQ=255)
                                                              tGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGt                                                                                        >  1:1139659/1‑150 (MQ=255)
                                                                aTCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTATTCAGGAGAGAAGCCATGAAGTgg                                                                                      >  1:654417/1‑150 (MQ=255)
                                                                   gggCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGc                                                                                     <  2:1354533/148‑1 (MQ=255)
                                                                   gggCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGc                                                                                     >  1:1354533/1‑148 (MQ=255)
                                                                    ggCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTag                                                                                                                        <  2:1048938/112‑1 (MQ=255)
                                                                    ggCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTag                                                                                                                        >  1:1048938/1‑112 (MQ=255)
                                                                        tAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGCACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACg                                                                              <  1:372286/150‑1 (MQ=255)
                                                                         aaTTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACgt                                                                             >  2:182378/1‑150 (MQ=255)
                                                                          aTTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACgtg                                                                            <  1:25795/150‑1 (MQ=255)
                                                                                cTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCgtg                                                                      >  1:562085/1‑150 (MQ=255)
                                                                                cTGTCTTTGCCTCTTATCTCAATGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCgtg                                                                      >  2:572523/1‑150 (MQ=255)
                                                                                 tGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCgtgt                                                                     <  1:1066739/150‑1 (MQ=255)
                                                                                    cTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGc                                                                  >  1:1422315/1‑150 (MQ=255)
                                                                                    cTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGc                                                                  >  2:1250497/1‑150 (MQ=255)
                                                                                       agCCACTTATCTCAATGAAAAACTTAATTAAGCTTCACAAAAACGTCCTGTCTTAAAGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAAc                                                               >  2:606287/2‑150 (MQ=255)
                                                                                        gCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCGCAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACg                                                              >  1:1438074/1‑150 (MQ=255)
                                                                                        gCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACg                                                              >  1:1237879/1‑150 (MQ=255)
                                                                                          ctctTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCgtgt                                                                     >  2:388142/1‑141 (MQ=255)
                                                                                          ctctTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCgtgt                                                                     <  1:388142/141‑1 (MQ=255)
                                                                                          ctctTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACggg                                                            <  1:330481/150‑1 (MQ=255)
                                                                                           tctTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGt                                                           >  1:751627/1‑150 (MQ=255)
                                                                                           tctTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCAGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGt                                                           >  2:1604978/1‑150 (MQ=255)
                                                                                                tctcATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAg                                                      >  1:598887/1‑150 (MQ=255)
                                                                                                 ctcATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGt                                                     <  1:1570451/150‑1 (MQ=255)
                                                                                                   cATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTg                                                   >  2:1579260/1‑150 (MQ=255)
                                                                                                     ttGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCt                                                 >  1:164327/1‑150 (MQ=255)
                                                                                                      tGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTg                                                >  2:381241/1‑150 (MQ=255)
                                                                                                       gAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTgg                                               >  2:825213/1‑150 (MQ=255)
                                                                                                         aaTAGTGTAGTCGGCGTCTCAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGGGTCGCAACGGGTCTAAGTTGCTGGCa                                             <  1:1423110/150‑1 (MQ=255)
                                                                                                          aTAGTGTAGTCGGTGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAg                                            <  1:273618/150‑1 (MQ=255)
                                                                                                           tagtgtagtgGGGGTCACATAAAGGTGCGGTCTTACGGAATTTTCCGTAATGTTCGGCACTCTGAGTAAGTTGAGATTAATTCTTCAGTAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGa                                           <  1:1628285/150‑1 (MQ=255)
                                                                                                             gtgtAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGAtt                                         <  1:55257/150‑1 (MQ=255)
                                                                                                              tgtAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTa                                        <  1:572523/150‑1 (MQ=255)
                                                                                                               gtAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTAt                                       >  2:1320973/1‑150 (MQ=255)
                                                                                                                tAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATg                                      <  1:177855/150‑1 (MQ=255)
                                                                                                                 aGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGt                                     <  1:1292828/150‑1 (MQ=255)
                                                                                                                       cGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAg                                                      <  2:744839/127‑1 (MQ=255)
                                                                                                                       cGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAg                                                      >  1:744839/1‑127 (MQ=255)
                                                                                                                          cacaAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTAc                            >  2:686423/1‑150 (MQ=255)
                                                                                                                           tcaAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACt                           >  2:186838/2‑150 (MQ=255)
                                                                                                                           acaAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACt                           >  1:1307844/1‑150 (MQ=255)
                                                                                                                            caAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTg                          >  1:1176181/1‑150 (MQ=255)
                                                                                                                            caAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTg                          >  1:1022957/1‑150 (MQ=255)
                                                                                                                             aaaaaGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTgg                         <  2:1139659/150‑1 (MQ=255)
                                                                                                                             aaaaaGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTgg                         <  2:455337/150‑1 (MQ=255)
                                                                                                                             aaaaaGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTgg                         <  1:182378/150‑1 (MQ=255)
                                                                                                                              aaaaGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGt                        <  1:1597774/150‑1 (MQ=255)
                                                                                                                              aaaaGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGt                        <  1:64306/150‑1 (MQ=255)
                                                                                                                                    ccgaactttACGGAAAATTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGgc                  >  2:857205/8‑150 (MQ=255)
                                                                                                                                         ctttACGGAAAATTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCgctgg             >  2:18759/3‑149 (MQ=255)
                                                                                                                                          cTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCgctgct            >  1:1372460/1‑150 (MQ=255)
                                                                                                                                          cTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCgctgct            <  2:1498617/150‑1 (MQ=255)
                                                                                                                                          cTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCgctgct            <  2:226114/150‑1 (MQ=255)
                                                                                                                                           ttACGGAAAATTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGctgctg           >  1:1262101/1‑150 (MQ=255)
                                                                                                                                            tACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCtgg          <  1:288417/150‑1 (MQ=255)
                                                                                                                                             aCGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCtggt         >  2:1014886/1‑150 (MQ=255)
                                                                                                                                              cGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAATCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCtggtg        >  2:1259649/1‑150 (MQ=255)
                                                                                                                                              cGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCtggtg        >  2:739648/1‑150 (MQ=255)
                                                                                                                                                 aaTTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCtggtggtg     <  1:606287/150‑1 (MQ=255)
                                                                                                                                                  aTTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGGCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCtggtggtgg    <  2:713608/150‑1 (MQ=255)
                                                                                                                                                   ttttCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCTGGTGGTggg   >  1:1592074/1‑150 (MQ=255)
                                                                                                                                                    tttCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCTGGTGGTGGGc  <  1:1579260/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TAAACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCTGGTGGTGGGC  >  NC_000913/824425‑824722

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: