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.

AAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCC  >  NC_000913/1850617‑1850814
                                                                                                  |                                                                                                    
aaagtgaCAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAg                                                                                                     >  1:207937/1‑100 (MQ=255)
 aagtgcCAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTCCTGCCTTACCTGCGCCAACACGCGTTACGGCTTGTACTGAGCGGCCGACAg                                                                                                     >  1:567996/1‑99 (MQ=255)
          gCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACGATCTGTCGGcc                                                                                          <  2:1843417/100‑1 (MQ=255)
           cTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCca                                                                                         <  2:2086789/100‑1 (MQ=255)
           cTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACGATCTGTCGGCca                                                                                         <  2:1884123/100‑1 (MQ=255)
            tGGCCGCGTTTAATGACCCACAAAACCGTTATGTGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACGATCTGTCGGCcac                                                                                        <  1:1097038/100‑1 (MQ=255)
            tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCcac                                                                                        <  2:2775758/100‑1 (MQ=255)
            tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCcac                                                                                        <  2:1516262/100‑1 (MQ=255)
            tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACGATCTGTCGGacac                                                                                        <  1:2230481/100‑1 (MQ=255)
            tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACGATCTGTCGGCcac                                                                                        <  2:2410111/100‑1 (MQ=255)
            tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACGATCTGTCGGCcac                                                                                        <  2:1222372/100‑1 (MQ=255)
            tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACGATCTGTCGGCcac                                                                                        <  1:217873/100‑1 (MQ=255)
            tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCGTGTACTGAGTGGCCGACAGATCGTCGGCcac                                                                                        <  1:2722936/100‑1 (MQ=255)
              gCCGCGTTTAATGACCCCCAAGGCCGGTATGCGTTTTGCCTGAACTTCGCCAACATGCGTTAAGTCTTGTACTGGGTGGCCGACGATCTGTCGGcc                                                                                          <  1:590881/96‑1 (MQ=255)
              gCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACat                                                                                      <  2:2255697/100‑1 (MQ=255)
              gCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACat                                                                                      >  2:2831857/1‑100 (MQ=255)
              gCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTGGCCTGACCTGAGCCAATATGCGTTAAGTCTTGTACTGAGTGGCCGACGATCTGTCGGGcactc                                                                                      <  1:1616025/100‑3 (MQ=255)
              gCCGCGTTTAATGACCCACAAAACCGTTAGGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACGATCTGTCGGCcactc                                                                                      <  1:329221/100‑3 (MQ=255)
                cgcgTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACatta                                                                                    >  2:1205863/1‑100 (MQ=255)
                cgcgTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACGATCTGTCGGCcactcag                                                                                    <  2:1182224/100‑5 (MQ=38)
                  cgTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACattatt                                                                                  <  2:2415477/100‑1 (MQ=255)
                  cgTTTAATGACCCACAAAACCGATATGCGTTTTGACTGACATGCGCCAACATGCGTTAAATCGTGTACTCAGTGGCCGACAGATCGTCGGGCACatta                                                                                    >  1:1370237/1‑98 (MQ=255)
                    tttcacgccccACAAAACCGTTATGCGTTTCGCCTGACCTACGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACATATCGTCGGCCACATTAtttt                                                                                <  1:1095405/92‑1 (MQ=255)
                      tAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTAtttttt                                                                              <  1:624428/100‑1 (MQ=255)
                       aaTGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTa                                                                             <  2:2916232/100‑1 (MQ=255)
                       aaTGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTa                                                                             <  2:2242628/100‑1 (MQ=255)
                        acgACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTAc                                                                            <  1:221301/98‑1 (MQ=255)
                        aTGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAATGGCCGACAGATCGTCGGCCACATTATTTTTTAc                                                                            >  1:1535774/1‑100 (MQ=255)
                            cccACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTcg                                                                        >  1:2633635/1‑100 (MQ=255)
                            cccACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTcg                                                                        >  1:2707572/1‑100 (MQ=255)
                            cccACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTcg                                                                        >  1:1165154/1‑100 (MQ=255)
                              cacaAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTcgac                                                                      <  2:2152774/100‑1 (MQ=255)
                              cacaAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTcgac                                                                      <  2:2222780/100‑1 (MQ=255)
                              cacaAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTcgac                                                                      >  2:1323264/1‑100 (MQ=255)
                              cacaAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTcgac                                                                      >  2:2691182/1‑100 (MQ=255)
                              cacaAAACCGTTATGCGTTGTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTcgac                                                                      <  1:891143/100‑1 (MQ=255)
                               acaAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTcgacg                                                                     <  1:2569297/100‑1 (MQ=255)
                                 aaaaCCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGaa                                                                   >  2:1903480/1‑100 (MQ=255)
                                  aaaCCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAAt                                                                  <  1:1121233/100‑1 (MQ=255)
                                     ccGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCct                                                               >  2:380368/1‑100 (MQ=255)
                                       nttATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCctct                                                             >  1:1942224/2‑100 (MQ=255)
                                       gTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCctct                                                             <  1:1501066/100‑1 (MQ=255)
                                       gTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCctct                                                             <  2:207938/100‑1 (MQ=255)
                                         tATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTc                                                           <  1:2715888/100‑1 (MQ=255)
                                           tGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTccc                                                         >  2:2049218/1‑100 (MQ=255)
                                             cGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCCGCCACATTATTTTTTACGTCGACGAATCCTCTTccc                                                         <  1:1203380/98‑1 (MQ=255)
                                              gTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCt                                                      >  1:896097/1‑100 (MQ=255)
                                                  tGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGttt                                                  >  1:261335/1‑100 (MQ=255)
                                                   gCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTc                                                 >  2:199296/1‑100 (MQ=255)
                                                   gCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTc                                                 >  1:134674/1‑100 (MQ=255)
                                                      tGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGcc                                              <  1:1579423/100‑1 (MQ=255)
                                                       gACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGccc                                             <  2:252175/100‑1 (MQ=255)
                                                         ccTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCa                                           >  2:2891732/1‑100 (MQ=255)
                                                         ccTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCa                                           >  2:1347774/1‑100 (MQ=255)
                                                          cTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCa                                           <  1:1324951/99‑1 (MQ=255)
                                                           tGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCata                                         >  2:1009772/1‑100 (MQ=255)
                                                            gcgcCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCatat                                        >  2:2324583/1‑100 (MQ=255)
                                                            gcgcCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCatat                                        >  2:1671906/1‑100 (MQ=255)
                                                              gtcAACATGCGTTCAGTCTTGTACTGAGTGGCCGACAGATCGACGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATAttt                                      >  2:1974170/3‑100 (MQ=255)
                                                              gcCAACATGCGGTAAGTCTTGCACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCGTCTTCCCGCTGTTTCGCTCCATAtt                                       >  1:2178048/1‑99 (MQ=255)
                                                               ccAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTc                                     <  2:1141155/100‑1 (MQ=255)
                                                               ccAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTc                                     >  1:573885/1‑100 (MQ=255)
                                                                 aaCATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCt                                   >  2:465096/1‑100 (MQ=255)
                                                                 aaCATGCGTGCAGTCTTGTACTGCGTGGTCGACAGATGGTCGGCCACATTATTTTTTACTTCGTCGAATCCTCTTCTAGCTGTTTCGCTCCGTATTTc                                     >  1:1208592/1‑98 (MQ=255)
                                                                   cATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTa                                 <  2:1872168/100‑1 (MQ=255)
                                                                     tGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCTGCCACATTATTTTTGACGTCGACGAATCCTCTTCCCGCTGGTTCGCCCCATATTTCCTTATc                               <  1:2712085/100‑1 (MQ=255)
                                                                      gCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATcc                              >  1:1504284/1‑100 (MQ=255)
                                                                        gTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGc                            >  2:603368/1‑100 (MQ=255)
                                                                        gTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGc                            >  2:575996/1‑100 (MQ=255)
                                                                          tAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCAt                          <  1:1395540/100‑1 (MQ=255)
                                                                           aaGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATa                         >  2:1302901/1‑100 (MQ=255)
                                                                             gTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGt                       >  1:2182717/1‑100 (MQ=255)
                                                                             gTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCGCTTCCCGCTGTTTCGCCCCATATTTCATTATCCGCATAGt                       >  1:699767/1‑100 (MQ=255)
                                                                               cTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCAATATTTCCTTATCCGCATAGTAt                     >  2:1616568/1‑100 (MQ=255)
                                                                                ttGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATc                    >  2:2693063/1‑100 (MQ=255)
                                                                                 tGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATcn                   <  2:861106/100‑2 (MQ=255)
                                                                                 tGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCa                   <  2:502534/100‑1 (MQ=255)
                                                                                 tGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCa                   <  2:573888/100‑1 (MQ=255)
                                                                                  gTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAg                  >  2:2388427/1‑100 (MQ=255)
                                                                                   tACTGAGTGGCCGACAGCTCGTCGGCCACATTTTTTTTTACGTCTACGAATGCTCTTCCCGCTGCGTCCCTCCATATTTCCTTATCCGCATCGTATCAgg                 >  1:405312/1‑100 (MQ=255)
                                                                                   aataatgTGGCCGACGATCTGTTGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAgg                 >  1:2537014/7‑100 (MQ=255)
                                                                                    ataatgTGGCCGACGATCTGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAgg                 >  2:109860/6‑99 (MQ=38)
                                                                                    aCTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGt                <  1:2497549/100‑1 (MQ=255)
                                                                                      aatgTGGCCGACGATCTGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTgc              >  2:1490694/4‑100 (MQ=39)
                                                                                       gAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTgcg             >  2:2325233/1‑100 (MQ=255)
                                                                                       atgTGGCCGACGATCTGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTgcg             >  2:828166/3‑100 (MQ=255)
                                                                                        aGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTgcgc            <  2:2340237/100‑1 (MQ=255)
                                                                                         gTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCt           >  2:640594/1‑100 (MQ=255)
                                                                                         gTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATACTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCt           >  1:1184674/1‑100 (MQ=255)
                                                                                         gTGGCCGACGATCTGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCt           >  2:37076/1‑100 (MQ=255)
                                                                                          tGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTtcgc            >  1:336207/1‑98 (MQ=255)
                                                                                          tGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTc          >  2:1762058/1‑100 (MQ=255)
                                                                                           ggCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTcc         >  2:2400071/1‑100 (MQ=255)
                                                                                           ggCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTcc         >  1:2852787/1‑100 (MQ=255)
                                                                                             ccGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTcccc       >  2:2137076/1‑100 (MQ=255)
                                                                                             ccGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTcccc       >  2:94282/1‑100 (MQ=255)
                                                                                              cGACAGATCGTCGGCCACATTATTTTTTACGTCGACTAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATGCGCATAGTATCAGGTGCGCTccccc      >  1:523492/1‑100 (MQ=255)
                                                                                                  aGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTTCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGcc  <  1:2049207/100‑1 (MQ=255)
                                                                                                  aGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGcc  >  1:1996575/1‑100 (MQ=255)
                                                                                                  aGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGcc  <  2:86999/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: