Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,720,364 G→T 7.5% 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 7.5% 1339.5 / 12.9 402intergenic (+21/+26)slyB/slyAouter membrane lipoprotein/global transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (264/90);  new base T (3/45);  total (267/135)
Fisher's exact test for biased strand distribution p-value = 1.11e-20
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 419 total aligned reads displayed.

CTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCT  >  NC_000913/1720215‑1720513
                                                                                                                                                     |                                                                                                                                                     
  ggAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGtt                                                                                                                                                     >  1:10068670/1‑150 (MQ=255)
    aaaTTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTAc                                                                                                                                                   <  2:9382277/150‑1 (MQ=255)
       ttCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGa                                                                                                                                                <  1:512650/150‑1 (MQ=255)
        tCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGAc                                                                                                                                               >  1:5231180/1‑150 (MQ=255)
        tCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGAc                                                                                                                                               >  1:9474078/1‑150 (MQ=255)
        tCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGAc                                                                                                                                               <  1:10577721/150‑1 (MQ=255)
        tCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGAc                                                                                                                                               <  1:8259877/150‑1 (MQ=255)
          gTACAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACca                                                                                                                                             <  2:2953084/150‑1 (MQ=255)
          gTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACca                                                                                                                                             <  1:8778127/150‑1 (MQ=255)
            aaaGACGATGGAAATACCATCATGGTGGTACAGAAACAAGCCAACACTCGATACTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACcaca                                                                                                                                           <  2:5219297/150‑1 (MQ=255)
              aGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACg                                                                                                                                         <  2:7418045/150‑1 (MQ=255)
                  gATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTCAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGcccc                                                                                                                                     >  1:4383199/1‑150 (MQ=255)
                  gATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTCCTGACCACACGcccc                                                                                                                                     >  1:2645514/1‑150 (MQ=255)
                  gATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGcccc                                                                                                                                     >  2:10671274/1‑150 (MQ=255)
                  gATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGcccc                                                                                                                                     >  1:3703940/1‑150 (MQ=255)
                   aTGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGccccc                                                                                                                                    >  2:10737959/1‑150 (MQ=255)
                    tGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCt                                                                                                                                   >  1:4842354/1‑150 (MQ=255)
                       tAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCa                                                                                                                                >  1:3376777/1‑150 (MQ=255)
                       tAATACCATCATGGTGGTACAGAAAAAAGGAAACACTCGTTTCTCTCCGGGCAAAAGTGTCGTACTGGCCAGCAATGGCAGTCAGTTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGTTTACTGACCACACGCCCCCTACa                                                                                                                                >  1:8077095/1‑150 (MQ=255)
                        aaTACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCAt                                                                                                                               >  1:5803804/1‑150 (MQ=255)
                         aTACCTTCATGGTGGTACAGAACCACGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCTTTTCTCCGCGCTAAATACGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCAtt                                                                                                                              <  2:6360753/150‑1 (MQ=255)
                             catcatGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCAc                                                                                                                          >  2:2266457/1‑150 (MQ=255)
                             catcatGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCAc                                                                                                                          >  2:2164920/1‑150 (MQ=255)
                                    gtggtACATAAACTAGCCAACACTCGTTTCTCTCGGGGCCACGTTGTCGTACGGGCCTGCAATGGCTTTCAGGTGAGGTTTTGTCGGCGCTATTTAATTTTGGGTGTGGTCTGTTTACTGACCACACGCCCCCTTCATTTCACCCTTTgg                                                                                                                   <  2:4898224/150‑1 (MQ=255)
                                    ggggtACATAAACACGGCAACCCGGGTTTCTCTCCGGGCCACGGTGTCTTAGTGGCCAGCATTGGCTGGCTGGTGACCGTTTCTCCGCGATAAATAAGTTTGGGTGTGTTCAGGTTACGGACCACACGCCCCCTTCATTTCACCCTTTgg                                                                                                                   <  2:8511027/148‑1 (MQ=255)
                                     tggtACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGc                                                                                                                  >  2:8897004/1‑150 (MQ=255)
                                           aGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTaa                                                                                                            >  1:1449354/1‑150 (MQ=255)
                                            gAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAAc                                                                                                           >  2:8277481/1‑150 (MQ=255)
                                              aaCAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTAGTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTc                                                                                                         >  1:10014827/1‑150 (MQ=255)
                                               aCAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCa                                                                                                        >  1:10281326/1‑150 (MQ=255)
                                                cAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCaa                                                                                                       >  2:1619014/1‑150 (MQ=255)
                                                 acgcaaaCACTCGTTTCTCTCCGGGCAACGGTGTGGTAGTGGCCAGCATGGGCAGTCTGGTGAGGGTTTCTCGGCGATAAATTCGTTTGGGTGTGTTCTGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAAt                                                                                                      <  2:7851858/145‑1 (MQ=255)
                                                 aaGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAAt                                                                                                      >  1:3123206/1‑150 (MQ=255)
                                                 aaGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAAt                                                                                                      >  2:7700478/1‑150 (MQ=255)
                                                   ggCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGa                                                                                                    >  1:3894044/1‑150 (MQ=255)
                                                         tctggTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCTGTCAGGTGACCGTTTCTCGGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACTCGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGAtattat                                                                                              >  2:9734135/5‑150 (MQ=255)
                                                         acTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGAtattat                                                                                              >  1:7364401/1‑150 (MQ=255)
                                                          cTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATTATATTATg                                                                                             >  1:5067289/1‑150 (MQ=255)
                                                             gTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTc                                                                                          >  1:7295200/1‑150 (MQ=255)
                                                             gTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTc                                                                                          >  2:1221607/1‑150 (MQ=255)
                                                             gTTTATCTCGGGGCCAACGTGTCGTCGTGGCCAGCATGGGCTGTCTGGTGACGGTTTCTCCGCGATAAATATGTTTGGGTGTGTTCTGTTTACTGACAACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTc                                                                                          <  2:4468074/150‑1 (MQ=255)
                                                              tttCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCa                                                                                         >  1:2331976/1‑150 (MQ=255)
                                                              tttCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCa                                                                                         >  1:4623079/1‑150 (MQ=255)
                                                               ttCTCTCCGGGCCACGGTGTCGTAGTGGCCAGCACTGGCTGTCGGGTGACCTTTTCACCGCGATAAATATGTTTGCGTTTGTTCTGTTTACTCACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCaa                                                                                        <  2:3046305/150‑1 (MQ=255)
                                                                 ctctcCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGt                                                                                      >  2:6784979/1‑150 (MQ=255)
                                                                  tctcCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGtt                                                                                     >  1:6983817/1‑150 (MQ=255)
                                                                    tcCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGtttt                                                                                   >  1:9877745/1‑150 (MQ=255)
                                                                    tcCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGtttt                                                                                   >  1:1820892/1‑150 (MQ=255)
                                                                                  tCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAAt                                                                     >  2:10902120/1‑150 (MQ=255)
                                                                                    gTACTGGCCAGCATTGGCTGTCAGGTGACCGTTTCTCCGCGATAAATAATTTTGGGTGTGGTCTGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCa                                                                   <  2:543781/150‑1 (MQ=255)
                                                                                    gTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCa                                                                   >  1:6338466/1‑150 (MQ=255)
                                                                                     tACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAg                                                                  >  2:7621244/1‑150 (MQ=255)
                                                                                     tACTGGCCAGCAATAGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAg                                                                  >  2:11136349/1‑150 (MQ=255)
                                                                                        tGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTg                                                               >  2:4108422/1‑150 (MQ=255)
                                                                                        tGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTg                                                               >  1:5371140/1‑150 (MQ=255)
                                                                                          gCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACTACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCt                                                             >  1:11047099/1‑150 (MQ=255)
                                                                                          gCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCt                                                             >  2:6160782/1‑150 (MQ=255)
                                                                                          gCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCt                                                             >  1:8513061/1‑150 (MQ=255)
                                                                                          gCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCt                                                             >  1:9687799/1‑150 (MQ=255)
                                                                                            cAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTTTGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTcc                                                           >  1:5592131/1‑150 (MQ=255)
                                                                                            cAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTcc                                                           >  2:1017742/1‑150 (MQ=255)
                                                                                             aGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCa                                                          >  2:419325/1‑150 (MQ=255)
                                                                                             aGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCa                                                          >  2:5604582/1‑150 (MQ=255)
                                                                                             aGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCa                                                          >  1:1301734/1‑150 (MQ=255)
                                                                                              gCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAg                                                         >  2:1945340/1‑150 (MQ=255)
                                                                                                aaTGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGtt                                                       >  2:1470767/1‑150 (MQ=255)
                                                                                                  tGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTcc                                                     >  2:6526221/1‑150 (MQ=255)
                                                                                                   ggCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCct                                                    >  1:8955410/1‑150 (MQ=255)
                                                                                                           ggTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGgag                                            >  1:1166551/1‑150 (MQ=255)
                                                                                                           ggTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGgag                                            >  2:6286953/1‑150 (MQ=255)
                                                                                                            gggACCGTTTCTCCGCGATAAATACGTTGGGGTGTGGACTTGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGgaga                                           <  1:930806/148‑1 (MQ=255)
                                                                                                            gTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGgaga                                           >  1:5997061/1‑150 (MQ=255)
                                                                                                            gTGACCGTTTCTCCGCGCTAAATAAGTTTGCGGGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGgaga                                           >  2:9003936/1‑150 (MQ=255)
                                                                                                               aCCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGc                                        >  2:3934379/1‑150 (MQ=255)
                                                                                                               aCCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGc                                        >  2:4957221/1‑150 (MQ=255)
                                                                                                               aCCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCGTTGGCCTGTAACTCAATGATATTATGCTCAAGGTTTGCGATGAGCGGAATCAGTTGCTCCAGTTCCTCTGCGGAGATGc                                        >  1:8857544/1‑150 (MQ=255)
                                                                                                                 ggttgcgccGCGCTAAATCTGTGGGGGGGGGGTCGGTTGACTCCCCCCACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCa                                      <  1:8628102/143‑1 (MQ=255)
                                                                                                                 cGTTTCCCCGCGATAAATAAGTTTGCGTGTGTACAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCa                                      <  1:10204036/150‑1 (MQ=255)
                                                                                                                  gTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCAt                                     >  2:11089082/1‑150 (MQ=255)
                                                                                                                   tttCTCCGCGCTAATTATGTGGGGGTGTGGTCTGTTTACTCACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATg                                    <  1:10737959/150‑1 (MQ=255)
                                                                                                                    gtctcCGCGATAATTTTGTTGGGGTGTGGTCTGTTTACCCACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCGAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGt                                   <  1:2164920/149‑1 (MQ=255)
                                                                                                                          gcgcTAAATAAGTTTGCGGGTGGTCAGGTTACTGACCACACGCCCCCATCAGTTCACCCATTGGCCTGTAACTCAAAGATATTATGCTAAAGTTTTGCGCTGAGCGTAATCAGGTGCGCCACTTCCTCTGCGGGGCTGCCAGGTAATAtt                             >  1:2906754/1‑150 (MQ=255)
                                                                                                                          gcgaTAAATATGTTTGCTTGTGGTCTGTTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATAtt                             <  1:3586399/150‑1 (MQ=255)
                                                                                                                            gcaaaaTAAGTTGGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTc                           <  1:6509333/147‑1 (MQ=255)
                                                                                                                             cTAAATTAGTGTGGGTGTGTTCCTTTTACTCACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTcc                          <  1:9614954/150‑1 (MQ=255)
                                                                                                                             cTAAATATGTTTGCGTGTGGTCTGGTTACGGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTcc                          <  2:5714781/150‑1 (MQ=255)
                                                                                                                              tAAATATGTTTGCTTGTGGTCTGTTTACTCACCCCACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCcg                         <  1:7759788/150‑1 (MQ=255)
                                                                                                                              tAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCcg                         <  2:696087/150‑1 (MQ=255)
                                                                                                                                     gTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGt                  >  2:2401274/1‑150 (MQ=255)
                                                                                                                                     gTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGt                  >  2:10418173/1‑150 (MQ=255)
                                                                                                                                          gttgtgGTCTGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGt             <  1:7473082/148‑1 (MQ=255)
                                                                                                                                          cGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGt             <  2:7295200/150‑1 (MQ=255)
                                                                                                                                           gtgtGGTCAGTAACCTGACCACACGCCCCCTTCATTTCACCCGTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGtt            >  2:4364128/1‑150 (MQ=255)
                                                                                                                                              tgGTCTGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTtaat         <  1:7051108/150‑1 (MQ=255)
                                                                                                                                              tgGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTtaat         >  1:10458700/1‑150 (MQ=255)
                                                                                                                                                gTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTtaataa       <  2:3703940/150‑1 (MQ=255)
                                                                                                                                                 tCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTtaat         <  1:4146058/147‑1 (MQ=255)
                                                                                                                                                 tCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAAc      >  2:2533182/1‑150 (MQ=255)
                                                                                                                                                     gTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCt                                                             >  1:6728960/1‑91 (MQ=255)
                                                                                                                                                     gTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCt  >  2:9533987/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 < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑