Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,984,542 T→A 6.3% 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 6.3% 1848.7 / 12.4 592Q489L (CAG→CTG) aslAputative Ser‑type periplasmic non‑aryl sulfatase
Reads supporting (aligned to +/- strand):  ref base T (164/351);  new base A (46/2);  total (239/353)
Fisher's exact test for biased strand distribution p-value = 3.42e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 598 total aligned reads displayed.

GGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTC  >  NC_000913/3984393‑3984691
                                                                                                                                                     |                                                                                                                                                     
   tctcCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGat                                                                                                                                                    <  1:571773/150‑1 (MQ=255)
      ccATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCCACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATc                                                                                                                                                 >  2:8028677/1‑150 (MQ=255)
        atatACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCg                                                                                                                                               <  2:6556726/150‑1 (MQ=255)
         tataCGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTTAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGc                                                                                                                                              <  1:8706282/150‑1 (MQ=255)
             cgcgTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTg                                                                                                                                          <  2:6114152/150‑1 (MQ=255)
             cgcgTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTg                                                                                                                                          <  2:5041166/150‑1 (MQ=255)
             cgcgTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTg                                                                                                                                          <  1:2599387/150‑1 (MQ=255)
              gcgTGCATTTTGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTgg                                                                                                                                         <  1:315066/150‑1 (MQ=255)
              gcgTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTgg                                                                                                                                         <  2:1443327/150‑1 (MQ=255)
              gcgTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTgg                                                                                                                                         <  2:3145671/150‑1 (MQ=255)
              gcgTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCAATAATGTGCCGGTGAATACACCCGGATATACCCTCTgg                                                                                                                                         >  1:5472735/1‑150 (MQ=255)
                gTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGt                                                                                                                                       <  1:4359716/150‑1 (MQ=255)
                    aTTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATaa                                                                                                                                   <  2:8211226/150‑1 (MQ=255)
                     tttCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGGATAAg                                                                                                                                  >  2:5386468/1‑150 (MQ=255)
                      ttCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGAACGGGGGAGAGGTTAAAAACCGACGATCCCGCCCTTTGCCAGAAGGTTGCGGTGAATACACCCAGATATACCCGCGGGGTATaacc                                                                                                                                 >  2:7092203/1‑148 (MQ=255)
                         ggTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGAACCCGCCGTTTGCATTACAGTGCCGGTGAATACACCCAGATATACGCTCTGGGTAAAAGCCTa                                                                                                                              >  2:7007234/1‑150 (MQ=255)
                          gTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTaa                                                                                                                             <  2:5739766/150‑1 (MQ=255)
                          gTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTaa                                                                                                                             <  1:3931593/150‑1 (MQ=255)
                          gTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTACTGCGGATCGGTGTATATGTTAAACACCGACGAACCCGCCCTTTTCATTAAAGTTCCCGTGAAAACAACCAGATATACCCTAGGGGTATAAGCGTaa                                                                                                                             >  2:3260359/1‑150 (MQ=255)
                               tAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTg                                                                                                                        <  1:9046161/150‑1 (MQ=255)
                                 gCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCt                                                                                                                      <  1:1209169/150‑1 (MQ=255)
                                 gCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCt                                                                                                                      <  1:4726415/150‑1 (MQ=255)
                                 gCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCt                                                                                                                      <  1:273117/150‑1 (MQ=255)
                                      acacCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATc                                                                                                                 <  2:8829512/150‑1 (MQ=255)
                                           cATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGTCGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGAc                                                                                                            <  2:6005879/150‑1 (MQ=255)
                                           cATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGAc                                                                                                            <  1:7790835/150‑1 (MQ=255)
                                            aTCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACg                                                                                                           <  1:420838/150‑1 (MQ=255)
                                            aTCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACg                                                                                                           >  2:644722/1‑150 (MQ=255)
                                            aTCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACg                                                                                                           <  1:6006198/150‑1 (MQ=255)
                                             tCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGt                                                                                                          <  1:5819725/150‑1 (MQ=255)
                                             tCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTAATGTGCCGGTGAAAACACCCTGATATACGCTCTGGGTATAAGCGTAAGGTTTCTTAATAAGGACGt                                                                                                          >  2:3627430/1‑150 (MQ=255)
                                             tCGGAATATGGCGCACGCCGATGGAGTCGCGTTCCGGCGGATCGGTGTAGAGGTTAAACACCCACCATCCCCCCGGTTGCAATATTTTGCCGGTGAATCCACCCAGAAATGCGCACGGGGTATAAGCCGAAGGGTGCTGGAACAGGAccg                                                                                                          >  2:912542/1‑148 (MQ=255)
                                                 actcgggCGCACGCCGAGGGAGTCGCGTTCCTGCGGATCGGGGTAGAGGGTAAACACCGACGCGCCCGCCGTGTGCATGACTGTGCCGGTGAATCCACCCTGAGATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGAtc                                                                                                      <  1:400596/145‑2 (MQ=255)
                                                 aaTATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATa                                                                                                      <  2:6650723/150‑1 (MQ=255)
                                                  atatGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATAc                                                                                                     <  1:8904405/150‑1 (MQ=255)
                                                     tGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTg                                                                                                  <  2:7423009/150‑1 (MQ=255)
                                                        cgcACGCCGATGGAGTCGCTTTACTGCGGATCGGTGGAGAGGTTAAAAACCGACGAACCCGCCGTTTGCATTAAATTGCCCGTGAATACACCCGGATATACCCGCGGGGTATAAGCGGAAGGTTGCTGAATAAGGAACTGATaaattaac                                                                                               >  2:1779335/1‑143 (MQ=255)
                                                         gcACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACt                                                                                              <  1:2674785/150‑1 (MQ=255)
                                                            cGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGTTATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCAt                                                                                           <  2:7258471/150‑1 (MQ=255)
                                                            cGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCAt                                                                                           <  2:783769/150‑1 (MQ=255)
                                                            cGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCAt                                                                                           <  2:1363002/150‑1 (MQ=255)
                                                            cGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCAt                                                                                           <  1:1439565/150‑1 (MQ=255)
                                                             gcagATGGAGTCGCTTTCCTGCGGATCTGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATc                                                                                          <  1:8160112/147‑1 (MQ=255)
                                                             gCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATc                                                                                          <  1:7456999/150‑1 (MQ=255)
                                                             gCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCCCCGTTTGCATTACTGTGCCGGTGAATCCACCCGGATATCCGCGCGGGGGATAAGCGGAAGGGTGCTGAATCAGGCCGTGATACTTGAActcctc                                                                                          >  2:522480/1‑150 (MQ=255)
                                                                 aTGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATa                                                                                      <  1:920747/150‑1 (MQ=255)
                                                                          cTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCAcagca                                                                             >  1:6160853/1‑150 (MQ=255)
                                                                            ttCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACagcagc                                                                           <  1:5435744/150‑1 (MQ=255)
                                                                            ttCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACagcagc                                                                           <  2:390173/150‑1 (MQ=255)
                                                                              ccTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGa                                                                         <  1:3749088/150‑1 (MQ=255)
                                                                               cTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGTATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGGTGCTGAATCAGGTCGTGATAATTGAACTCATACATAAGCACAGCAGCGAg                                                                        >  1:4921404/1‑150 (MQ=255)
                                                                               cTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAg                                                                        >  2:2810113/1‑150 (MQ=255)
                                                                                tGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCGGATATACGCGCTGGGTATAAGCGTAAGGTTGCGGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGt                                                                       >  1:1630917/1‑150 (MQ=255)
                                                                                    gATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTAc                                                                   <  1:5232045/150‑1 (MQ=255)
                                                                                    gATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTAc                                                                   <  2:5316304/150‑1 (MQ=255)
                                                                                    gATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTAc                                                                   <  1:6859734/150‑1 (MQ=255)
                                                                                         gtgtAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCAGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATACAAAAGCACAGCAGCGAGTTTACCGTTg                                                              >  1:8482997/1‑150 (MQ=255)
                                                                                         gtgtAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCCGTGAATCCACCCAGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCATCGAGTTTAACGTTg                                                              >  2:896668/1‑150 (MQ=255)
                                                                                          tgtAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGa                                                             <  1:3677031/150‑1 (MQ=255)
                                                                                          tgtAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGa                                                             <  2:8698487/150‑1 (MQ=255)
                                                                                          tgtAGAGGATAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGa                                                             <  2:6916050/150‑1 (MQ=255)
                                                                                           gtAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAg                                                            <  2:3986476/150‑1 (MQ=255)
                                                                                             agagGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGa                                                          >  2:4574415/1‑150 (MQ=255)
                                                                                               agGTTACACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAag                                                        <  1:1825439/150‑1 (MQ=255)
                                                                                               agGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAag                                                        <  2:1158444/150‑1 (MQ=255)
                                                                                                  ttAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAagtag                                                     >  2:3820575/1‑150 (MQ=255)
                                                                                                      acacCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCt                                                 >  2:2923696/1‑150 (MQ=255)
                                                                                                          cgacgaTCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGc                                             <  1:705257/150‑1 (MQ=255)
                                                                                                          cgacgaTCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGc                                             <  1:5153372/150‑1 (MQ=255)
                                                                                                          cgacgaTCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGc                                             <  1:6318070/150‑1 (MQ=255)
                                                                                                          cgacgaTCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTAACGTTGAGGAAGTAGTGCTCGGc                                             >  2:2055465/1‑150 (MQ=255)
                                                                                                           gacgaTCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGcc                                            <  2:5527669/150‑1 (MQ=255)
                                                                                                               aTCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTAc                                        <  2:4339854/150‑1 (MQ=255)
                                                                                                               aTCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTAc                                        <  1:3422455/150‑1 (MQ=255)
                                                                                                               aTCCCGCCGTGTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTAc                                        <  1:6530187/150‑1 (MQ=255)
                                                                                                                   cgccgTTTGCATTACTGTGCCGGTGAATCCACCCAGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGtt                                    >  1:7218785/1‑150 (MQ=255)
                                                                                                                    gccgTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTa                                   <  2:4735193/150‑1 (MQ=255)
                                                                                                                          tGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGa                             >  2:8610712/1‑150 (MQ=255)
                                                                                                                          tGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGa                             >  2:4118332/1‑150 (MQ=255)
                                                                                                                           gCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGAc                            <  1:1041628/150‑1 (MQ=255)
                                                                                                                           gCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGAc                            <  2:2840379/150‑1 (MQ=255)
                                                                                                                                aCTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCAttt                       <  1:2351665/150‑1 (MQ=255)
                                                                                                                                aCTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCAttt                       <  2:4912085/150‑1 (MQ=255)
                                                                                                                                aCTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCAttt                       <  1:6723693/150‑1 (MQ=255)
                                                                                                                                  tgtgCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGt                     >  2:12548/1‑150 (MQ=255)
                                                                                                                                    tgCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTc                   <  1:2473316/150‑1 (MQ=255)
                                                                                                                                       cGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCa                <  1:1980948/150‑1 (MQ=255)
                                                                                                                                       cGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCa                <  1:975012/150‑1 (MQ=255)
                                                                                                                                       cGGTGAAACCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCa                <  1:9016164/150‑1 (MQ=255)
                                                                                                                                           gAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGgaa            <  1:3799613/150‑1 (MQ=255)
                                                                                                                                              tCCACCCTGATGTCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGgaagaa         >  1:7688724/1‑150 (MQ=255)
                                                                                                                                              tCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGgaagaa         >  2:2281135/1‑150 (MQ=255)
                                                                                                                                              tCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGgaagaa         >  1:2070200/1‑150 (MQ=255)
                                                                                                                                              tCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGgaagaa         >  2:8505412/1‑150 (MQ=255)
                                                                                                                                                 aCCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGa      >  1:2219604/1‑150 (MQ=255)
                                                                                                                                                  cccTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGAt     <  2:9284638/150‑1 (MQ=255)
                                                                                                                                                  cccTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGAt     <  2:1652703/150‑1 (MQ=255)
                                                                                                                                                  cccTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGAt     >  1:8596112/1‑150 (MQ=255)
                                                                                                                                                     tGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTc  <  1:3591534/150‑1 (MQ=255)
                                                                                                                                                     tGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTc  <  1:4530685/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: