Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,720,364 G→T 9.6% intergenic (+21/+26) slyB → / ← slyA outer membrane lipoprotein/global transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,720,3640GT 9.6% 534.4 / 8.3 178intergenic (+21/+26)slyB/slyAouter membrane lipoprotein/global transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (92/61);  new base T (2/23);  total (94/84)
Fisher's exact test for biased strand distribution p-value = 6.02e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 185 total aligned reads displayed.

CTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCT  >  NC_000913/1720215‑1720513
                                                                                                                                                       |                                                                                                                                                     
cTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAgg                                                                                                                                                       >  2:165587/1‑150 (MQ=255)
 tGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGt                                                                                                                                                      >  2:3614590/1‑150 (MQ=255)
  ggAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGtt                                                                                                                                                     >  1:918369/1‑150 (MQ=255)
   gAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTa                                                                                                                                                    >  1:2905601/1‑150 (MQ=255)
     aaTTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACt                                                                                                                                                  <  1:3339053/150‑1 (MQ=255)
     aaTTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACt                                                                                                                                                  <  1:826284/150‑1 (MQ=255)
        tCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGAc                                                                                                                                               <  1:316141/150‑1 (MQ=255)
        tCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGAc                                                                                                                                               <  1:316140/150‑1 (MQ=255)
        tCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGAc                                                                                                                                               <  1:3063084/150‑1 (MQ=255)
          gTAAAGACGTTGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTGTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCTGGTTACTGTCca                                                                                                                                             <  2:1257728/150‑1 (MQ=255)
             aaGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCacac                                                                                                                                          >  1:3819264/1‑150 (MQ=255)
             aaGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCacac                                                                                                                                          <  2:500793/150‑1 (MQ=255)
             aaGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCacac                                                                                                                                          <  2:500787/150‑1 (MQ=255)
               gACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGc                                                                                                                                        >  1:72645/1‑150 (MQ=255)
               gACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGc                                                                                                                                        >  1:2301034/1‑150 (MQ=255)
                 ctttgttAATACAATTTTGGTGTTACAGAAACATGGCAACACGGTTTTCTCTCGGGGCCAACGTGTCGTACTGGCCAGCATTGGCTTTCTGGTGACGTTTTCTCCGCGATAAATATTTTTGGGTGTGGTCTGTTTACTCACCACACGccc                                                                                                                                      <  2:604988/144‑1 (MQ=255)
                             catcatGGTGTTACTGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCTGTCTGGTGACGTTTTCTCCGCGCTAAATAAGTTTGCGTTTGGTCTGGTTACTGACCACACGCCCCCTTCATTTCAc                                                                                                                          <  2:3092542/150‑1 (MQ=255)
                             catcatGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCAc                                                                                                                          >  1:2462809/1‑150 (MQ=255)
                             aatattggtgtTACAGAAACAAGGCAACACGCTTTTCTCTCGGGGCAACGGTGTGGTACTGGCCTGAAAGGGCTGCCTGGTGCCCTTTTCTCCGCGCTAAATATTTTTGGGTTTGTTCTGTTTACTGACCACACGCCCCCTTCATTTCAc                                                                                                                          <  2:2789582/145‑1 (MQ=255)
                                     tggtACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGc                                                                                                                  >  2:689709/1‑150 (MQ=255)
                                     tggtACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGc                                                                                                                  >  1:595694/1‑150 (MQ=255)
                                      ggtACATAACCAAGGAAACACTCTTTTCACTCGGGGCAAAGGTGTGGTACTGGCCAGCAAGGGCTTTCTGGGGCCGGTTTCGCCGCGCTAAATAAGTTTGTGTGTGGTCTGTTTCGTGACCACACGCCCCCTTCATTTCACCCTTTGGcc                                                                                                                 <  2:1915585/150‑1 (MQ=255)
                                       gtACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCt                                                                                                                >  1:1830186/1‑150 (MQ=255)
                                        tACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTg                                                                                                               >  1:2463177/1‑150 (MQ=255)
                                         aCAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGt                                                                                                              >  1:3106427/1‑150 (MQ=255)
                                          cAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGTCCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACGACTCGGCCGCTTCATTTCTCCCTTTGGCCTGTa                                                                                                             >  1:3077458/1‑150 (MQ=255)
                                          cAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTa                                                                                                             >  1:2892091/1‑150 (MQ=255)
                                          cAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTa                                                                                                             >  1:4004294/1‑150 (MQ=255)
                                          cAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTa                                                                                                             >  2:1770989/1‑150 (MQ=255)
                                           aGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTaa                                                                                                            >  1:2161194/1‑150 (MQ=255)
                                             aaaCAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACt                                                                                                          >  1:911566/1‑150 (MQ=255)
                                                cAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCaa                                                                                                       >  1:1277034/1‑150 (MQ=255)
                                                 aaGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGAGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCAAATCACCCTTACGCCAGTAACTCAAt                                                                                                      >  2:1257543/1‑150 (MQ=255)
                                                   ggCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTTCGTGTGGTCAGGTTAATGACCACAAGACACCTTCATTTCACACTTTGGCCTGTAACTCAATGa                                                                                                    >  1:4022055/1‑150 (MQ=255)
                                                    gCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGat                                                                                                   >  1:743649/1‑150 (MQ=255)
                                                     cAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGa                                                                                                                                                >  2:3988739/1‑104 (MQ=255)
                                                        caggggTTTCTCTCGGGGCCAACGTGTCTTAGTGGCAAGCTTTGGCTGTCTGGTGACCGTTTCTCCGCGCTAATTATGTTTGCTTGTTTTCTTGTTCCGGACCCCACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGAtatta                                                                                               <  2:802748/145‑1 (MQ=255)
                                                        cacgggTTTCTCTCGGGGCCACGGTGTCGTCGTGGCCAGCACTGGCAGTCGGGTGCCGTGTTCTCGGCGATAAATAATTTTGGGTTTGGGCTGTTTACTCACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGAtatta                                                                                               <  2:1514191/145‑1 (MQ=255)
                                                         acTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTCACTCAATGAtattat                                                                                              >  2:3323062/1‑150 (MQ=255)
                                                           tCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGc                                                                                            >  1:2220153/1‑150 (MQ=255)
                                                           tCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCAGTTCACCCTTTGGCCTGTAACTCAATGATATTATGc                                                                                            >  2:2717242/1‑150 (MQ=255)
                                                               ttCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCaa                                                                                        >  1:341816/1‑150 (MQ=255)
                                                               ttCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCaa                                                                                        >  1:342098/1‑150 (MQ=255)
                                                                    tcCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGtttt                                                                                   >  2:810073/1‑150 (MQ=255)
                                                                    tcCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGAGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGattt                                                                                   >  1:2488554/1‑150 (MQ=255)
                                                                     cggggCCAAGGTGTCGTAGTGGCCAGCATTGGATGTCTGGTGACCGTTTCTCCGCGATAAATACGTTTGGGTGTGGTCTGTTTACTCACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTg                                                                                  <  2:87552/148‑1 (MQ=255)
                                                                      cGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGc                                                                                 >  1:3979164/1‑150 (MQ=255)
                                                                      cGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGc                                                                                 >  2:4099857/1‑150 (MQ=255)
                                                                      cGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACAACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGATCAAGTTTTGc                                                                                 >  2:4100109/1‑150 (MQ=255)
                                                                       gggCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCg                                                                                >  2:484382/1‑150 (MQ=255)
                                                                         gCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACAGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATATCACCCTTTGGCCTGTAAATCAAAGATATTATGCTCAAGTTTTGCGAt                                                                              >  1:1307717/1‑150 (MQ=255)
                                                                                tgtCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTa                                                                       >  1:3038706/1‑150 (MQ=255)
                                                                                  tCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAAt                                                                     >  2:2366993/1‑150 (MQ=255)
                                                                                   cttcggggCCCGCACGGGCTGTCCGTTGACGTTTTCGCCGCGATAAATTTTTTGGGGTGTGTTCTGTTTACTCACCCCACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATc                                                                    <  2:1086675/144‑1 (MQ=255)
                                                                                   cggacgggCCAGCAATGGCAGTCCGGTGACCGTTTCTCCGCGCTAAATCAGTTTGCGTGTGGGCTGGCTACTGACCACACGCCCCCTTCATTTCCCCCTTTGGCCTGTAACTCACTGATATAATGCTCAAGTTTTGCGATGAGCGTAATc                                                                    <  2:1735052/144‑1 (MQ=255)
                                                                                     ttcTGGCCAGCATTGGCTGTCAGGTGACCGTTTCTCCGGGCTATATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAg                                                                  <  2:2367574/148‑1 (MQ=255)
                                                                                        tGGCCACCATTGGCTGTCCGGTGACCGTTTCTCCGCTATAAATATTTTGGGGTGTGTCCAGGTTACTCACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTg                                                               <  1:3189338/150‑1 (MQ=255)
                                                                                         ggCCAGCATGGGCAGTCTGGTGACCGTTTCGCCGCGATAAATATGTTTGGGTGTGGTCTGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGc                                                              <  2:2161194/150‑1 (MQ=255)
                                                                                           ccAGAATTGGCTGTCTGGGGACCTTTTCTCCGCGCTAAATAAGTTTGCTTGTGTTCTGGTAACTCACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTc                                                            <  2:2905601/150‑1 (MQ=255)
                                                                                           caagaaaTGGCTGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGGCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTc                                                            <  2:1305129/145‑1 (MQ=255)
                                                                                            cAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTcc                                                           >  2:3444209/1‑150 (MQ=255)
                                                                                                     ctgTCGGGTCACCTTTTCCCCGCTCTAAATTCTTTGGGGTGTGGTATGGTTACTCCCCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCctct                                                  <  1:1203242/148‑1 (MQ=255)
                                                                                                      tgTCTGGTGCCCGTTTCCCCGCTATAAATATGTTTGCGTGTGTTCTGTTTCCTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTg                                                 <  1:1335071/149‑1 (MQ=255)
                                                                                                       gTCTGGTCACGGTTTCTCCGCCAAAATATTGTTGGGGTGTGGTCTGGTTACTCACCCCACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGc                                                <  1:3834933/150‑1 (MQ=255)
                                                                                                        tctggTGACCTTTTCCCCGCGCTAAATATGTTGGCGTGTGGTCTGTTTACTCACCCCACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCg                                               <  1:2800591/147‑1 (MQ=255)
                                                                                                        tCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCg                                               <  2:341863/150‑1 (MQ=255)
                                                                                                         cAGGGGACCGTTTCTCCGCGATAAATATGTTTGCGTGTGGTCTGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCgg                                              <  2:2874541/150‑1 (MQ=255)
                                                                                                            gTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGGTCCTCTGCGgagc                                           >  2:2601879/1‑149 (MQ=255)
                                                                                                               aCCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGc                                        >  1:257496/1‑150 (MQ=255)
                                                                                                                cggTTTATCCGGGCTAAATAAGTTTGGGTGGTGTCATTTTACTTACCACTCGCCCCCTTCATTTCACCCTTTGGCATGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGcc                                       <  2:3671056/148‑1 (MQ=255)
                                                                                                                 cGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGCTGCCa                                      >  2:3949924/1‑150 (MQ=255)
                                                                                                                  ctttcccccGCGAAAAATCAGTTTGGGTGGGGCCTGGTAACGGACCCCACGCCCCCTTCATGTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCa                                      <  1:1613111/143‑1 (MQ=255)
                                                                                                                   tttCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATg                                    >  1:2534630/1‑150 (MQ=255)
                                                                                                                     ttCCCCGCGCTAAATATGTTGGCGTGTGGTCTGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGt                                   <  1:4113621/150‑1 (MQ=255)
                                                                                                                      tctcCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTa                                  <  2:2201483/150‑1 (MQ=255)
                                                                                                                      tctcCGCGATAAATTTGTTTGTGTGTGGTCTGGTTACTCACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTa                                  <  1:165587/150‑1 (MQ=255)
                                                                                                                        tcCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAat                                >  2:3819579/1‑150 (MQ=255)
                                                                                                                          cgcgCTAATTATTTTTGCGTGTGGTCAGTTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAatat                              <  2:2462809/150‑1 (MQ=255)
                                                                                                                            cgcTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCCTGTAATAttt                            >  2:3797652/1‑150 (MQ=255)
                                                                                                                            cgcTAAATAAGTTTGCGTGTGGTCAGGTTAATGACCACACGCCCCCTTCATTTCACCCATTGGCCTGTAAATCAATGATATTATGCTCAAGTTTTGCGATGTGCGGAATAAGTTGCTCCAGTTCCTATGCGGGGATGCCCTGGACTAttt                            >  1:3910513/1‑150 (MQ=255)
                                                                                                                              atAAATAAGTTGGCGTGTGGTCTGGTTAATGACAACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTcc                          <  2:595694/149‑1 (MQ=255)
                                                                                                                                  tttttgTTTGGGTGGGTTCTGTTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCcgcg                       <  2:1431218/145‑1 (MQ=255)
                                                                                                                                  aTAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCcgcgc                      >  2:3495063/1‑150 (MQ=255)
                                                                                                                                    ttaGTTGGGGTGTGGTCCGGTTCCTCACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCcgcgcg                     <  2:3140885/148‑1 (MQ=255)
                                                                                                                                    tatgTTTGCGTGTGGTCTGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCcgcgcg                     <  1:3442005/147‑1 (MQ=255)
                                                                                                                                     atttttGCGTGGGTTCTGTTTACGGACCCCACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCgg                    <  1:67966/147‑1 (MQ=255)
                                                                                                                                         ttGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGttt                <  2:2646373/150‑1 (MQ=255)
                                                                                                                                         tggggtgtGGTCTGTTTACTCACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGttt                <  1:1929931/146‑1 (MQ=255)
                                                                                                                                           gCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTg              <  2:3414468/150‑1 (MQ=255)
                                                                                                                                            cGTGTGGTCAGTTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGt             <  2:233385/150‑1 (MQ=255)
                                                                                                                                              tgtgGTCAGTAACCTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTta           >  2:1201938/1‑150 (MQ=255)
                                                                                                                                               gtgGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTtaa          >  1:3740326/1‑150 (MQ=255)
                                                                                                                                                ttgacAGGATACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTtaat         <  2:944216/146‑1 (MQ=255)
                                                                                                                                                tgGTCTGTTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTtaat         <  1:3365626/150‑1 (MQ=255)
                                                                                                                                                 ggTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTtaata        >  1:2677331/1‑150 (MQ=255)
                                                                                                                                                    cAGTTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACa     <  1:1018121/150‑1 (MQ=255)
                                                                                                                                                    cAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACa     >  2:2020891/1‑150 (MQ=255)
                                                                                                                                                      ggTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGc   >  2:1809035/1‑150 (MQ=255)
                                                                                                                                                      agTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGc   >  1:2013401/2‑150 (MQ=255)
                                                                                                                                                       gTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCt  >  1:362327/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
CTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCT  >  NC_000913/1720215‑1720513

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