Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,720,364 G→T 7.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 7.2% 217.0 / 2.4 83intergenic (+21/+26)slyB/slyAouter membrane lipoprotein/global transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (56/19);  new base T (2/6);  total (58/25)
Fisher's exact test for biased strand distribution p-value = 8.16e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGC  >  NC_000913/1720219‑1720512
                                                                                                                                                 |                                                                                                                                                    
aaaTTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCTCGCTAAATAAGTTTGCGTGTGGTCAGGTTaa                                                                                                                                                  >  1:7041/1‑149 (MQ=255)
 aaTTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACt                                                                                                                                                 >  1:736468/1‑150 (MQ=255)
  aTTCGTAAAGACGATGGTAATACCATCATGGTTATACAGCACCAAGGCAACACTCGTTTATCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTg                                                                                                                                                <  2:1206996/150‑1 (MQ=255)
  aTTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGGTTCTCTCCGGGCCTACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCTCTAAATAACTTTGCGTGTGGTCAGGTTACTg                                                                                                                                                <  2:518238/150‑1 (MQ=255)
   ttCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGATTCTCTCCGGGCCAACGAGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGATTGTGGTCAGGTTACTGa                                                                                                                                               >  1:49044/1‑150 (MQ=255)
     cGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGAcc                                                                                                                                             <  1:5080/150‑1 (MQ=255)
     cGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTAAcc                                                                                                                                             >  2:1158144/1‑150 (MQ=255)
     cGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTAAGGTTAATGAcc                                                                                                                                             >  1:1115255/1‑150 (MQ=255)
      gTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGAccc                                                                                                                                            >  2:1387237/1‑149 (MQ=255)
       tAACGACGATGGTAATACCATCATGGTGGCACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGGACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACcac                                                                                                                                           <  2:866953/150‑1 (MQ=255)
       tAAAGACGACGGTAATACCATCATGGTGGTACAGAACCAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACcac                                                                                                                                           <  2:168792/150‑1 (MQ=255)
        aaaGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGGCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACcaca                                                                                                                                          <  2:350602/150‑1 (MQ=255)
         aaGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGTAACCTGACCacac                                                                                                                                         <  2:781299/150‑1 (MQ=255)
         aaGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCacac                                                                                                                                         <  1:1556736/150‑1 (MQ=255)
          aGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCCACGTGTCGTACTGGCCAGCCATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACg                                                                                                                                        >  1:333158/1‑150 (MQ=255)
          aGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACg                                                                                                                                        <  2:1266781/150‑1 (MQ=255)
          aGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGACAGGTGCCCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTTACCACACg                                                                                                                                        >  2:1129449/1‑150 (MQ=255)
                   tAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCa                                                                                                                               >  2:523435/1‑150 (MQ=255)
                     aTACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCAtt                                                                                                                             >  1:1399877/1‑150 (MQ=255)
                     aTACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCAtt                                                                                                                             >  1:495431/1‑150 (MQ=255)
                       aCCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGCCCACACGCCCCCGTCATTTc                                                                                                                           >  2:217690/1‑150 (MQ=255)
                        ccATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACCCGCCCCCTTCATTTCa                                                                                                                          >  1:194583/1‑150 (MQ=255)
                                 tggtACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGCCCACACGCCCCCTTCATTTCACCCTTGGGc                                                                                                                 >  2:385544/1‑150 (MQ=255)
                                  ggtACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCCCCCGCCCCCTTCATTTCACCCTTTGGcc                                                                                                                >  2:1399664/1‑150 (MQ=255)
                                   gtACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCt                                                                                                               >  1:1715216/1‑150 (MQ=255)
                                     aCAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTTCTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGt                                                                                                             >  2:1093961/1‑150 (MQ=255)
                                     aCAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGt                                                                                                             >  1:1718428/1‑150 (MQ=255)
                                     aCAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGt                                                                                                             >  2:20059/1‑150 (MQ=255)
                                     aCAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGt                                                                                                             >  2:610158/1‑150 (MQ=255)
                                     aCAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCATTTGGCCTGt                                                                                                             >  1:790593/1‑150 (MQ=255)
                                      cAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGCTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTgga                                                                                                            >  1:1180609/1‑148 (MQ=255)
                                        gAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCCTTTCACCCTTTGGCCTGTAAc                                                                                                          >  2:1553117/1‑150 (MQ=255)
                                            cAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTAACCCTTTGGCCTGTAACTCaa                                                                                                      >  2:601895/1‑150 (MQ=255)
                                              aGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCGTCATTTCACCCTTTGGCCTGTAACTCAATg                                                                                                    >  2:1609261/1‑150 (MQ=255)
                                                      cTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATg                                                                                            >  2:43084/1‑150 (MQ=255)
                                                              tctcCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGtt                                                                                    >  1:352260/1‑150 (MQ=255)
                                                               ctcCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGCCCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGGTACTGACGACACGCCCCCGTCATTTCACCCTTTGGCCTGTAACTCAATGAGATTATGCTCAAGgtt                                                                                   >  1:1158291/1‑150 (MQ=255)
                                                               ctcCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGttt                                                                                   >  2:1306001/1‑150 (MQ=255)
                                                                   gggCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCCCACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCg                                                                               >  1:1716801/1‑150 (MQ=255)
                                                                    ggCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCGATATAAGTTTGCGTGTGGTGAGGTTACTGACCACACGCCCCCTTCCTTTCACCCTATGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGa                                                                              >  1:171863/1‑150 (MQ=255)
                                                                        aaCGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTAACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGac                                                                          >  1:941476/1‑149 (MQ=255)
                                                                             gtCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGATAAATAAGTTTGCGTGTGGTCAGTTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTaa                                                                     <  2:7614/150‑1 (MQ=255)
                                                                              tCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCCCCCGCCCCCTTCATTTCACCCTTTGGCCTGTAACTGAATGATATTATGCTCAAGTTTTGCGATGCGCGTAAt                                                                    >  2:166800/1‑150 (MQ=255)
                                                                                   cTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCCCACGCCCCCGTCATTTCTCCCATTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGtt                                                               >  1:366022/1‑150 (MQ=255)
                                                                                   cTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGtt                                                               >  1:617001/1‑150 (MQ=255)
                                                                                     ggCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGc                                                             >  2:454779/1‑150 (MQ=255)
                                                                                      gCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCt                                                            >  2:1081407/1‑150 (MQ=255)
                                                                                       cctccAATGGAAGTCAGGTGACCGTTTCACCGCGCTAATTAAGTTTGGGTGGGGTCAGGTTACTCACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTc                                                           <  2:790593/146‑1 (MQ=255)
                                                                                       ccAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTc                                                           >  2:54252/1‑150 (MQ=255)
                                                                                          gCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAg                                                        >  1:1854625/1‑150 (MQ=255)
                                                                                          gCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAg                                                        >  2:236471/1‑150 (MQ=255)
                                                                                           cAATGGCAGTCAGGGGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACAACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGt                                                       <  1:236017/150‑1 (MQ=255)
                                                                                                gCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCctc                                                  <  2:110437/150‑1 (MQ=255)
                                                                                                   gTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGc                                               >  2:1641050/1‑150 (MQ=255)
                                                                                                       ggTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTACTCTCCGcag                                           >  1:402091/1‑150 (MQ=255)
                                                                                                          gACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGCCCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATg                                        >  2:926777/1‑150 (MQ=255)
                                                                                                          gACCGTTTATCCGCGCTAAATAAGTTTGCGGGTGGTCTGGTTACTTACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATg                                        <  2:1279957/150‑1 (MQ=255)
                                                                                                           aCCGTTTCTCGGCTGTATATCTGTTTGCGGGTGGGCTGTTTAGTTACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGc                                       <  2:1837170/150‑1 (MQ=255)
                                                                                                                ttCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTAGTGACCACACGACCCGTTCATTTCACCCTTTGGCCTGTATCTCTATGATATTATGCTCTAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGt                                  <  1:72187/150‑1 (MQ=255)
                                                                                                                   tcCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTAACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAat                               >  2:1831172/1‑150 (MQ=255)
                                                                                                                    ccGCGAAAAATAATTTTGCGTGTGGTCTGTTTACTCACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAata                              <  2:681451/150‑1 (MQ=255)
                                                                                                                      gcgcGAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTATGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATGAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATAtt                            >  2:257894/1‑150 (MQ=255)
                                                                                                                      ccgcAAAATCAGTTTGCGTGCGGTCAGGTTACTCACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATAtt                            <  2:1294862/149‑1 (MQ=255)
                                                                                                                       cgcTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATAttt                           >  1:1080186/1‑150 (MQ=255)
                                                                                                                        gcTAAATATGTTTGCGTGTGGTATGGTTACTTACAACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTc                          <  2:438488/150‑1 (MQ=255)
                                                                                                                         cTAAATAAGTTTGCGTGGGGCCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATAATATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTcc                         <  2:1117819/150‑1 (MQ=255)
                                                                                                                           aaaTAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAATTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCcgc                       >  2:234718/1‑150 (MQ=255)
                                                                                                                           aaaTAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCcgc                       >  2:1797444/1‑150 (MQ=255)
                                                                                                                            aataatTTGGGGTGTGGTTTGTTTATTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCcgcg                      <  2:816024/150‑1 (MQ=255)
                                                                                                                             aTAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCcgcgc                     <  2:7041/150‑1 (MQ=255)
                                                                                                                                 gTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGt                 >  2:837840/1‑150 (MQ=255)
                                                                                                                                       gtgtGGTCAGTAACCTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTAAATGATATTATGCTCAAGTTTAGCGATGAGCGTTATCTGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATTTTTCCGCGCGGGTTTTGtt           >  1:279355/1‑150 (MQ=255)
                                                                                                                                       gggggggcAGGTTACGGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGtt           <  2:352260/143‑1 (MQ=255)
                                                                                                                                        tgtgGTCAGTTTACTTACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTta          <  1:1245168/150‑1 (MQ=255)
                                                                                                                                        tgtgGTCAGTAACCTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTta          >  1:448328/1‑150 (MQ=255)
                                                                                                                                           ggTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTAAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCCCGCGGGTTTTGTtaata       >  2:848565/1‑150 (MQ=255)
                                                                                                                                           ggTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCGGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTtaata       <  2:333158/150‑1 (MQ=255)
                                                                                                                                              cAGGTTACTGACCACTCGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTTCTCCAGTTCCTCTGCGGTGAGGCCCTGTTATATTTCCGCGCGGGTTTTATTAATAACa    >  2:301904/1‑150 (MQ=255)
                                                                                                                                              cAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACa    >  1:985665/1‑150 (MQ=255)
                                                                                                                                              cAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACa    >  2:1466457/1‑150 (MQ=255)
                                                                                                                                              cAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACa    >  2:1512700/1‑150 (MQ=255)
                                                                                                                                              cAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACa    >  1:1229868/1‑150 (MQ=255)
                                                                                                                                               aGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAg   >  1:1077255/1‑150 (MQ=255)
                                                                                                                                                ggTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGc  >  1:1274390/1‑150 (MQ=255)
                                                                                                                                                ggTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGc  >  2:1867593/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
AAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGC  >  NC_000913/1720219‑1720512

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑