Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,849,711 A→C E292A (GAA→GCA)  sppA → protease IV (signal peptide peptidase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,849,7110AC99.5% 690.7 / NA 194E292A (GAA→GCA) sppAprotease IV (signal peptide peptidase)
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (93/100);  total (93/101)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 198 total aligned reads displayed.

GGGTTGCTTGAGGGTTTAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGAAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTATTACGATTACGCATTGAAAACGCCGGCAGATA  >  NC_000913/1849614‑1849806
                                                                                                 |                                                                                               
 ggTTGCTTGAGGGTTTAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGcaaa                                                                                              >  2:3296980/1‑100 (MQ=255)
 ggTTGCTTGAGGGTTTAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGcaaa                                                                                              >  1:1232711/1‑100 (MQ=255)
   ttGCTTGAGGGTTTAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAg                                                                                            <  2:1252201/100‑1 (MQ=255)
        tGAGGGTTTAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTg                                                                                       <  1:300273/100‑1 (MQ=255)
        tGAGGGTTTAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTg                                                                                       <  2:3018238/100‑1 (MQ=255)
         gAGGGTTTAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGa                                                                                      >  2:583960/1‑100 (MQ=255)
          aGGGTTTAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGAc                                                                                     >  1:2544979/1‑100 (MQ=255)
           gggTTTAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCTATGCACTGGCATCGAGTGCGGAAATCGCAAAAGAAGTGAcc                                                                                    >  1:4725233/1‑100 (MQ=255)
           gggTTTAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGAcc                                                                                    <  2:87105/100‑1 (MQ=255)
            ggTTTAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCa                                                                                   >  1:2787613/1‑100 (MQ=255)
            ggTTTAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCa                                                                                   >  2:3913682/1‑100 (MQ=255)
             gTTTAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCaa                                                                                  >  2:3373719/1‑100 (MQ=255)
             gTTTAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCaa                                                                                  >  2:4469139/1‑100 (MQ=255)
             gTTTAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCaa                                                                                  >  2:3642496/1‑100 (MQ=255)
             gTTTAACCAAAACCGGTGGCGATACCGAGAAATATGCACTGGAAAACAAGCGGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCaa                                                                                  <  2:1437463/100‑1 (MQ=255)
               ttAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAg                                                                                >  2:2298691/1‑100 (MQ=255)
               ttAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAg                                                                                <  2:2452250/100‑1 (MQ=255)
                tAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGa                                                                               <  2:4418588/100‑1 (MQ=255)
                   ccAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTTGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAAtt                                                                            >  2:571696/1‑100 (MQ=255)
                   ccAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAAtt                                                                            >  1:854386/1‑100 (MQ=255)
                    cAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTc                                                                           >  2:224052/1‑100 (MQ=255)
                     aaaaCCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCg                                                                          <  1:426101/100‑1 (MQ=255)
                     aaaaCCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCg                                                                          <  2:3767791/100‑1 (MQ=255)
                       aaCCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGc                                                                        <  1:4828355/100‑1 (MQ=255)
                           ggTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTgag                                                                    >  1:2901308/1‑98 (MQ=255)
                           ggTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGa                                                                    >  1:2549939/1‑100 (MQ=255)
                           ggTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGa                                                                    <  2:4792913/100‑1 (MQ=255)
                           ggTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGa                                                                    <  2:1735656/100‑1 (MQ=255)
                           ggTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGa                                                                    >  1:2210371/1‑100 (MQ=255)
                           ggTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGa                                                                    >  1:5482230/1‑100 (MQ=255)
                            gatctgATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAg                                                                   <  1:1913448/95‑1 (MQ=255)
                            gTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAg                                                                   <  1:1027763/100‑1 (MQ=255)
                            gTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAg                                                                   >  2:3864176/1‑100 (MQ=255)
                            gTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAg                                                                   <  1:1277475/100‑1 (MQ=255)
                            gTGGCGATACCGCGAAAGATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAg                                                                   <  1:1197793/100‑1 (MQ=255)
                              ggCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTa                                                                 >  1:4777992/1‑100 (MQ=255)
                              ggCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTa                                                                 <  1:1742330/100‑1 (MQ=255)
                                cGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAg                                                               >  2:2859256/1‑100 (MQ=255)
                                cGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAg                                                               >  1:406217/1‑100 (MQ=255)
                                 gATACCGCGAGATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGa                                                              >  2:5059117/1‑100 (MQ=255)
                                 gATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAgagatc                                                              >  2:1892377/1‑95 (MQ=255)
                                  aTACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGAc                                                             <  1:525307/100‑1 (MQ=255)
                                   tACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACt                                                            <  2:3003211/100‑1 (MQ=255)
                                     ccGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTAGGCTGGCGTAAGACTGa                                                          >  2:1078484/1‑100 (MQ=255)
                                       gcgAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATa                                                        >  1:2831002/1‑100 (MQ=255)
                                           aaTATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATaaaaa                                                    <  2:3205793/100‑1 (MQ=255)
                                           aaTATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATaaaaa                                                    <  1:5212662/100‑1 (MQ=255)
                                             tatGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAAtt                                                  >  1:3881190/1‑100 (MQ=255)
                                               tGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTAt                                                >  1:1136040/1‑100 (MQ=255)
                                               tGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTAt                                                <  2:4368262/100‑1 (MQ=255)
                                                 cACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATcn                                              <  1:4418669/100‑2 (MQ=255)
                                                  aCTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATcgc                                             >  1:1530444/1‑100 (MQ=255)
                                                  aCTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATcgc                                             <  1:1922817/100‑1 (MQ=255)
                                                  aCTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATcgc                                             <  2:1738115/100‑1 (MQ=255)
                                                   cTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATcgcg                                            <  2:2782208/100‑1 (MQ=255)
                                                   cTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATcgcg                                            <  1:2322283/100‑1 (MQ=255)
                                                   cTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATcgcg                                            <  1:4609707/100‑1 (MQ=255)
                                                   cTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAACAATTCGGCTGGAGTAAGACTGATAAAAATTATcgcg                                            >  1:2526496/1‑100 (MQ=255)
                                                   cTGGAAAACAAGCTGGTCGATGCACTGGCACCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATAATcgcg                                            >  2:889197/1‑100 (MQ=255)
                                                     ggAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGcc                                          >  2:240619/1‑100 (MQ=255)
                                                     ggAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGcc                                          <  2:843187/100‑1 (MQ=255)
                                                       aaaaCAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCAt                                        >  1:4551954/1‑100 (MQ=255)
                                                        aaaCAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATc                                       >  1:966707/1‑100 (MQ=255)
                                                           cAAGCTGGTCGATGCACTGGCATCGAGTGTGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGt                                    >  2:5015120/1‑100 (MQ=255)
                                                           cAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGt                                    >  1:3704011/1‑100 (MQ=255)
                                                           cAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACGGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGt                                    >  2:705009/1‑100 (MQ=255)
                                                               cTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGtagat                                >  2:4278439/1‑96 (MQ=255)
                                                                tGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGttatta                               >  1:2447267/1‑100 (MQ=255)
                                                                tGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGttatta                               >  1:4119751/1‑100 (MQ=255)
                                                                tGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGttatta                               >  1:3621489/1‑100 (MQ=255)
                                                                tGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGttatta                               <  2:3649158/100‑1 (MQ=255)
                                                                 ggTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTattac                              <  1:1970069/100‑1 (MQ=255)
                                                                  gTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTattacg                             <  2:4712020/100‑1 (MQ=255)
                                                                    cGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGGTAAAAATTATCGCGCCATCAGTTattacgat                           <  1:5355704/100‑1 (MQ=255)
                                                                    cGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTattacgat                           >  1:4548007/1‑100 (MQ=255)
                                                                     gATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTattacgatt                          <  1:3548162/100‑1 (MQ=255)
                                                                     gATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTattacgatt                          >  2:2067681/1‑100 (MQ=255)
                                                                       tGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTattacgattac                        >  1:744031/1‑100 (MQ=255)
                                                                        gCACTGGCATCGAGTGTGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTattacgattacg                       <  1:4395301/100‑1 (MQ=255)
                                                                         ccagtgCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTATTACGATTACGc                      >  2:1644136/6‑100 (MQ=255)
                                                                         cACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCATTTATTACGATTACGc                      <  2:4171378/100‑1 (MQ=255)
                                                                         cACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTATTACGATTACGc                      <  2:4193311/100‑1 (MQ=255)
                                                                         cACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTATTACGATTACGc                      >  2:465640/1‑100 (MQ=255)
                                                                             ggCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTATTACGATTACGCATTg                  <  1:590938/100‑1 (MQ=255)
                                                                             ggCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTATTACGATTACGCATTg                  <  1:5168102/100‑1 (MQ=255)
                                                                             ggCATCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTATTACGATTACGCATTg                  >  2:2557922/1‑100 (MQ=255)
                                                                                aTCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGTAGTAAGACTGATAAAAATTATCGCGCCATCAGTTATTACGATTACGCATTGaaa               >  2:451738/1‑100 (MQ=255)
                                                                                aTCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTATTACGATTACGCATTGaaa               >  1:4321598/1‑100 (MQ=255)
                                                                                aTCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTATTACGATTACGCATTGaaa               <  1:5052813/100‑1 (MQ=255)
                                                                                 tCGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTATTACGATTACGCATTGaaaa              <  1:3212266/100‑1 (MQ=255)
                                                                                  cGAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTATTACGATTACGCATTGAAAAc             <  2:293347/100‑1 (MQ=255)
                                                                                   gAGTGCGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTATTACGATTACGCATTGAAAACg            <  2:4479744/100‑1 (MQ=255)
                                                                                        cGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTATTACGATTACGCATTGAAAACGCCGGc       <  1:3455826/100‑1 (MQ=255)
                                                                                        cGGAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTATTACGATTACGCATTGAAAACGCCGGc       <  2:2894468/100‑1 (MQ=255)
                                                                                         ggAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTATTACGATTACGCATTGAAAACGCCgaga      >  1:3310905/1‑97 (MQ=255)
                                                                                          naaaTCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTATTACGATTACGCATTGAAAACGCCGGCAg     >  2:2681786/2‑100 (MQ=255)
                                                                                          gAAATCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTATTACGATTACGCATTGAAAACGCCGGCAg     >  1:3443895/1‑100 (MQ=255)
                                                                                            aaTCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTATTACGATTACGCATTGAAAACGCCGGCAGAt   <  1:2724794/100‑1 (MQ=255)
                                                                                            aaTCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTATTACGATTACGCATTGAAAACGCCGGCAGAt   <  1:3088442/100‑1 (MQ=255)
                                                                                             aTCGCAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTATTACGATTACGCATTGAAAACGCCGGCAGATa  >  1:85802/1‑100 (MQ=255)
                                                                                                 |                                                                                               
GGGTTGCTTGAGGGTTTAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGAAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTATTACGATTACGCATTGAAAACGCCGGCAGATA  >  NC_000913/1849614‑1849806

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