Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,720,364 G→T 8.2% 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 8.2% 871.0 / 15.1 275intergenic (+21/+26)slyB/slyAouter membrane lipoprotein/global transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (168/73);  new base T (3/31);  total (171/104)
Fisher's exact test for biased strand distribution p-value = 6.85e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 290 total aligned reads displayed.

CTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCT  >  NC_000913/1720215‑1720513
                                                                                                                                                     |                                                                                                                                                     
cTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAgg                                                                                                                                                       <  2:5875685/150‑1 (MQ=255)
 tGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGt                                                                                                                                                      <  1:4637347/150‑1 (MQ=255)
  ggAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGtt                                                                                                                                                     <  1:6858536/150‑1 (MQ=255)
   gAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTCCTGGCCCGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTa                                                                                                                                                    >  1:5151453/1‑150 (MQ=255)
   gAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTa                                                                                                                                                    >  1:6326680/1‑150 (MQ=255)
   gAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTa                                                                                                                                                    >  1:859030/1‑150 (MQ=255)
    aaaTTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTAc                                                                                                                                                   <  1:6514694/150‑1 (MQ=255)
     aaTTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACt                                                                                                                                                  <  2:1344666/150‑1 (MQ=255)
       ttCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGATCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGa                                                                                                                                                <  2:707134/150‑1 (MQ=255)
       ttCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGa                                                                                                                                                <  1:5582657/150‑1 (MQ=255)
       ttCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGa                                                                                                                                                <  2:6635425/150‑1 (MQ=255)
       ttCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGa                                                                                                                                                >  1:3460449/1‑150 (MQ=255)
       ttCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGa                                                                                                                                                <  2:1530148/150‑1 (MQ=255)
         cGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGACAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGAcc                                                                                                                                              <  1:1472950/150‑1 (MQ=255)
           tAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGTAACCTGACcac                                                                                                                                            <  2:26418/150‑1 (MQ=255)
               gACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGc                                                                                                                                        >  2:3350143/1‑150 (MQ=255)
                   aTGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGccccc                                                                                                                                    >  1:494767/1‑150 (MQ=255)
                     ggTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCtt                                                                                                                                  >  1:6078397/1‑150 (MQ=255)
                       tAATACCATCTTGGTGTTACAGAAACAAGGCAACACTCGTTTCTCTCGGGGCAAACGTGTGGTACTGGCCTGCAATGGCTGTCTGGTGACCGTTTCTCCGCGCTAAATATTTTTGCGTGTGTTCTGGTTAGTGACCACACGCCCCCTTCa                                                                                                                                <  2:5445586/150‑1 (MQ=255)
                       tAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTCCTGGCCAGCAATGGCAGTCCGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCa                                                                                                                                >  1:3241050/1‑150 (MQ=255)
                       ataTACAATCTTGGTGTTACATAAACAAGGCAACACTCGTTTCTCCCCGGGCGAACGTGTCGTACTGGCCAGCAATGGCTGTCAGTTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCa                                                                                                                                <  2:4666509/148‑1 (MQ=255)
                          tACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCAttt                                                                                                                             >  2:1961960/1‑150 (MQ=255)
                          aaacatcttGGTGTAACATAAACAAGCCAACCCTCGTTTCTCTCCGGGCCAAGGTGTAGTCGTGGCCAGCATTGTATTTGCGTTCACGTTTTCGCCGCGATATATTTTTTTGCGTGTGGTCATTTTACTCACCACACGCCCCCTTCAttt                                                                                                                             <  2:2534612/147‑1 (MQ=255)
                           aCCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTc                                                                                                                            >  2:3684186/1‑150 (MQ=255)
                             catcatGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCAc                                                                                                                          >  1:425436/1‑150 (MQ=255)
                                 atGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCtt                                                                                                                      >  1:3867149/1‑150 (MQ=255)
                                 atGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCtt                                                                                                                      >  1:1518258/1‑150 (MQ=255)
                                  tggtggTACAGAAACACGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCCGGTGACCGTTTCTCCGCGCTAAATAAGTTTGGGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCttt                                                                                                                     <  2:4100443/150‑1 (MQ=255)
                                  tggtggTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCttt                                                                                                                     <  2:6008445/150‑1 (MQ=255)
                                  tggtggAACAGAAACCCGGCAACACTGGTTTCTCTCGGGGCCAACGTGTGGTACTGGCCAGCATTGGCTGTCTGGTGCCCTTTTCTCGGCGCTAAATATTTTGGGGTGTGGTCTGGTTACTGACCACACGCCCCCTTCATTTCACCCttt                                                                                                                     <  2:100768/150‑1 (MQ=255)
                                         aCAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGt                                                                                                              >  1:2812198/1‑150 (MQ=255)
                                         aCAGAAACAAGGCAACACTCGTATCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGATACTCCTCGCTAAATAAGTTAGCGTGTGGTCAGGTTACTGACCAAACGCCCCCTGCATTTCACCCATTGGCCTGt                                                                                                              >  1:6744136/1‑150 (MQ=255)
                                          cAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACTCCCCCTTCATTTCACCCTTTGGCCTGTa                                                                                                             >  1:1006725/1‑150 (MQ=255)
                                          cAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTa                                                                                                             >  1:4716101/1‑150 (MQ=255)
                                          cAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTa                                                                                                             >  2:2929346/1‑150 (MQ=255)
                                          cAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTACATTTCACCCTTTGGCCTGTa                                                                                                             >  2:4529100/1‑150 (MQ=255)
                                           aGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTaa                                                                                                            >  2:1652826/1‑150 (MQ=255)
                                           aGAAAAAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCACTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCCTCAATTAACCCATTTGCCTGTaa                                                                                                            >  2:1513126/1‑150 (MQ=255)
                                                     cAACACTGTTTTCTCTCCGGGCCAACGTGTCTTAGGGGCCAGCACGGGCTGTCTGGTGACCGTTTCTCCGCGCTAAATATGTTTGGGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGata                                                                                                  <  2:4986906/150‑1 (MQ=255)
                                                          cTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAAGGGCTGTCAGGTGACCGTTTCTCCGCGATAAATAATTTTGCGTGTGGTCTGTTTACTTACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATg                                                                                             <  2:6370306/150‑1 (MQ=255)
                                                              tttCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCa                                                                                         >  1:2020065/1‑150 (MQ=255)
                                                               ttCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCaa                                                                                        >  2:5855710/1‑150 (MQ=255)
                                                                tctctcCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAg                                                                                       >  2:6888586/1‑150 (MQ=255)
                                                                    tcCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGtttt                                                                                   >  2:5618267/1‑150 (MQ=255)
                                                                    tcCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGtttt                                                                                   >  2:5947727/1‑150 (MQ=255)
                                                                     ccGGCCCAACGTGTCGAACTGCCCAGCAATGGCAGACAGGACACCGTTACACCCCGATAAAAAAGTTTGCGTGTGGACAGGTTACTCACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTg                                                                                  <  2:2391071/150‑1 (MQ=255)
                                                                         gCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATCTTATGCTCAAGTTTTGCGAt                                                                              >  1:5796970/1‑150 (MQ=255)
                                                                          ccAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATg                                                                             >  1:5825850/1‑150 (MQ=255)
                                                                           cAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGa                                                                            >  1:3012887/1‑150 (MQ=255)
                                                                           cAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCATTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGa                                                                            >  2:1267123/1‑150 (MQ=255)
                                                                           cAACGTGTCGTACTGGCCAGCAATGGCAGGCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTACCTCAATGATATTATGCTCCAGTTTTGCGATGa                                                                            >  2:3768773/1‑150 (MQ=255)
                                                                              cGTGTCTTTCGGGCCAGCATGGGCAGTCCGGTGACGTTTTCCCCGCGCTAAATATTTGTGGGTGTGGTCTGTTGCCTGCCCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCg                                                                         <  2:6114007/150‑1 (MQ=255)
                                                                                  tCGTAGTGGCCAGCTTTGGCAGTCTGGTGACCGTTGCTCCGCGCTAAATAAGTTTGCGTGTGGTCTGGTTACTGACCACAGGCCACCTTCATTTCACCCTTTGGCCTGTAACTCAATTATATTATGCTCAAGTTTTGCGATGAGCGTAAt                                                                     <  2:5390311/150‑1 (MQ=255)
                                                                                  tCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAAt                                                                     >  2:1212966/1‑150 (MQ=255)
                                                                                  tCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAAt                                                                     >  2:4328218/1‑150 (MQ=255)
                                                                                     tACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGTCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAg                                                                  >  2:508810/1‑150 (MQ=255)
                                                                                        tGGCCTGCATTGGCAGTCTGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCTGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTg                                                               <  2:3030275/150‑1 (MQ=255)
                                                                                        tGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTg                                                               >  1:371756/1‑150 (MQ=255)
                                                                                        tGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTg                                                               >  1:5526230/1‑150 (MQ=255)
                                                                                         ggCCAGCATTGGCTGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTGGCGTGTGGTCAGTTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGc                                                              <  2:6557571/150‑1 (MQ=255)
                                                                                          gCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCt                                                             >  2:4185261/1‑150 (MQ=255)
                                                                                           ccAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTc                                                            >  1:5202252/1‑150 (MQ=255)
                                                                                              gCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAg                                                         >  2:4386868/1‑150 (MQ=255)
                                                                                                aaTGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGtt                                                       >  1:5635227/1‑150 (MQ=255)
                                                                                                     cAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCctct                                                  >  2:4865859/1‑150 (MQ=255)
                                                                                                        tCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCg                                               >  1:6982864/1‑150 (MQ=255)
                                                                                                        tCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCATTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCg                                               >  1:36371/1‑150 (MQ=255)
                                                                                                        tCAGGAGACCGTTTCTACGAGCTAAATAAGTTTCCGTGTGCTCAGGTTACTGACCACACCCCCCCTTCATTTCACCATGTGGCCTGTAACTAAATGATATTATGCTCAAGTTTAGCGATGAGCGTCATCAGTTGCTCCAGTACCTCTGCg                                               >  2:3465711/1‑150 (MQ=255)
                                                                                                          tggTGACCTTTTCTCCGCGCTAAATAAGTTTGCGTGGGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGga                                             <  1:1079476/149‑1 (MQ=255)
                                                                                                          tggTGACCGTTGCGCCGCGATAATTATGTTTGGGTGTGGTCTGTTTACTCACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGga                                             <  1:6186792/149‑1 (MQ=255)
                                                                                                          tggTGACCGTTACCCCGCTATAAATTAGTTTGCTTGTGGACTGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGga                                             <  1:1079438/149‑1 (MQ=255)
                                                                                                          aGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCGTTGGCCTGTAACTAAATGATATTATGCTCAAGGTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCggt                                             >  1:2493896/1‑149 (MQ=255)
                                                                                                            gTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGgaga                                           >  2:3870102/1‑150 (MQ=255)
                                                                                                            gTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGgaga                                           >  1:3570869/1‑150 (MQ=255)
                                                                                                             tgacgttttCGCCGCTATAAATATGTTTGGGTGTGTTCTGTTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGAt                                          <  1:1698814/144‑1 (MQ=255)
                                                                                                                  gTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCGTTGGCCTGTTACTCAATGTTATTTTGCTAAAGTTTTGCGATGAGCGTAAGCAGTTGCTCCAGGTCCACTGCGGTGATGCCAt                                     >  1:7210760/1‑150 (MQ=255)
                                                                                                                    ttCTCGGCTCTAAATATGTTTGGTTGTGTTCAGTTTACTCACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGt                                   <  1:188910/150‑1 (MQ=255)
                                                                                                                       tcCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCGTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGGGATGAGCGTAATCACTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAat                                >  2:5226678/1‑150 (MQ=255)
                                                                                                                         cgcgCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGCTGCCATGTAatat                              >  2:5737051/1‑150 (MQ=255)
                                                                                                                         cgcgCTAAATAAGTTTGCGTGTGGTCAGGGTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTGTTGCGATGAGCGGAAGCAGTTGCTCCAGTGCCTCTGCGGAGATGCCATggcctag                              >  2:6174951/1‑144 (MQ=255)
                                                                                                                           cgcTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATAttt                            <  2:3548434/150‑1 (MQ=255)
                                                                                                                                   aaGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTATGCGGAGATGCCATGTAATATTTCCGCGCgg                    >  1:4411887/1‑150 (MQ=255)
                                                                                                                                      tttGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATAAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGtt                 >  2:460978/1‑150 (MQ=255)
                                                                                                                                         gCGTGTGGTCTTGTTACTCCCCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTg              <  1:350031/150‑1 (MQ=255)
                                                                                                                                          cGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGt             <  2:2985862/150‑1 (MQ=255)
                                                                                                                                           gtgtGGTCTGTTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGtt            <  2:1936485/150‑1 (MQ=255)
                                                                                                                                           gtgtGGTCAGTAACCTGACCACTCGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGtt            >  2:2163170/1‑150 (MQ=255)
                                                                                                                                           gtggGGTCTGGTTACTGCCCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGtt            <  1:4497219/150‑1 (MQ=255)
                                                                                                                                            tgtgGTCAGTAACCTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTta           >  2:5398113/1‑150 (MQ=255)
                                                                                                                                            tgtgGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATGTCCGCGCGGGTTTTGTta           >  1:5226881/1‑150 (MQ=255)
                                                                                                                                            tggggTCAGGTTACTGACCCCACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTta           <  2:6413037/147‑1 (MQ=255)
                                                                                                                                             gtgGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTtaa          <  2:6033920/150‑1 (MQ=255)
                                                                                                                                              tgttctggTTACTCACCACACGCCCCCTTCATGTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTtaat         <  1:1656932/144‑1 (MQ=255)
                                                                                                                                                gTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTtaataa       >  1:2331732/1‑150 (MQ=255)
                                                                                                                                                gTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTtaataa       >  1:2729675/1‑150 (MQ=255)
                                                                                                                                                 tCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAAc      >  2:1739636/1‑150 (MQ=255)
                                                                                                                                                   ttgTTCCTTACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTTATATTTCCGCGCGGGTTTTGTTAATAACAg    <  2:4440613/148‑1 (MQ=255)
                                                                                                                                                     gTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCt  >  2:2764752/1‑150 (MQ=255)
                                                                                                                                                     gTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCt  >  1:5869415/1‑150 (MQ=255)
                                                                                                                                                     gTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGAAACACAAAGATATTATGATAAAGTTTTGCGATGAGCGTAATAAGTTGCTCCAGTTCCTCTGCTGAGATGCCATGTAATATTTCAGCGCCGGTTTTGTTAATAACAgat  >  1:5382200/1‑148 (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: