Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,850,715 Δ1 bp 15.8% 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.15.8% 396.7 / 69.7 133intergenic (+22/‑145)sppA/ansAprotease IV (signal peptide peptidase)/cytoplasmic L‑asparaginase 1
Reads supporting (aligned to +/- strand):  ref base A (59/53);  new base . (7/14);  total (66/67)
Fisher's exact test for biased strand distribution p-value = 1.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01

Only 100 of 141 total aligned reads displayed.

GAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCC  >  NC_000913/1850616‑1850814
                                                                                                    |                                                                                                    
gaaagtgaCAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑a                                                                                                      >  1:148590/1‑100 (MQ=255)
gaaagtgaCAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑a                                                                                                      <  1:2066412/100‑1 (MQ=255)
gaaagtgaCAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATCCGTTAAGTCTTGTACTGAGTGGCCGAC‑a                                                                                                      >  2:2789221/1‑100 (MQ=255)
 aaagtgaCAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑Ag                                                                                                     <  1:272724/100‑1 (MQ=255)
 aaagtgaCAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑Ag                                                                                                     >  2:3561795/1‑100 (MQ=255)
    gtgaCAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGAtc                                                                                                  <  2:2141889/100‑1 (MQ=255)
        cAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGAtc‑gtcg                                                                                             >  2:1912038/1‑100 (MQ=255)
         aaGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCgg                                                                                            <  2:1695726/100‑1 (MQ=255)
           gCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGcc                                                                                          <  2:1547685/100‑1 (MQ=255)
           gCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑‑GATCTGTCGGcc                                                                                          <  2:2567333/100‑1 (MQ=255)
            cTGGCGGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑‑GATCTGTCGGCca                                                                                         <  2:3398078/100‑1 (MQ=255)
            cTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAAGAtc‑ggaagag                                                                                          >  1:3176805/1‑94 (MQ=255)
             tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCcac                                                                                        <  2:1064900/100‑1 (MQ=255)
             tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑‑GATCTGTCGGCcac                                                                                        <  2:317944/100‑1 (MQ=255)
             tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑‑GATCTGTCGGCcac                                                                                        <  1:3551784/100‑1 (MQ=255)
             tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑‑GATCTGTCGGCcac                                                                                        <  1:386230/100‑1 (MQ=255)
             tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑‑GATCTGTCGGCcac                                                                                        <  1:708488/100‑1 (MQ=255)
             tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACCGAGTGGCCGAC‑‑GATCTGTCGGCcac                                                                                        <  2:2624773/100‑1 (MQ=255)
               gCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACat                                                                                      >  2:1783693/1‑100 (MQ=255)
                ccGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACatt                                                                                     >  2:2640407/1‑100 (MQ=255)
                 cgcgTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCTCCAACATGCGCTAAGTCTTGTACTGAGTGGCCGAC‑‑GATCTGTCGGCcac                                                                                        <  2:2380646/96‑1 (MQ=255)
                 cgcgTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑‑GATCTGTCGGCcactcag                                                                                    <  1:660612/100‑5 (MQ=38)
                  gcgTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCCACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑‑GATCTGTCGGCcactcagt                                                                                   <  2:1884024/100‑6 (MQ=38)
                  gcgTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACattat                                                                                   >  1:2615423/1‑100 (MQ=255)
                    gtctattgtgccACGAATCCGTTGTGCGTGTCGCCTGTCGAGCGCCAGCTTGGGTGAACTCTTGTACTGAGGCGGCGAG‑ACATC‑GTCGGCCACATTAttt                                                                                 <  2:2964568/90‑1 (MQ=255)
                       tAATGACCCACAAGACCGTTATGCGTTTTGCCTGGCCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTAtttttt                                                                              <  2:3379388/100‑1 (MQ=255)
                       tAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTAtttttt                                                                              <  1:3228097/100‑1 (MQ=255)
                        aaTGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTa                                                                             <  1:3572713/100‑1 (MQ=255)
                        aaTGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTa                                                                             <  2:3623140/100‑1 (MQ=255)
                        aaTGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTa                                                                             <  1:2866851/100‑1 (MQ=255)
                         aTGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTAc                                                                            >  2:2034438/1‑100 (MQ=255)
                         aTGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTAc                                                                            >  1:1689243/1‑100 (MQ=255)
                         aTGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTAc                                                                            <  1:1912048/100‑1 (MQ=255)
                          tGACCCAAAAAACCGTTATGCGTTTTGCCTGACCTGCCCCAACATGCGTTAAGTGTTGTACTGAGTTGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTTCg                                                                           >  2:3366515/1‑100 (MQ=255)
                             cccACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTcg                                                                        >  1:3075363/1‑100 (MQ=255)
                             cccACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTcg                                                                        >  1:3321668/1‑100 (MQ=255)
                               cacaAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTcgac                                                                      >  1:3092489/1‑100 (MQ=255)
                               cacaAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTcgac                                                                      >  2:468503/1‑100 (MQ=255)
                               cacaAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTcgac                                                                      <  1:2829646/100‑1 (MQ=255)
                               cacaAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTcgac                                                                      >  1:210073/1‑100 (MQ=255)
                                acaAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTcgacg                                                                     <  1:1655443/100‑1 (MQ=255)
                                 caAAACCGTTATGCGTTTTGCCTGACCTGCTCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTcgacg                                                                     <  2:2464312/99‑1 (MQ=255)
                                 caAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTcgacga                                                                    >  2:3610111/1‑100 (MQ=255)
                                      ccGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTACGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATcc                                                                <  1:3523254/99‑1 (MQ=255)
                                       cGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCctc                                                              <  1:1101360/100‑1 (MQ=255)
                                       cGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCctc                                                              <  1:2951905/100‑1 (MQ=255)
                                       cGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCctc                                                              >  1:3081139/1‑100 (MQ=255)
                                        gTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCctct                                                             >  1:303780/1‑100 (MQ=255)
                                          tATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTc                                                           >  1:940895/1‑100 (MQ=255)
                                           aTGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTcc                                                          <  1:3105015/100‑1 (MQ=255)
                                            tGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTccc                                                         >  1:2655229/1‑100 (MQ=255)
                                              cGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTcctg                                                        >  1:3460852/1‑97 (MQ=255)
                                              cGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGc                                                       <  2:729224/100‑1 (MQ=255)
                                              cGTTTGGCCTGACCTGCGCCAACAGGCGCTAAGTCTTGTACAGAGTGGCCGAC‑AGATC‑GTCGGCGACATTATTTTTTACGTCGACGAATCCTCTTCCCGc                                                       >  2:1994216/1‑100 (MQ=255)
                                               gTTTTGCCTGACCTGCGCCAACATGCGTTAATTGTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCt                                                      <  2:1801080/100‑1 (MQ=255)
                                               gTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCt                                                      <  1:2764466/100‑1 (MQ=255)
                                                 tttGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCt                                                      <  1:1984145/98‑1 (MQ=255)
                                                 tttGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGt                                                    <  2:1516088/100‑1 (MQ=255)
                                                    gCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTc                                                 <  1:813343/100‑1 (MQ=255)
                                                    gCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTc                                                 >  1:101096/1‑100 (MQ=255)
                                                     ccTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCg                                                <  1:1036939/100‑1 (MQ=255)
                                                       tGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGcc                                              >  1:3018447/1‑100 (MQ=255)
                                                       tGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGcc                                              >  1:2725267/1‑100 (MQ=255)
                                                       tGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGcc                                              >  1:250229/1‑100 (MQ=255)
                                                       tGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGcc                                              >  1:1735145/1‑100 (MQ=255)
                                                       tGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGcc                                              <  1:1323232/100‑1 (MQ=255)
                                                          ccTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCa                                           >  1:1315428/1‑100 (MQ=255)
                                                              cgcCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATAtt                                       >  1:19591/1‑100 (MQ=255)
                                                              cgcCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATAtt                                       <  2:940888/100‑1 (MQ=255)
                                                                   aCATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCtt                                  >  2:993898/1‑100 (MQ=255)
                                                                      tGCGTTAAGTGTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTGCCCGCTGTTTCGCCCCATATTTCCTTATc                               <  2:3629425/100‑1 (MQ=255)
                                                                      tGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTAt                                >  2:1184735/1‑99 (MQ=255)
                                                                        cGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCg                             <  1:3370777/100‑1 (MQ=255)
                                                                         gTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCTTTTAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGc                            >  2:574128/1‑100 (MQ=255)
                                                                         gTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGc                            <  2:1299852/100‑1 (MQ=255)
                                                                         gTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGc                            >  2:2601148/1‑100 (MQ=255)
                                                                         gTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGc                            >  2:773762/1‑100 (MQ=255)
                                                                           tAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCAt                          <  1:2296824/100‑1 (MQ=255)
                                                                              gTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAga                       >  1:1224398/1‑99 (MQ=255)
                                                                              gTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGt                       >  2:1701680/1‑100 (MQ=255)
                                                                              gTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGt                       >  1:1480548/1‑100 (MQ=255)
                                                                                cTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTAt                     >  2:1284218/1‑100 (MQ=255)
                                                                                aaaaaataatgTGGCCGAC‑‑GATCTGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTAt                     >  2:419895/11‑100 (MQ=37)
                                                                                     aCTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGt                <  1:2614512/100‑1 (MQ=255)
                                                                                      cTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTg               <  2:667853/100‑1 (MQ=255)
                                                                                       tGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTgc              >  1:2282456/1‑100 (MQ=255)
                                                                                       aatgTGGCCGAC‑‑GATCTGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTgc              >  2:3540864/4‑100 (MQ=39)
                                                                                          gTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACTTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCt           >  1:638445/1‑100 (MQ=255)
                                                                                          gTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCt           <  2:2807620/100‑1 (MQ=255)
                                                                                          gTGGCCGAC‑‑GATCTGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCt           >  1:586872/1‑100 (MQ=255)
                                                                                          gTGGCCGAC‑‑GATCTGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCt           >  2:605523/1‑100 (MQ=255)
                                                                                           tGGCCGAC‑‑GATCTGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTc          >  1:2365484/1‑100 (MQ=255)
                                                                                            ggCCGAC‑‑GATCTGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTcc         >  2:1924644/1‑100 (MQ=255)
                                                                                             gCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTccc        >  1:391188/1‑100 (MQ=255)
                                                                                              ccGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTcccc       >  1:3445333/1‑100 (MQ=255)
                                                                                              ccGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTcccc       >  1:2042969/1‑100 (MQ=255)
                                                                                              ccGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTcccc       <  1:1561979/100‑1 (MQ=255)
                                                                                                  c‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGc   <  2:2538009/100‑1 (MQ=255)
                                                                                                    aGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGcc  <  1:3610138/100‑1 (MQ=255)
                                                                                                    aGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGcc  <  2:1821485/100‑1 (MQ=255)
                                                                                                    |                                                                                                    
GAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGAC‑AGATC‑GTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCC  >  NC_000913/1850616‑1850814

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