Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,850,719:1 +T 16.7% intergenic (+26/‑141) 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,7191.T16.7% 302.6 / 70.6 102intergenic (+26/‑141)sppA/ansAprotease IV (signal peptide peptidase)/cytoplasmic L‑asparaginase 1
Reads supporting (aligned to +/- strand):  ref base . (50/35);  new base T (6/11);  total (56/46)
Fisher's exact test for biased strand distribution p-value = 1.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.69e-01

Only 100 of 102 total aligned reads displayed.

GCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTC  >  NC_000913/1850627‑1850816
                                                                                             |                                                                                                 
gCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑GATCTGTCGGcc                                                                                            <  2:1843417/100‑1 (MQ=255)
 cTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCca                                                                                           <  2:2086789/100‑1 (MQ=255)
 cTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑GATCTGTCGGCca                                                                                           <  2:1884123/100‑1 (MQ=255)
  tGGCCGCGTTTAATGACCCACAAAACCGTTATGTGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑GATCTGTCGGCcac                                                                                          <  1:1097038/100‑1 (MQ=255)
  tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCcac                                                                                          <  2:2775758/100‑1 (MQ=255)
  tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCcac                                                                                          <  2:1516262/100‑1 (MQ=255)
  tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑GATCTGTCGGacac                                                                                          <  1:2230481/100‑1 (MQ=255)
  tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑GATCTGTCGGCcac                                                                                          <  2:2410111/100‑1 (MQ=255)
  tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑GATCTGTCGGCcac                                                                                          <  2:1222372/100‑1 (MQ=255)
  tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑GATCTGTCGGCcac                                                                                          <  1:217873/100‑1 (MQ=255)
  tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCGTGTACTGAGTGGCCGACAGATC‑GTCGGCcac                                                                                          <  1:2722936/100‑1 (MQ=255)
    gCCGCGTTTAATGACCCCCAAGGCCGGTATGCGTTTTGCCTGAACTTCGCCAACATGCGTTAAGTCTTGTACTGGGTGGCCGAC‑GATCTGTCGGcc                                                                                            <  1:590881/96‑1 (MQ=255)
    gCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACat                                                                                        <  2:2255697/100‑1 (MQ=255)
    gCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACat                                                                                        >  2:2831857/1‑100 (MQ=255)
    gCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTGGCCTGACCTGAGCCAATATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑GATCTGTCGGGcactc                                                                                        <  1:1616025/100‑3 (MQ=255)
    gCCGCGTTTAATGACCCACAAAACCGTTAGGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑GATCTGTCGGCcactc                                                                                        <  1:329221/100‑3 (MQ=255)
      cgcgTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACatta                                                                                      >  2:1205863/1‑100 (MQ=255)
        cgTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACattatt                                                                                    <  2:2415477/100‑1 (MQ=255)
        cgTTTAATGACCCACAAAACCGATATGCGTTTTGACTGACATGCGCCAACATGCGTTAAATCGTGTACTCAGTGGCCGACAGATC‑GTCGGGCACatta                                                                                      >  1:1370237/1‑98 (MQ=255)
          tttcacgccccACAAAACCGTTATGCGTTTCGCCTGACCTACGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACATATC‑GTCGGCCACATTAtttt                                                                                  <  1:1095405/92‑1 (MQ=255)
            tAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTAtttttt                                                                                <  1:624428/100‑1 (MQ=255)
             aaTGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTa                                                                               <  2:2916232/100‑1 (MQ=255)
             aaTGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTa                                                                               <  2:2242628/100‑1 (MQ=255)
              acgACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTAc                                                                              <  1:221301/98‑1 (MQ=255)
              aTGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAATGGCCGACAGATC‑GTCGGCCACATTATTTTTTAc                                                                              >  1:1535774/1‑100 (MQ=255)
                  cccACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTcg                                                                          >  1:1165154/1‑100 (MQ=255)
                  cccACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTcg                                                                          >  1:2707572/1‑100 (MQ=255)
                  cccACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTcg                                                                          >  1:2633635/1‑100 (MQ=255)
                    cacaAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTcgac                                                                        >  2:1323264/1‑100 (MQ=255)
                    cacaAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTcgac                                                                        <  2:2152774/100‑1 (MQ=255)
                    cacaAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTcgac                                                                        <  2:2222780/100‑1 (MQ=255)
                    cacaAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTcgac                                                                        >  2:2691182/1‑100 (MQ=255)
                    cacaAAACCGTTATGCGTTGTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTcgac                                                                        <  1:891143/100‑1 (MQ=255)
                     acaAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTcgacg                                                                       <  1:2569297/100‑1 (MQ=255)
                       aaaaCCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGaa                                                                     >  2:1903480/1‑100 (MQ=255)
                        aaaCCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAAt                                                                    <  1:1121233/100‑1 (MQ=255)
                           ccGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCct                                                                 >  2:380368/1‑100 (MQ=255)
                             nttATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCctct                                                               >  1:1942224/2‑100 (MQ=255)
                             gTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCctct                                                               <  2:207938/100‑1 (MQ=255)
                             gTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCctct                                                               <  1:1501066/100‑1 (MQ=255)
                               tATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTc                                                             <  1:2715888/100‑1 (MQ=255)
                                 tGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTccc                                                           >  2:2049218/1‑100 (MQ=255)
                                   cGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCCGCCACATTATTTTTTACGTCGACGAATCCTCTTccc                                                           <  1:1203380/98‑1 (MQ=255)
                                    gTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCt                                                        >  1:896097/1‑100 (MQ=255)
                                        tGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGttt                                                    >  1:261335/1‑100 (MQ=255)
                                         gCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTc                                                   >  2:199296/1‑100 (MQ=255)
                                         gCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTc                                                   >  1:134674/1‑100 (MQ=255)
                                            tGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGcc                                                <  1:1579423/100‑1 (MQ=255)
                                             gACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGccc                                               <  2:252175/100‑1 (MQ=255)
                                               ccTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCa                                             >  2:1347774/1‑100 (MQ=255)
                                               ccTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCa                                             >  2:2891732/1‑100 (MQ=255)
                                                cTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCa                                             <  1:1324951/99‑1 (MQ=255)
                                                 tGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCata                                           >  2:1009772/1‑100 (MQ=255)
                                                  gcgcCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCatat                                          >  2:2324583/1‑100 (MQ=255)
                                                  gcgcCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCatat                                          >  2:1671906/1‑100 (MQ=255)
                                                    gtcAACATGCGTTCAGTCTTGTACTGAGTGGCCGACAGATC‑GACGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATAttt                                        >  2:1974170/3‑100 (MQ=255)
                                                    gcCAACATGCGGTAAGTCTTGCACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCGTCTTCCCGCTGTTTCGCTCCATAtt                                         >  1:2178048/1‑99 (MQ=255)
                                                     ccAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTc                                       >  1:573885/1‑100 (MQ=255)
                                                     ccAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTc                                       <  2:1141155/100‑1 (MQ=255)
                                                       aaCATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCt                                     >  2:465096/1‑100 (MQ=255)
                                                       aaCATGCGTGCAGTCTTGTACTGCGTGGTCGACAGATG‑GTCGGCCACATTATTTTTTACTTCGTCGAATCCTCTTCTAGCTGTTTCGCTCCGTATTTc                                       >  1:1208592/1‑98 (MQ=255)
                                                         cATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTa                                   <  2:1872168/100‑1 (MQ=255)
                                                           tGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCTGCCACATTATTTTTGACGTCGACGAATCCTCTTCCCGCTGGTTCGCCCCATATTTCCTTATc                                 <  1:2712085/100‑1 (MQ=255)
                                                            gCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATcc                                >  1:1504284/1‑100 (MQ=255)
                                                              gTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGc                              >  2:575996/1‑100 (MQ=255)
                                                              gTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGc                              >  2:603368/1‑100 (MQ=255)
                                                                tAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCAt                            <  1:1395540/100‑1 (MQ=255)
                                                                 aaGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATa                           >  2:1302901/1‑100 (MQ=255)
                                                                   gTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGt                         >  1:2182717/1‑100 (MQ=255)
                                                                   gTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCGCTTCCCGCTGTTTCGCCCCATATTTCATTATCCGCATAGt                         >  1:699767/1‑100 (MQ=255)
                                                                     cTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCAATATTTCCTTATCCGCATAGTAt                       >  2:1616568/1‑100 (MQ=255)
                                                                      ttGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATc                      >  2:2693063/1‑100 (MQ=255)
                                                                       tGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATcn                     <  2:861106/100‑2 (MQ=255)
                                                                       tGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCa                     <  2:502534/100‑1 (MQ=255)
                                                                       tGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCa                     <  2:573888/100‑1 (MQ=255)
                                                                        gTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAg                    >  2:2388427/1‑100 (MQ=255)
                                                                         tACTGAGTGGCCGACAGCTC‑GTCGGCCACATTTTTTTTTACGTCTACGAATGCTCTTCCCGCTGCGTCCCTCCATATTTCCTTATCCGCATCGTATCAgg                   >  1:405312/1‑100 (MQ=255)
                                                                         aataatgTGGCCGAC‑GATCTGTTGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAgg                   >  1:2537014/7‑100 (MQ=255)
                                                                          ataatgTGGCCGAC‑GATCTGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAgg                   >  2:109860/6‑99 (MQ=38)
                                                                          aCTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGt                  <  1:2497549/100‑1 (MQ=255)
                                                                            aatgTGGCCGAC‑GATCTGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTgc                >  2:1490694/4‑100 (MQ=39)
                                                                             gAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTgcg               >  2:2325233/1‑100 (MQ=255)
                                                                             atgTGGCCGAC‑GATCTGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTgcg               >  2:828166/3‑100 (MQ=255)
                                                                              aGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTgcgc              <  2:2340237/100‑1 (MQ=255)
                                                                               gTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCt             >  2:640594/1‑100 (MQ=255)
                                                                               gTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATACTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCt             >  1:1184674/1‑100 (MQ=255)
                                                                               gTGGCCGAC‑GATCTGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCt             >  2:37076/1‑100 (MQ=255)
                                                                                tGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTtcgc              >  1:336207/1‑98 (MQ=255)
                                                                                tGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTc            >  2:1762058/1‑100 (MQ=255)
                                                                                 ggCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTcc           >  1:2852787/1‑100 (MQ=255)
                                                                                 ggCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTcc           >  2:2400071/1‑100 (MQ=255)
                                                                                   ccGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTcccc         >  2:94282/1‑100 (MQ=255)
                                                                                   ccGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTcccc         >  2:2137076/1‑100 (MQ=255)
                                                                                    cGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACTAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATGCGCATAGTATCAGGTGCGCTccccc        >  1:523492/1‑100 (MQ=255)
                                                                                        aGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTTCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGcc    <  1:2049207/100‑1 (MQ=255)
                                                                                        aGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGcc    <  2:86999/100‑1 (MQ=255)
                                                                                        aGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGcc    >  1:1996575/1‑100 (MQ=255)
                                                                                         gATT‑GTCGGCCACATTATTTCTTGCGACGAAGCATCCTCTTCCTGCTGCGTCGCCGCATCTTCCGTTAGCCGCGTAGTATCGGGGGCGCTCCGCCTGcc    <  1:2007368/99‑1 (MQ=255)
                                                                                         gATC‑GTCGGCCACATTAGTTTTTACGGGGACGAATCCTCTTCCCGCTGTCTCGCCCCACATTTCCGTATCCCCATAGTATGAGGTGGGCCCGCCCTccc    >  1:970270/1‑99 (MQ=255)
                                                                                          aTC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCcctccctc  >  1:1453501/1‑100 (MQ=255)
                                                                                             |                                                                                                 
GCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTC  >  NC_000913/1850627‑1850816

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