Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,850,715 Δ1 bp 17.3% 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.17.3% 279.1 / 55.7 98intergenic (+22/‑145)sppA/ansAprotease IV (signal peptide peptidase)/cytoplasmic L‑asparaginase 1
Reads supporting (aligned to +/- strand):  ref base A (49/32);  new base . (6/11);  total (55/43)
Fisher's exact test for biased strand distribution p-value = 6.58e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 101 total aligned reads displayed.

AAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCC  >  NC_000913/1850617‑1850814
                                                                                                  |                                                                                                    
aaagtgaCAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAg                                                                                                     >  1:207937/1‑100 (MQ=255)
 aagtgcCAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTCCTGCCTTACCTGCGCCAACACGCGTTACGGCTTGTACTGAGCGGCCGACAg                                                                                                     >  1:567996/1‑99 (MQ=255)
          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)
                cgcgTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑GATCTGTCGGCcactcag                                                                                    <  2:1182224/100‑5 (MQ=38)
                  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:2633635/1‑100 (MQ=255)
                            cccACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTcg                                                                        >  1:2707572/1‑100 (MQ=255)
                            cccACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTcg                                                                        >  1:1165154/1‑100 (MQ=255)
                              cacaAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTcgac                                                                      <  2:2152774/100‑1 (MQ=255)
                              cacaAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTcgac                                                                      <  2:2222780/100‑1 (MQ=255)
                              cacaAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTcgac                                                                      >  2:1323264/1‑100 (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                                                             <  1:1501066/100‑1 (MQ=255)
                                       gTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCctct                                                             <  2:207938/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:2891732/1‑100 (MQ=255)
                                                         ccTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCa                                           >  2:1347774/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                                     <  2:1141155/100‑1 (MQ=255)
                                                               ccAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTc                                     >  1:573885/1‑100 (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:603368/1‑100 (MQ=255)
                                                                        gTTAAGTCTTGTACTGAGTGGCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGc                            >  2:575996/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         >  2:2400071/1‑100 (MQ=255)
                                                                                           ggCCGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTcc         >  1:2852787/1‑100 (MQ=255)
                                                                                             ccGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTcccc       >  2:2137076/1‑100 (MQ=255)
                                                                                             ccGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTcccc       >  2:94282/1‑100 (MQ=255)
                                                                                              cGACAGATC‑GTCGGCCACATTATTTTTTACGTCGACTAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATGCGCATAGTATCAGGTGCGCTccccc      >  1:523492/1‑100 (MQ=255)
                                                                                                  aGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTTCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGcc  <  1:2049207/100‑1 (MQ=255)
                                                                                                  aGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGcc  >  1:1996575/1‑100 (MQ=255)
                                                                                                  aGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGcc  <  2:86999/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: ‑