Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,850,715 Δ1 bp 19.6% intergenic (+22/‑145) sppA → / → ansA protease IV (signal peptide peptidase)/cytoplasmic L‑asparaginase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,850,7150A.19.6% 302.3 / 75.1 112intergenic (+22/‑145)sppA/ansAprotease IV (signal peptide peptidase)/cytoplasmic L‑asparaginase 1
Reads supporting (aligned to +/- strand):  ref base A (48/42);  new base . (6/16);  total (54/58)
Fisher's exact test for biased strand distribution p-value = 3.36e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

Only 100 of 115 total aligned reads displayed.

AAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCC  >  NC_000913/1850617‑1850814
                                                                                                  |                                                                                                    
aaagtgaCAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAg                                                                                                     <  1:1747614/100‑1 (MQ=255)
  agtgaCGAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGAt                                                                                                   <  2:1600802/100‑1 (MQ=255)
       cAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGAtc‑gtcg                                                                                             >  2:2897571/1‑100 (MQ=255)
       cAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGAtc‑gtc                                                                                              >  2:3179152/1‑99 (MQ=255)
        aaGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCgg                                                                                            <  1:2087593/100‑1 (MQ=255)
         aGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGc                                                                                           <  1:2482297/100‑1 (MQ=255)
         aGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGc                                                                                           <  1:1421453/100‑1 (MQ=255)
          gCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGcc                                                                                          >  1:63099/1‑100 (MQ=255)
          gCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑GATCTGTCGGcc                                                                                          <  1:1767507/100‑1 (MQ=255)
           cTGGCTGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑GATCTGTCGGCca                                                                                         <  2:3182285/100‑1 (MQ=255)
           cTGGCCGCGTTTAATTACCCACAAAACCGTTATGCGTTTTGCTTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGAtcttagtcac                                                                                          >  2:2901288/1‑92 (MQ=255)
           cTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑GATCTGTCGGCca                                                                                         <  2:2389435/100‑1 (MQ=255)
            tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGcccc                                                                                        >  2:2316799/1‑98 (MQ=255)
            tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCcac                                                                                        <  1:3294805/100‑1 (MQ=255)
            tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCcac                                                                                        >  1:3578101/1‑100 (MQ=255)
            tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑GATCTGTCGGCcac                                                                                        <  1:1508675/100‑1 (MQ=255)
            tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑GATCTGTCGGCcac                                                                                        <  1:127103/100‑1 (MQ=255)
            tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑GATCTGTCGGCcac                                                                                        <  1:2059382/100‑1 (MQ=255)
            tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑GATCTGTCGGCcac                                                                                        <  1:3424761/100‑1 (MQ=255)
            tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑GATCTGTCGGCcac                                                                                        <  1:1059986/100‑1 (MQ=255)
            tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑GATCTGTCGGCcac                                                                                        <  1:2731218/100‑1 (MQ=255)
            tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTATTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCcac                                                                                        >  2:635531/1‑100 (MQ=255)
              gCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTGAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACat                                                                                      >  2:2738906/1‑100 (MQ=255)
              gCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑GATCTGTCGGCcactc                                                                                      <  2:1733579/100‑3 (MQ=255)
              gCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGACTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACat                                                                                      >  2:1876445/1‑100 (MQ=255)
                ctcgTTTAATGACCCACAAAACCGGTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGGACTGAGTGGCCGAC‑GATCTGTCGGCcactc                                                                                      <  2:3530416/96‑3 (MQ=38)
                 gcgTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACattat                                                                                   >  1:1532230/1‑100 (MQ=255)
                 gcgTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑GATCTGTCGGCcactcagt                                                                                   <  1:890587/100‑6 (MQ=38)
                    tttAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTAtttt                                                                                <  1:1109043/100‑1 (MQ=255)
                    tttAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTAtttt                                                                                <  1:1108957/100‑1 (MQ=255)
                     ttAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTAttttt                                                                               <  2:2692796/100‑1 (MQ=255)
                      tagtGACCCACAGAACCGTTATGCGTTTTGCCTGTCCTGCGCCAACGTGCGTTAAGGCTTGTACGGGGTGGCCGAC‑GATCTGTCGGCcactc                                                                                      <  2:1339724/89‑3 (MQ=21)
                      tAATGACCCACAATAGCGTTATGTGTTTTGCCTGTCCTGCGCGAACATGCGTTAAGTCTTGTACTGTGTGGCCGTCAGATC‑GTCGGCCACATTAtttttt                                                                              <  2:1669663/100‑1 (MQ=255)
                      tAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTAtttttt                                                                              <  2:3428415/100‑1 (MQ=255)
                      tAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTAtttttt                                                                              <  2:586215/100‑1 (MQ=255)
                      tAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTAtttttt                                                                              >  2:2547133/1‑100 (MQ=255)
                       aaTGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTa                                                                             >  1:3782783/1‑100 (MQ=255)
                        aTGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTAc                                                                            <  1:1045575/100‑1 (MQ=255)
                        aTGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTAc                                                                            >  2:2820714/1‑100 (MQ=255)
                         tGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTCCTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACg                                                                           >  2:1796182/1‑100 (MQ=255)
                           aCCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTc                                                                         <  2:378639/100‑1 (MQ=255)
                           aCCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTc                                                                         <  2:642347/100‑1 (MQ=255)
                            cccACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAAAGATC‑GTCGGCCACATTATTTTTTACGTcg                                                                        >  1:1657397/1‑100 (MQ=255)
                               acaAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTAGGCGACATTATTTTTTACGTcgac                                                                      <  2:382190/99‑1 (MQ=255)
                                caAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTcgacga                                                                    <  1:3289059/100‑1 (MQ=255)
                                caAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAACTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTcgacga                                                                    >  2:103182/1‑100 (MQ=255)
                                    aCCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATcc                                                                >  1:3527638/1‑100 (MQ=255)
                                     ccGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCct                                                               <  2:2665908/100‑1 (MQ=255)
                                      cGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCctc                                                              <  1:960245/100‑1 (MQ=255)
                                       gTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCctct                                                             >  2:128034/1‑100 (MQ=255)
                                          aTGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTcc                                                          >  1:1602937/1‑100 (MQ=255)
                                           tGCGTTTTGCCTGATCTGCGCCAACAGGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTcc                                                          <  2:1635045/99‑1 (MQ=255)
                                             cGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGCGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGc                                                       <  2:2440905/100‑1 (MQ=255)
                                              gTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCt                                                      >  2:1056461/1‑100 (MQ=255)
                                               ttttGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCtn                                                     <  1:3220116/100‑2 (MQ=255)
                                                 ttGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGtt                                                   <  2:1823003/100‑1 (MQ=255)
                                                  tGCCTGACCTGCGACAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGttt                                                  <  2:1431854/100‑1 (MQ=255)
                                                   gCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTc                                                 >  1:177262/1‑100 (MQ=255)
                                                     cTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGc                                               >  1:1444883/1‑100 (MQ=255)
                                                     cTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGc                                               >  2:3253257/1‑100 (MQ=255)
                                                      tGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGcc                                              >  1:1109724/1‑100 (MQ=255)
                                                       gACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCTACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGcc                                              <  2:1558417/99‑1 (MQ=255)
                                                       gACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGccc                                             <  1:2969568/100‑1 (MQ=255)
                                                       gACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGccc                                             >  1:3457779/1‑100 (MQ=255)
                                                                 aaCATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCt                                   <  2:777265/100‑1 (MQ=255)
                                                                     tGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATc                               <  2:1602927/100‑1 (MQ=255)
                                                                     tGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATc                               >  2:535709/1‑100 (MQ=255)
                                                                       cGTTAAGTCTTGTACTGAGTGGCCGACAGGTC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCg                             <  1:1796192/100‑1 (MQ=255)
                                                                       cGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCg                             >  1:1972116/1‑100 (MQ=255)
                                                                       cGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCg                             >  2:2695686/1‑100 (MQ=255)
                                                                         ttAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCa                           <  1:40377/100‑1 (MQ=255)
                                                                          ttccgatcttACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCAt                          <  2:2783169/91‑1 (MQ=255)
                                                                          tAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCAt                          <  1:2896574/100‑1 (MQ=255)
                                                                          tAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCAt                          >  1:1250242/1‑100 (MQ=255)
                                                                           aaGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATa                         <  1:1529114/100‑1 (MQ=255)
                                                                             gTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCACTGTTTCGCCCCATATTTCCTTATCCGCATAGt                       >  1:1108049/1‑100 (MQ=255)
                                                                               cTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTAt                     >  2:3438671/1‑100 (MQ=255)
                                                                               cTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTAt                     >  1:2850187/1‑100 (MQ=255)
                                                                               cTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTAt                     >  1:2348928/1‑100 (MQ=255)
                                                                                aaaaataatgTGGCCGAC‑GATCTGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATc                    >  1:3389609/10‑100 (MQ=37)
                                                                                  gTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAg                  >  1:701934/1‑100 (MQ=255)
                                                                                  gTACTGAGTGGCCGACAGATC‑GTCGGCAACATCATTATTTACGTCAACAAATCCTATTCCATCCGCTTCGCCGCCTATTTCCTTAGCCCCACAg                        >  2:2319059/1‑94 (MQ=255)
                                                                                   tACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAgatcggaa                 >  1:2783189/1‑93 (MQ=255)
                                                                                     cTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTg               >  2:3367361/1‑100 (MQ=255)
                                                                                      tGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTgc              <  2:488601/100‑1 (MQ=255)
                                                                                      aatgTGGCCGAC‑GATCTGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTgc              >  1:3120427/4‑100 (MQ=39)
                                                                                       gAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTgcg             >  2:3296771/1‑100 (MQ=255)
                                                                                       gAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTgcg             >  1:3738791/1‑100 (MQ=255)
                                                                                        tgTGGCCGAC‑GATCTGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTgcgc            >  1:2378450/2‑100 (MQ=255)
                                                                                        aGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTgcgc            <  1:3253283/100‑1 (MQ=255)
                                                                                         gTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGATTCCTTTTCCCGCCGTTTCGCCCCATATGTCCTTATCCGCATAGTATCTGGTACGCt           >  1:606229/1‑100 (MQ=255)
                                                                                         gTGGCCGAC‑GATCTGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCt           >  1:2753819/1‑100 (MQ=255)
                                                                                           ggCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTcc         >  1:510401/1‑100 (MQ=255)
                                                                                           ggCCGAC‑GATCTGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTCCGCTcc         >  2:1200297/1‑100 (MQ=255)
                                                                                            gCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTccc        >  1:2899459/1‑100 (MQ=255)
                                                                                              cGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTccccc      <  1:1062761/100‑1 (MQ=255)
                                                                                                 cAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGc   <  2:1532220/100‑1 (MQ=255)
                                                                                                  aGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGcc  <  2:3457753/100‑1 (MQ=255)
                                                                                                  aGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGcc  <  2:3612376/100‑1 (MQ=255)
                                                                                                  aGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGcc  <  2:1728653/100‑1 (MQ=255)
                                                                                                  |                                                                                                    
AAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCC  >  NC_000913/1850617‑1850814

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