Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,047,222 A→G 5.9% D86G (GAC→GGC)  polA → 5' to 3' DNA polymerase and 3' to 5'/5' to 3' exonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,047,2220AG 5.9% 307.4 / 7.1 102D86G (GAC→GGC) polA5' to 3' DNA polymerase and 3' to 5'/5' to 3' exonuclease
Reads supporting (aligned to +/- strand):  ref base A (45/51);  new base G (4/2);  total (49/53)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.57e-01

Only 100 of 104 total aligned reads displayed.

CGGTGCGATGTATGGTGTCCTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGC  >  NC_000913/4047073‑4047370
                                                                                                                                                     |                                                                                                                                                    
cGGTGCGATGTATGGTGTCCTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTCTGAACATTACAAATCACATCGCCCGCCAATGCCGGa                                                                                                                                                      <  1:905336/150‑1 (MQ=255)
 ggTGCGATGTATGGTGTCCTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAATCCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGAc                                                                                                                                                     >  1:720876/1‑150 (MQ=255)
 ggTGCGATGTATGGTGTCCTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGAc                                                                                                                                                     <  1:231569/150‑1 (MQ=255)
  gTGCGATGTATGGTGTCCTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATgccggccg                                                                                                                                                    >  1:1058515/1‑150 (MQ=255)
     cGATGTATGGTGTCCTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATc                                                                                                                                                 <  1:1662757/150‑1 (MQ=255)
     cGATGTATGGTGTCCTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATc                                                                                                                                                 <  2:1121754/150‑1 (MQ=255)
      gATGTATGGTGTCCTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGTCCGGCATTGGCGGGCGAtgt                                                                                                                                                <  1:1610957/150‑3 (MQ=255)
      gATGTATGGTGTCCTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCt                                                                                                                                                <  2:1612498/150‑1 (MQ=255)
         gTATGGTGTCCTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGTCGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGTCCGGCATTGGCGGGCGAtgtgat                                                                                                                                             <  1:477580/150‑3 (MQ=255)
         gTATGGTGTCCTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGCCGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCg                                                                                                                                             >  1:57655/1‑150 (MQ=255)
          tATGGTGTCCTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGt                                                                                                                                            >  1:1553912/1‑150 (MQ=255)
           aTGGTGTCCTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTg                                                                                                                                           <  2:1430690/150‑1 (MQ=255)
                  ccTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGCCGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAAt                                                                                                                                    >  1:206145/1‑150 (MQ=255)
                  ccTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAAt                                                                                                                                    >  2:547863/1‑150 (MQ=255)
                  ccTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAAt                                                                                                                                    >  2:1525449/1‑150 (MQ=255)
                   cTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATc                                                                                                                                   <  1:1022106/150‑1 (MQ=255)
                      aaCATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCATGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGaa                                                                                                                                >  2:1383549/1‑150 (MQ=255)
                      aaCATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGaa                                                                                                                                >  1:764481/1‑150 (MQ=255)
                       ccATGCCGCGCACTCAGCTCCTGCACTATCACCCGACCCCTGCAGCGGTGCGCATCGACGCCCAGGGAAAAACCTTCCGTGATGAACTGTTTGAACAATACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGGACAAATCGAAc                                                                                                                               <  2:1169769/149‑1 (MQ=255)
                         ttgctgcGCTGTCTGATCATTCAATATAAACCGACGCATGCAGCGGTGGTATTTGACTCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAAccc                                                                                                                             <  1:1422572/149‑1 (MQ=255)
                         aTGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGGGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAAccc                                                                                                                             <  2:1091653/150‑1 (MQ=255)
                            ctgcGCAGTCTGATCATGCTATATAAACCGACGCATGCAGCGGTGGTCATTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTg                                                                                                                          <  1:254360/150‑1 (MQ=255)
                             tgcGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGc                                                                                                                         >  1:712003/1‑150 (MQ=255)
                                gcAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCAcg                                                                                                                      <  1:766021/150‑1 (MQ=255)
                                 cAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCAcgc                                                                                                                     <  1:983857/150‑1 (MQ=255)
                                 cAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCAcgc                                                                                                                     <  1:396964/150‑1 (MQ=255)
                                    tCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGAt                                                                                                                  <  1:9586/150‑1 (MQ=255)
                                     cTGATCATGCAATATAAACCGACGGATGCAGCGGTGGGCTTTGACGCCAAGCGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATg                                                                                                                 <  2:903523/150‑1 (MQ=255)
                                      tGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGCACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATgg                                                                                                                >  2:1100003/1‑150 (MQ=255)
                                      tGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATgg                                                                                                                >  1:845299/1‑150 (MQ=255)
                                      tGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATgg                                                                                                                >  1:978267/1‑150 (MQ=255)
                                      tGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATgg                                                                                                                >  2:1684860/1‑150 (MQ=255)
                                       gATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGt                                                                                                               <  2:295699/150‑1 (MQ=255)
                                        aTCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGtt                                                                                                              <  2:451557/150‑1 (MQ=255)
                                             gCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGAGGGTTACAGc                                                                                                         >  1:135493/1‑150 (MQ=255)
                                                atatAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGAt                                                                                                      <  1:420846/150‑1 (MQ=255)
                                                    aaaCCGCCGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGACCATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGGACCATTGCACGCGATGGTTAAAGCGATgggc                                                                                                  >  2:225723/1‑149 (MQ=255)
                                                     aaCCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGAc                                                                                                 >  1:944281/1‑150 (MQ=255)
                                                     aaCCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGAc                                                                                                 >  1:1082584/1‑150 (MQ=255)
                                                         gACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGcc                                                                                             >  2:723824/1‑150 (MQ=255)
                                                           cGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCgc                                                                                           <  2:666349/150‑1 (MQ=255)
                                                            gCATGCAGCGGTTGTCTTTTACGCCAACGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCTTGCACAAATCGAACCCTTGCACGCGTTGGTTAAAGCGATGGGACTGCCgct                                                                                          <  2:752684/150‑1 (MQ=255)
                                                                  aGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGc                                                                                    <  1:1325266/150‑1 (MQ=255)
                                                                    cGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCgg                                                                                  <  1:443376/150‑1 (MQ=255)
                                                                     ggtggtCTTTGACGCCAACGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGt                                                                                 <  1:1573692/150‑1 (MQ=255)
                                                                      gtggtCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGtt                                                                                <  2:1436729/150‑1 (MQ=255)
                                                                        ggtCTTTGACGCCAAGGGATAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTc                                                                              >  2:1135351/1‑150 (MQ=255)
                                                                        ggtCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTc                                                                              <  2:620609/150‑1 (MQ=255)
                                                                        ggtCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTc                                                                              <  1:404903/150‑1 (MQ=255)
                                                                          tCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATTGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTg                                                                            <  2:1137820/150‑1 (MQ=255)
                                                                          tCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTg                                                                            >  1:1464739/1‑150 (MQ=255)
                                                                           cTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATTGTTAAAGCGATGGTACTGCCGCTGCTGGCGGTTTCTgg                                                                           >  2:259332/1‑150 (MQ=255)
                                                                           cTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTgg                                                                           <  1:1081279/150‑1 (MQ=255)
                                                                           cTTTGACGCCAAGGGAAAAACCTTTCGTGAGGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGAGGGGACTGCCGCTGCTGGCGGTTTCTgg                                                                           <  1:1071431/150‑1 (MQ=255)
                                                                             ttGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCTCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCg                                                                         <  2:773949/150‑1 (MQ=255)
                                                                             ttGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCg                                                                         <  1:26690/150‑1 (MQ=255)
                                                                             ttGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCg                                                                         <  1:1065373/150‑1 (MQ=255)
                                                                             ttGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCATGCACAAATCGAACCCTTGCACTCGATGGTTAAAGCGATCGGACTGCCGCTGCTGGAGGTTTCTgact                                                                         >  2:288403/1‑147 (MQ=255)
                                                                                    cAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGc                                                                  >  1:814515/1‑150 (MQ=255)
                                                                                    cAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGc                                                                  >  2:1315831/1‑150 (MQ=255)
                                                                                      aGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCgg                                                                >  2:1235472/1‑150 (MQ=255)
                                                                                       gggAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGga                                                               >  1:1435125/1‑150 (MQ=255)
                                                                                        ggAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGCACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGgac                                                              >  2:1149499/1‑150 (MQ=255)
                                                                                        ggAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGgac                                                              <  1:1525449/150‑1 (MQ=255)
                                                                                           aaaaCCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGCACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGgacgac                                                           >  2:727813/1‑150 (MQ=255)
                                                                                           aaaaCCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGgacgac                                                           <  1:61229/150‑1 (MQ=255)
                                                                                            aaaCCTTTCGTGATGAACTGTTTGAACATTACAAATCCCATCGCCCGACAATGACGGACGAACTGCGTGCAAAAATAGAAACCTTGAACGCGATGCTTAAAGCGATGGGAATGCCGCAGCTGGCGTTTTCTGGCGTAGAAGCGGacgccg                                                          >  2:1573625/1‑150 (MQ=255)
                                                                                             aaCCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGt                                                         <  1:1572303/150‑1 (MQ=255)
                                                                                             aaCCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGt                                                         <  1:1565317/150‑1 (MQ=255)
                                                                                             aaCCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGt                                                         <  1:1135351/150‑1 (MQ=255)
                                                                                                cTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTTGAAGCGGACGACGTTAt                                                      >  1:1586652/1‑150 (MQ=255)
                                                                                                  ttCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCg                                                    >  2:583743/1‑150 (MQ=255)
                                                                                                    cGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGt                                                  >  2:1210790/1‑150 (MQ=255)
                                                                                                     gTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTa                                                 <  2:944281/150‑1 (MQ=255)
                                                                                                         tgaACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTActct                                             >  2:1387465/1‑150 (MQ=255)
                                                                                                          gaACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTg                                            <  1:775413/150‑1 (MQ=255)
                                                                                                             cTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGgcg                                         <  1:723824/150‑1 (MQ=255)
                                                                                                              tGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGgcgc                                        <  2:656126/150‑1 (MQ=255)
                                                                                                               gTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGgcgcg                                       >  1:156297/1‑150 (MQ=255)
                                                                                                                  tGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTTTCGGTACTCTGGCGCGCGa                                    >  2:702304/1‑150 (MQ=255)
                                                                                                                  tGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGa                                    >  2:197328/1‑150 (MQ=255)
                                                                                                                   gAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGaa                                   >  2:1072754/1‑150 (MQ=255)
                                                                                                                      cATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGcc                                >  1:571337/1‑150 (MQ=255)
                                                                                                                        ttACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGa                              <  1:1313353/150‑1 (MQ=255)
                                                                                                                            aaaTCACCTCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGaaaaa                          <  1:266547/150‑1 (MQ=255)
                                                                                                                             aaTCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAg                         >  1:810171/1‑150 (MQ=255)
                                                                                                                               gcagaTCGTCCGGCATTGGCGGGCGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGcc                       >  1:32230/5‑150 (MQ=255)
                                                                                                                                cagaTCGTCCGGCATTGGCGGGCGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGCCGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCg                      >  2:1411352/4‑150 (MQ=255)
                                                                                                                                cagaTCGTCCGGCATTGGCGGGCGATCTGCGTGCACAAATCGAACCCTTGCACCCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGTAATAGCCg                      >  2:144095/4‑150 (MQ=255)
                                                                                                                                  caTCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCggg                    >  1:51236/1‑150 (MQ=255)
                                                                                                                                   aTCGTCCGGCATTGGCGGGCGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGc                   >  2:732635/1‑150 (MQ=255)
                                                                                                                                       accgccAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTcc               >  1:333186/2‑150 (MQ=255)
                                                                                                                                              aTGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTg        <  2:1464739/150‑1 (MQ=255)
                                                                                                                                               tGCCGGACGATCTGCGTGCACAAATCGACCCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGa       <  1:1614619/150‑1 (MQ=255)
                                                                                                                                               tGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGa       >  2:1279817/1‑150 (MQ=255)
                                                                                                                                                 ccGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATc     <  1:1235472/150‑1 (MQ=255)
                                                                                                                                                   ggACGATCTGCGGGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCGGGCGCGCGAGGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAg   <  1:574039/150‑1 (MQ=255)
                                                                                                                                                   cgACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAg   <  1:1420469/149‑1 (MQ=255)
                                                                                                                                                    gACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGc  >  2:911960/1‑150 (MQ=255)
                                                                                                                                                    gACGAGCTGCGTGCACAAATCGCACCCTGGCACGCGATTGTTAAAGCGAGGGGACTGCCGCTGCTGGCGGTTGCTGGCGTAGAAGCTGACGACGTGATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGc  <  1:89445/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGGTGCGATGTATGGTGTCCTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGC  >  NC_000913/4047073‑4047370

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: