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.

AAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCC  >  NC_000913/1850617‑1850814
                                                                                                  |                                                                                                    
aaagtgaCAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAg                                                                                                     <  1:1747614/100‑1 (MQ=255)
  agtgaCGAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGAt                                                                                                   <  2:1600802/100‑1 (MQ=255)
       cAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGAtcgtcg                                                                                             >  2:2897571/1‑100 (MQ=255)
       cAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGAtcgtc                                                                                              >  2:3179152/1‑99 (MQ=255)
        aaGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCgg                                                                                            <  1:2087593/100‑1 (MQ=255)
         aGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGc                                                                                           <  1:2482297/100‑1 (MQ=255)
         aGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGc                                                                                           <  1:1421453/100‑1 (MQ=255)
          gCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGcc                                                                                          >  1:63099/1‑100 (MQ=255)
          gCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACGATCTGTCGGcc                                                                                          <  1:1767507/100‑1 (MQ=255)
           cTGGCTGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACGATCTGTCGGCca                                                                                         <  2:3182285/100‑1 (MQ=255)
           cTGGCCGCGTTTAATTACCCACAAAACCGTTATGCGTTTTGCTTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGAtcttagtcac                                                                                          >  2:2901288/1‑92 (MQ=255)
           cTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACGATCTGTCGGCca                                                                                         <  2:2389435/100‑1 (MQ=255)
            tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGcccc                                                                                        >  2:2316799/1‑98 (MQ=255)
            tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCcac                                                                                        <  1:3294805/100‑1 (MQ=255)
            tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCcac                                                                                        >  1:3578101/1‑100 (MQ=255)
            tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACGATCTGTCGGCcac                                                                                        <  1:1508675/100‑1 (MQ=255)
            tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACGATCTGTCGGCcac                                                                                        <  1:127103/100‑1 (MQ=255)
            tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACGATCTGTCGGCcac                                                                                        <  1:2059382/100‑1 (MQ=255)
            tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACGATCTGTCGGCcac                                                                                        <  1:3424761/100‑1 (MQ=255)
            tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACGATCTGTCGGCcac                                                                                        <  1:1059986/100‑1 (MQ=255)
            tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACGATCTGTCGGCcac                                                                                        <  1:2731218/100‑1 (MQ=255)
            tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTATTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCcac                                                                                        >  2:635531/1‑100 (MQ=255)
              gCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTGAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACat                                                                                      >  2:2738906/1‑100 (MQ=255)
              gCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACGATCTGTCGGCcactc                                                                                      <  2:1733579/100‑3 (MQ=255)
              gCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGACTTGTACTGAGTGGCCGACAGATCGTCGGCCACat                                                                                      >  2:1876445/1‑100 (MQ=255)
                ctcgTTTAATGACCCACAAAACCGGTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGGACTGAGTGGCCGACGATCTGTCGGCcactc                                                                                      <  2:3530416/96‑3 (MQ=38)
                 gcgTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACattat                                                                                   >  1:1532230/1‑100 (MQ=255)
                 gcgTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACGATCTGTCGGCcactcagt                                                                                   <  1:890587/100‑6 (MQ=38)
                    tttAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTAtttt                                                                                <  1:1109043/100‑1 (MQ=255)
                    tttAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTAtttt                                                                                <  1:1108957/100‑1 (MQ=255)
                     ttAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTAttttt                                                                               <  2:2692796/100‑1 (MQ=255)
                      tagtGACCCACAGAACCGTTATGCGTTTTGCCTGTCCTGCGCCAACGTGCGTTAAGGCTTGTACGGGGTGGCCGACGATCTGTCGGCcactc                                                                                      <  2:1339724/89‑3 (MQ=21)
                      tAATGACCCACAATAGCGTTATGTGTTTTGCCTGTCCTGCGCGAACATGCGTTAAGTCTTGTACTGTGTGGCCGTCAGATCGTCGGCCACATTAtttttt                                                                              <  2:1669663/100‑1 (MQ=255)
                      tAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTAtttttt                                                                              <  2:3428415/100‑1 (MQ=255)
                      tAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTAtttttt                                                                              <  2:586215/100‑1 (MQ=255)
                      tAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTAtttttt                                                                              >  2:2547133/1‑100 (MQ=255)
                       aaTGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTa                                                                             >  1:3782783/1‑100 (MQ=255)
                        aTGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTAc                                                                            <  1:1045575/100‑1 (MQ=255)
                        aTGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTAc                                                                            >  2:2820714/1‑100 (MQ=255)
                         tGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTCCTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACg                                                                           >  2:1796182/1‑100 (MQ=255)
                           aCCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTc                                                                         <  2:378639/100‑1 (MQ=255)
                           aCCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTc                                                                         <  2:642347/100‑1 (MQ=255)
                            cccACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAAAGATCGTCGGCCACATTATTTTTTACGTcg                                                                        >  1:1657397/1‑100 (MQ=255)
                               acaAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTAGGCGACATTATTTTTTACGTcgac                                                                      <  2:382190/99‑1 (MQ=255)
                                caAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTcgacga                                                                    <  1:3289059/100‑1 (MQ=255)
                                caAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAACTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTcgacga                                                                    >  2:103182/1‑100 (MQ=255)
                                    aCCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATcc                                                                >  1:3527638/1‑100 (MQ=255)
                                     ccGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCct                                                               <  2:2665908/100‑1 (MQ=255)
                                      cGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCctc                                                              <  1:960245/100‑1 (MQ=255)
                                       gTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCctct                                                             >  2:128034/1‑100 (MQ=255)
                                          aTGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTcc                                                          >  1:1602937/1‑100 (MQ=255)
                                           tGCGTTTTGCCTGATCTGCGCCAACAGGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTcc                                                          <  2:1635045/99‑1 (MQ=255)
                                             cGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGCGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGc                                                       <  2:2440905/100‑1 (MQ=255)
                                              gTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCt                                                      >  2:1056461/1‑100 (MQ=255)
                                               ttttGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCtn                                                     <  1:3220116/100‑2 (MQ=255)
                                                 ttGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGtt                                                   <  2:1823003/100‑1 (MQ=255)
                                                  tGCCTGACCTGCGACAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGttt                                                  <  2:1431854/100‑1 (MQ=255)
                                                   gCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTc                                                 >  1:177262/1‑100 (MQ=255)
                                                     cTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGc                                               >  1:1444883/1‑100 (MQ=255)
                                                     cTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGc                                               >  2:3253257/1‑100 (MQ=255)
                                                      tGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGcc                                              >  1:1109724/1‑100 (MQ=255)
                                                       gACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCTACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGcc                                              <  2:1558417/99‑1 (MQ=255)
                                                       gACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGccc                                             <  1:2969568/100‑1 (MQ=255)
                                                       gACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGccc                                             >  1:3457779/1‑100 (MQ=255)
                                                                 aaCATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCt                                   <  2:777265/100‑1 (MQ=255)
                                                                     tGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATc                               <  2:1602927/100‑1 (MQ=255)
                                                                     tGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATc                               >  2:535709/1‑100 (MQ=255)
                                                                       cGTTAAGTCTTGTACTGAGTGGCCGACAGGTCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCg                             <  1:1796192/100‑1 (MQ=255)
                                                                       cGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCg                             >  1:1972116/1‑100 (MQ=255)
                                                                       cGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCg                             >  2:2695686/1‑100 (MQ=255)
                                                                         ttAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCa                           <  1:40377/100‑1 (MQ=255)
                                                                          ttccgatcttACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCAt                          <  2:2783169/91‑1 (MQ=255)
                                                                          tAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCAt                          <  1:2896574/100‑1 (MQ=255)
                                                                          tAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCAt                          >  1:1250242/1‑100 (MQ=255)
                                                                           aaGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATa                         <  1:1529114/100‑1 (MQ=255)
                                                                             gTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCACTGTTTCGCCCCATATTTCCTTATCCGCATAGt                       >  1:1108049/1‑100 (MQ=255)
                                                                               cTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTAt                     >  2:3438671/1‑100 (MQ=255)
                                                                               cTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTAt                     >  1:2850187/1‑100 (MQ=255)
                                                                               cTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTAt                     >  1:2348928/1‑100 (MQ=255)
                                                                                aaaaataatgTGGCCGACGATCTGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATc                    >  1:3389609/10‑100 (MQ=37)
                                                                                  gTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAg                  >  1:701934/1‑100 (MQ=255)
                                                                                  gTACTGAGTGGCCGACAGATCGTCGGCAACATCATTATTTACGTCAACAAATCCTATTCCATCCGCTTCGCCGCCTATTTCCTTAGCCCCACAg                        >  2:2319059/1‑94 (MQ=255)
                                                                                   tACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAgatcggaa                 >  1:2783189/1‑93 (MQ=255)
                                                                                     cTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTg               >  2:3367361/1‑100 (MQ=255)
                                                                                      tGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTgc              <  2:488601/100‑1 (MQ=255)
                                                                                      aatgTGGCCGACGATCTGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTgc              >  1:3120427/4‑100 (MQ=39)
                                                                                       gAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTgcg             >  2:3296771/1‑100 (MQ=255)
                                                                                       gAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTgcg             >  1:3738791/1‑100 (MQ=255)
                                                                                        tgTGGCCGACGATCTGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTgcgc            >  1:2378450/2‑100 (MQ=255)
                                                                                        aGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTgcgc            <  1:3253283/100‑1 (MQ=255)
                                                                                         gTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGATTCCTTTTCCCGCCGTTTCGCCCCATATGTCCTTATCCGCATAGTATCTGGTACGCt           >  1:606229/1‑100 (MQ=255)
                                                                                         gTGGCCGACGATCTGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCt           >  1:2753819/1‑100 (MQ=255)
                                                                                           ggCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTcc         >  1:510401/1‑100 (MQ=255)
                                                                                           ggCCGACGATCTGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTCCGCTcc         >  2:1200297/1‑100 (MQ=255)
                                                                                            gCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTccc        >  1:2899459/1‑100 (MQ=255)
                                                                                              cGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTccccc      <  1:1062761/100‑1 (MQ=255)
                                                                                                 cAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGc   <  2:1532220/100‑1 (MQ=255)
                                                                                                  aGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGcc  <  2:3457753/100‑1 (MQ=255)
                                                                                                  aGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGcc  <  2:3612376/100‑1 (MQ=255)
                                                                                                  aGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGcc  <  2:1728653/100‑1 (MQ=255)
                                                                                                  |                                                                                                    
AAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCC  >  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: