Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,950,292 C→G 5.4% A77A (GCG→GCC)  aspS ← aspartyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,950,2920CG 5.4% 379.4 / 7.6 123A77A (GCG→GCC) aspSaspartyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (62/54);  new base G (5/2);  total (67/56)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 131 total aligned reads displayed.

AGGTAGCGGTATTTCAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCA  >  NC_000913/1950143‑1950439
                                                                                                                                                     |                                                                                                                                                    
aGGTAGCGGTATTTCAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCAcgc                                                                                                                                                      <  2:1580668/150‑1 (MQ=255)
 ggTAGCGGTATTTCAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCAcgcg                                                                                                                                                     >  1:1940533/1‑150 (MQ=255)
   tAGCGGTATTTCAGACGCGCTTCTTCGGTGTTGTCGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCTGCGTGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACgcgct                                                                                                                                                   <  1:514833/150‑2 (MQ=255)
   tAGCGGTATTTCAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCa                                                                                                                                                   <  1:1047237/150‑1 (MQ=255)
   tAGCGGTATTTCAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCa                                                                                                                                                   <  1:1054435/150‑1 (MQ=255)
   gaGCGGTATTTAATTCGAGCTTCTTCGGGTTGGACGTGGGTAGAGTCAAGCGGCAGAACATCTGCGCTGGTGATGATAGTCAGCTAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCa                                                                                                                                                   <  2:187661/149‑1 (MQ=255)
    aGCGGTATTTCAGCCGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑c                                                                                                                                                 >  2:1181075/1‑150 (MQ=255)
     gCGGTATTTCAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑Cg                                                                                                                                                >  1:816237/1‑150 (MQ=255)
      cGGTATTTCAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGt                                                                                                                                               >  2:1065111/1‑150 (MQ=255)
       ggTATTTCAGCCGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTa                                                                                                                                              >  1:1374370/1‑150 (MQ=255)
       ggTATTTCAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTa                                                                                                                                              >  2:464143/1‑150 (MQ=255)
        gTATTTCAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTAc                                                                                                                                             <  1:2078916/150‑1 (MQ=255)
         tacttCAGACGCGCTTCTTCGGTGTTGACGTGGCTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATTGTCAGCGATGACGCCAGCACTTCGATTTCGCCGGTCGCCAAATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACg                                                                                                                                            <  1:1324570/147‑1 (MQ=255)
          aTTTCAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACgg                                                                                                                                           >  2:739719/1‑150 (MQ=255)
            ttCAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTg                                                                                                                                         >  1:1423814/1‑150 (MQ=255)
            ttCAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCAAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTg                                                                                                                                         >  1:1320243/1‑150 (MQ=255)
               aGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGccc                                                                                                                                      <  1:1896375/150‑1 (MQ=255)
                gACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCg                                                                                                                                     <  2:2168691/150‑1 (MQ=255)
                  cgcgCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTg                                                                                                                                   <  2:82241/150‑1 (MQ=255)
                    cgcTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGAc                                                                                                                                 <  2:1113437/150‑1 (MQ=255)
                    cgcTGCTTCGGCGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAc                                                                                                                                 <  2:294493/150‑1 (MQ=255)
                       ttcttcGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAcgaa                                                                                                                              <  1:1259472/150‑4 (MQ=255)
                          ttcGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAAt                                                                                                                           >  2:1754027/1‑150 (MQ=255)
                               tgtTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGa                                                                                                                      <  1:120757/150‑1 (MQ=255)
                               ggtTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTATTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGa                                                                                                                      <  1:1812657/149‑1 (MQ=255)
                                 ttGACGTGGTTAGAGTCACGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAAc                                                                                                                    <  2:2100901/150‑1 (MQ=255)
                                 ttGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAAc                                                                                                                    <  2:1798684/150‑1 (MQ=255)
                                   gACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTc                                                                                                                  >  1:780583/1‑150 (MQ=255)
                                     cGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCAt                                                                                                                >  1:898767/1‑150 (MQ=255)
                                     cGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCAt                                                                                                                >  1:1745198/1‑150 (MQ=255)
                                        gggtAGAGTCAAGCGGCAGAACAGCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTAc                                                                                                             <  1:888123/147‑1 (MQ=255)
                                        ggTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTAc                                                                                                             >  1:517355/1‑150 (MQ=255)
                                        ggTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATAGCCGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCAATAc                                                                                                             >  2:1170019/1‑150 (MQ=255)
                                         gTTAGAGTCAAGCGGCAGACCATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACg                                                                                                            <  1:1441491/150‑1 (MQ=255)
                                          ttAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGc                                                                                                           >  1:421887/1‑150 (MQ=255)
                                           tAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCa                                                                                                          >  2:557181/1‑150 (MQ=255)
                                             gagTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGt                                                                                                        >  1:1099908/1‑150 (MQ=255)
                                              agTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGtt                                                                                                       <  2:1744678/150‑1 (MQ=255)
                                                 cAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGCACTCATTACGCAGTTCag                                                                                                    >  2:1935773/1‑150 (MQ=255)
                                                 cAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCag                                                                                                    >  2:818814/1‑150 (MQ=255)
                                                 cAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCag                                                                                                    >  2:1922683/1‑150 (MQ=255)
                                                 aaccgcgCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCag                                                                                                    >  2:1441676/7‑150 (MQ=255)
                                                    gCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAgg                                                                                                 <  1:809650/150‑1 (MQ=255)
                                                      ggCAGAACTTCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGGCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGcc                                                                                               <  1:441214/150‑1 (MQ=255)
                                                       gCAGAACTTCTGCGCGGTTGATGATAGTGAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCTGTTAATAATTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCa                                                                                              <  1:652828/150‑1 (MQ=255)
                                                       gCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCa                                                                                              >  1:1742139/1‑150 (MQ=255)
                                                             cATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTa                                                                                        <  2:816788/150‑1 (MQ=255)
                                                                cTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGTCCAGCTTTAAcg                                                                                     >  1:240320/1‑150 (MQ=255)
                                                                cTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAcg                                                                                     <  1:1665757/150‑1 (MQ=255)
                                                                  gcgcGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAcgcg                                                                                   >  2:1994573/1‑150 (MQ=255)
                                                                    gcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTc                                                                                 <  2:509998/150‑1 (MQ=255)
                                                                      ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCg                                                                               <  1:988856/150‑1 (MQ=255)
                                                                      ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCg                                                                               <  2:891579/150‑1 (MQ=255)
                                                                      ggTTGATGATAGTCAGCGAGGACACCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGGTCAGAGGCCAGCTTTAACGCGTCCg                                                                               <  2:156972/150‑1 (MQ=255)
                                                                         tgatgaTAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCAc                                                                            >  2:1383976/1‑150 (MQ=255)
                                                                         tgatgaTAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCAc                                                                            >  1:1069989/1‑150 (MQ=255)
                                                                          gatgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACg                                                                           >  2:2180195/1‑150 (MQ=255)
                                                                            tgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTACCGCGTCCGCACGAt                                                                         >  2:1401526/1‑150 (MQ=255)
                                                                            tgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGAt                                                                         <  1:1710032/150‑1 (MQ=255)
                                                                             gatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAATTCAGAGGCCAGCTTTAACGCGTCCGCACGATc                                                                        >  2:1543118/1‑150 (MQ=255)
                                                                              atAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATcc                                                                       <  1:2137897/150‑1 (MQ=255)
                                                                               tAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGCGGCCAGCTTTAACGCGTCCGCACGATCCg                                                                      >  1:1616036/1‑150 (MQ=255)
                                                                                aGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCgg                                                                     <  2:1099023/150‑1 (MQ=255)
                                                                                aGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCgg                                                                     <  1:709780/150‑1 (MQ=255)
                                                                                 gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGa                                                                    >  2:1816390/1‑150 (MQ=255)
                                                                                  tCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGAt                                                                   >  1:721559/1‑150 (MQ=255)
                                                                                   cAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATc                                                                  >  1:1268222/1‑150 (MQ=255)
                                                                                   cAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATc                                                                  <  2:603679/150‑1 (MQ=255)
                                                                                   cAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTAAGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATc                                                                  >  2:890401/1‑150 (MQ=255)
                                                                                    aGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCg                                                                 <  1:509112/150‑1 (MQ=255)
                                                                                     gCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGa                                                                >  2:2076567/1‑150 (MQ=255)
                                                                                         ggACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaaaaa                                                            >  1:1457652/1‑150 (MQ=255)
                                                                                            cGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc                                                         <  2:925883/150‑1 (MQ=255)
                                                                                              ccAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCt                                                       >  1:946314/1‑150 (MQ=255)
                                                                                                    cTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCCTTACGCTGGTCAGAGGCCAGCTTTAACGCGTCCGCACGAGCCGGAGCGAATAATTACTGCACGa                                                 >  2:1509788/1‑150 (MQ=255)
                                                                                                    cTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCCTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTTCGCAGTTCAGAGGCCAGCTTTAACGCGTCCGAACGATCCGTATCGAAAAATACCTGCACGa                                                 >  1:356523/1‑150 (MQ=255)
                                                                                                       cGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAGATACCTGCACGATAc                                              >  2:734510/1‑150 (MQ=255)
                                                                                                       cGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATAc                                              >  1:1040811/1‑150 (MQ=255)
                                                                                                             cGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTcgc                                        >  2:475003/1‑150 (MQ=255)
                                                                                                              gCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTcgcg                                       <  1:1608134/150‑1 (MQ=255)
                                                                                                                 ggTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGATCAGAGGCCAGCTTAAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTc                                    >  1:1124349/1‑150 (MQ=255)
                                                                                                                  gTCGCCATAACGAGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCTTTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTcg                                   <  1:166989/150‑1 (MQ=255)
                                                                                                                     gCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCgcgc                                >  1:2192597/1‑150 (MQ=255)
                                                                                                                     gCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCgcgc                                <  2:529291/150‑1 (MQ=255)
                                                                                                                       cATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCat                              >  1:1864485/1‑150 (MQ=255)
                                                                                                                          atCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATc                           >  2:1153439/1‑150 (MQ=255)
                                                                                                                            cgcgGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGa                         >  2:1289050/1‑150 (MQ=255)
                                                                                                                                gTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGaa                     <  2:1099908/150‑1 (MQ=255)
                                                                                                                                 ttAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCAATACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAg                    <  1:1894923/150‑1 (MQ=255)
                                                                                                                                  tAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGa                   <  2:1372131/150‑1 (MQ=255)
                                                                                                                                    atatTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAGATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATc                 <  1:734510/150‑1 (MQ=255)
                                                                                                                                      atTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAATAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAg               >  1:1611867/1‑150 (MQ=255)
                                                                                                                                        ttttttCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGc             <  1:1441676/150‑1 (MQ=255)
                                                                                                                                         tttttCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCTGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAACTACCTGCACGATACCTTCGCGGTCGCGCAGATCGATGAAGATCAGGcg            >  1:97214/1‑149 (MQ=255)
                                                                                                                                         tttttCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCt            >  2:1451894/1‑150 (MQ=255)
                                                                                                                                               gTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAAATCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCCCGATCCGGATCGAAAAATACCAGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAg      >  2:737401/1‑150 (MQ=255)
                                                                                                                                                tCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGa     >  2:2072543/1‑150 (MQ=255)
                                                                                                                                                 cACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGCGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAgct    >  1:44844/1‑148 (MQ=255)
                                                                                                                                                   cgcgCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCa  >  1:1298001/1‑150 (MQ=255)
                                                                                                                                                   agcgcA‑CGTACGGTGCCCGTGACGTGTATGAAGAACTCATTAAGCAGTTCACAGGCCAGCTTTAACGCGTCCGCACGATGAGGATCGAAAAATACATGCACGAAACCATCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCa  >  1:366130/2‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGGTAGCGGTATTTCAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCA  >  NC_000913/1950143‑1950439

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