Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,984,542 T→A 7.5% Q489L (CAG→CTG)  aslA ← putative Ser‑type periplasmic non‑aryl sulfatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,984,5420TA 7.5% 948.7 / 11.9 298Q489L (CAG→CTG) aslAputative Ser‑type periplasmic non‑aryl sulfatase
Reads supporting (aligned to +/- strand):  ref base T (76/192);  new base A (23/2);  total (104/194)
Fisher's exact test for biased strand distribution p-value = 3.51e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 302 total aligned reads displayed.

GGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTC  >  NC_000913/3984393‑3984691
                                                                                                                                                     |                                                                                                                                                     
 gATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTg                                                                                                                                                      <  2:3094280/150‑1 (MQ=255)
  aTCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGa                                                                                                                                                     <  1:3138696/150‑1 (MQ=255)
  aTCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGa                                                                                                                                                     <  1:3772559/150‑1 (MQ=255)
      ccATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGGAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATc                                                                                                                                                 >  2:2705727/1‑150 (MQ=255)
           tccgcgTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGctc                                                                                                                                            <  2:3365980/148‑1 (MQ=255)
            aCGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGctct                                                                                                                                           >  2:2854236/1‑150 (MQ=255)
             cgcgTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTg                                                                                                                                          <  1:750309/150‑1 (MQ=255)
             cgcgTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTg                                                                                                                                          <  2:4097148/150‑1 (MQ=255)
              gcgTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTgg                                                                                                                                         <  1:4184916/150‑1 (MQ=255)
                 tGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTTCGGATCGGTGTAGAGGTTAAACACCCACCATCCCGCCGTTTTCATTAAAGTGCCGGTGAAAACACCCAGAAATCCGCGCTGGGta                                                                                                                                      >  2:487925/1‑150 (MQ=255)
                 tGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGta                                                                                                                                      <  2:1687419/150‑1 (MQ=255)
                   cATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTTTAGGTTAAACACCCACGATACCGCCGGTTTCAATTATTTTCCGGTGTAAACAACCAGATATACGCTCAGGGgata                                                                                                                                    >  1:4098563/1‑150 (MQ=255)
                    aTTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATaa                                                                                                                                   <  1:2463148/150‑1 (MQ=255)
                          gTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTaa                                                                                                                             <  1:3026365/150‑1 (MQ=255)
                          gTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTaa                                                                                                                             <  1:897632/150‑1 (MQ=255)
                             tGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTACTTCGGATAGGTGTATTGGTTAAAAACAGACGATACCGCCGTTTTGATTTATTTGGCGGTGGATACACCCTGATATAAGCTCTGGGTTTAATCGTAATGt                                                                                                                          >  1:2051908/1‑150 (MQ=255)
                             tGTAGCGGTAAACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTACTGCGGATCGGTGTAGAGGTTAAAAACCCACCAAACCCCCGGTTTGATTAAAGTGCCGGTGGAAACAACCCCATAACCGCGCGGGGGGTAACCGGAAGGt                                                                                                                          >  1:2092369/1‑150 (MQ=255)
                              gTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGtt                                                                                                                         <  2:2783309/150‑1 (MQ=255)
                              gTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGtt                                                                                                                         <  2:2564825/150‑1 (MQ=255)
                              gTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGtt                                                                                                                         <  1:346974/150‑1 (MQ=255)
                               tAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTg                                                                                                                        <  1:494777/150‑1 (MQ=255)
                               tAGCGGTACACCCATCGGAATATGGCACACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTg                                                                                                                        <  1:2752212/150‑1 (MQ=255)
                                 gCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCt                                                                                                                      <  2:2453582/150‑1 (MQ=255)
                                   ggTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGa                                                                                                                    <  1:1539343/150‑1 (MQ=255)
                                      acacCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATc                                                                                                                 <  2:1363500/150‑1 (MQ=255)
                                       cacCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCa                                                                                                                <  1:564257/150‑1 (MQ=255)
                                       cacCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCa                                                                                                                <  1:391917/150‑1 (MQ=255)
                                            aTCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACg                                                                                                           <  1:4030797/150‑1 (MQ=255)
                                             tCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGt                                                                                                          <  1:2517675/150‑1 (MQ=255)
                                                 aaTATGGCGCACGCCGATGGAGTCGCTTTACTGCGGATCGGTGTAGAGGTTAAAAACCCACGATCCCCCCCTTTGCAGTACTGTGCCGGGGAATCCACCCAGATATACCCGCTGGGTATAAGCGGAAGGTTGGTGAATCAGGACGTGATa                                                                                                      >  2:425327/1‑150 (MQ=255)
                                                 aaTACGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATa                                                                                                      <  1:1291345/150‑1 (MQ=255)
                                                 aaTACGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGAGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATa                                                                                                      <  1:2595378/150‑1 (MQ=255)
                                                     tGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTg                                                                                                  <  2:4108113/150‑1 (MQ=255)
                                                      tgcgcACGCCGATGGAGTCGCTTTCCTGCGGATCTGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGa                                                                                                 <  1:3587634/149‑1 (MQ=255)
                                                       gcgcACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTTAGGTTGCTGAATCAGGACGTGATACTTGaa                                                                                                <  1:731750/150‑1 (MQ=255)
                                                       gcgcACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGaa                                                                                                <  1:2117095/150‑1 (MQ=255)
                                                             gCCGATGGAGTCGCTTTCCTGCGGATAGGTGTAGAGGTTAAACACCGACCATCCCGCCCTTTTGATTACAGTGCCCGTGAATCCAACCAGTTATACGCACTGGGTATAAGCCTAAGGTTGCTTAAACAGGACGTGATAAATGAAATCaaa                                                                                          >  2:792495/1‑148 (MQ=255)
                                                               cGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCa                                                                                        <  1:1019347/150‑1 (MQ=255)
                                                                gATGGAGTCGCTTTCCTGGGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCAt                                                                                       <  2:1854699/150‑1 (MQ=255)
                                                                    gAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCCACGATCCCGCCGTTTGCATTACATTGCCGGTGAATCCACCCTGATATCCGCACAGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATAATTGAAATCATCCATACGc                                                                                   >  2:3533600/1‑150 (MQ=255)
                                                                         gCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCAcagc                                                                              >  2:1345539/1‑150 (MQ=255)
                                                                           tttCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCAcagcag                                                                            <  2:3141900/150‑1 (MQ=255)
                                                                            ttCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACagcagc                                                                           <  1:227455/150‑1 (MQ=255)
                                                                            ttCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACagcagc                                                                           <  1:2405325/150‑1 (MQ=255)
                                                                                tGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGt                                                                       <  1:3336892/150‑1 (MQ=255)
                                                                                tGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGt                                                                       >  2:1069630/1‑150 (MQ=255)
                                                                                 gCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGtt                                                                      <  1:1252941/150‑1 (MQ=255)
                                                                                          tgtAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGa                                                             <  2:4235819/150‑1 (MQ=255)
                                                                                          tgtAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGa                                                             <  1:3218627/150‑1 (MQ=255)
                                                                                            tAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAgg                                                           >  2:3849577/1‑150 (MQ=255)
                                                                                            tAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAgg                                                           >  2:172337/1‑150 (MQ=255)
                                                                                             agagGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGa                                                          >  2:2236911/1‑150 (MQ=255)
                                                                                             agagGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGa                                                          >  2:821017/1‑150 (MQ=255)
                                                                                             agagGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGGGAATCCACCCTGATATACGCTCTGGGTATAAGCGTAAGGTTGCTGAAACAGGACGTGATACTTGAAATCATCCATAAGCACAGCAGCGAGTTTACCGTTGAGGa                                                          >  1:2723115/1‑150 (MQ=255)
                                                                                              gagGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGaa                                                         <  1:1445467/150‑1 (MQ=255)
                                                                                               agGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAag                                                        <  1:3923093/150‑1 (MQ=255)
                                                                                               agGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAag                                                        >  2:1209751/1‑150 (MQ=255)
                                                                                                  ttAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCAGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAagtag                                                     >  2:1890398/1‑150 (MQ=255)
                                                                                                      acacCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCt                                                 >  2:3940875/1‑150 (MQ=255)
                                                                                                         ccGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCgg                                              <  2:4165764/150‑1 (MQ=255)
                                                                                                         ccGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCgg                                              >  2:3056340/1‑150 (MQ=255)
                                                                                                         ccGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCgg                                              <  1:532936/150‑1 (MQ=255)
                                                                                                         ccGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCgg                                              >  2:3592571/1‑150 (MQ=255)
                                                                                                                tCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACg                                       <  2:3459874/150‑1 (MQ=255)
                                                                                                                tCCCGCCGTTTGCATGACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACg                                       <  1:4095836/150‑1 (MQ=255)
                                                                                                                   cgccgTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGtt                                    >  2:3970636/1‑150 (MQ=255)
                                                                                                                   cgccgTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGtt                                    >  2:102772/1‑150 (MQ=255)
                                                                                                                   cgccgTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGtt                                    >  2:1091172/1‑150 (MQ=255)
                                                                                                                   cgccgTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGtt                                    >  2:3264213/1‑150 (MQ=255)
                                                                                                                    gccgTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTa                                   <  1:2147385/150‑1 (MQ=255)
                                                                                                                    gccgTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTa                                   <  2:3678277/150‑1 (MQ=255)
                                                                                                                    gccgTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTa                                   <  2:3907998/150‑1 (MQ=255)
                                                                                                                      cgTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGa                                 <  1:998455/150‑1 (MQ=255)
                                                                                                                          tGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGa                             <  2:3184176/150‑1 (MQ=255)
                                                                                                                          tGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGa                             <  1:833856/150‑1 (MQ=255)
                                                                                                                             aTTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGAccc                          >  2:1521580/1‑149 (MQ=255)
                                                                                                                               tACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCAtt                        <  1:3808220/150‑1 (MQ=255)
                                                                                                                                aCTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCAttt                       <  1:1077325/150‑1 (MQ=255)
                                                                                                                                aCTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCAttt                       <  2:2742616/150‑1 (MQ=255)
                                                                                                                                aCTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTGTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCAttt                       <  1:425327/150‑1 (MQ=255)
                                                                                                                                  tgtgCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGt                     >  1:1183511/1‑150 (MQ=255)
                                                                                                                                     gCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTcc                  <  2:2951139/150‑1 (MQ=255)
                                                                                                                                     gCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTcc                  <  2:2275183/150‑1 (MQ=255)
                                                                                                                                     gCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACAGTTAGACTGACCATTTGTTcc                  <  2:558080/150‑1 (MQ=255)
                                                                                                                                      ccGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTccc                 <  2:954733/150‑1 (MQ=255)
                                                                                                                                      ccGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTccc                 <  1:1973075/150‑1 (MQ=255)
                                                                                                                                      ccGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTccc                 <  1:2870707/150‑1 (MQ=255)
                                                                                                                                      ccGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTccc                 <  1:3091229/150‑1 (MQ=255)
                                                                                                                                       cGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCa                <  1:2023869/150‑1 (MQ=255)
                                                                                                                                            aaTCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGgaag           >  1:1705731/1‑150 (MQ=255)
                                                                                                                                            aaTCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGCTCCCAGgaag           >  1:68464/1‑150 (MQ=255)
                                                                                                                                              tCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGgaagaa         <  2:2060583/150‑1 (MQ=255)
                                                                                                                                              tCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGgaagaa         <  2:46525/150‑1 (MQ=255)
                                                                                                                                                 aCCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGa      >  1:384949/1‑150 (MQ=255)
                                                                                                                                                 aCCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGa      <  2:3602942/150‑1 (MQ=255)
                                                                                                                                                 aCCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGa      >  2:2621824/1‑150 (MQ=255)
                                                                                                                                                   ccTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACTCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATg    >  1:3625824/1‑150 (MQ=255)
                                                                                                                                                   ccTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATg    >  1:1835988/1‑150 (MQ=255)
                                                                                                                                                    cTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGt   <  1:2704389/150‑1 (MQ=255)
                                                                                                                                                     tGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTc  <  2:2092495/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTC  >  NC_000913/3984393‑3984691

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