Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,382,6130AT23.6% 120.9 / 13.9 74L138F (TTA→TTTdegSserine endoprotease
Reads supporting (aligned to +/- strand):  ref base A (18/37);  new base T (14/3);  total (33/41)
Fisher's exact test for biased strand distribution p-value = 5.33e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGCGGTGGTTAACGTTTACAACCGTGGTTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTACCTACCATTCCAATTAATGCACGTCGCGTACCGCACATTGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCAAAACTTCCTCCAAACCGATGCTTCCATTAACCACGGTAACTCTGGCGGCGCGCTGGTGAACTCGCTGGGCGAACTGATGGGCATTAATACGCTGTCGTT  >  NC_000913/3382367‑3382858
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                
cgCGGTTGTTAACGTTAACAATCGTGCTTGGAACTCCAACTCTCACAACCAGCTTGACATCCTCACCCTGGGAACTTGTGCAATTATGGATCAACGCGGTTATTTCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTTCAGGATTGACGTGTATTTGAAGCATTGCTGGTGGGATCTGTCTCTCTAACCAATCCGGCGGTCCTTAAAATTAAGGCCACTCGCGGTTtaacta                                                                                                                                                                                                                                                             <  2:106561/250‑1 (MQ=255)
   ggtggtTAACGTTTACAACCGTGGTTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTtaactaccc                                                                                                                                                                                                                                                         >  1:686775/1‑250 (MQ=255)
     tggtTAACGTTTACAACCGTGGTTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGGTTTCCCACCAtt                                                                                                                                                                                                                                                      >  2:384643/1‑251 (MQ=255)
            cGTTTACATCCGGTGTTTGAACACCAACTCTCACAACCATCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTACCTACCATTCCAATTa                                                                                                                                                                                                                                               <  1:250053/251‑1 (MQ=255)
                  cAACCGTGGTTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTTCCTTCCCATTCCAATTattcac                                                                                                                                                                                                                                          >  2:296188/1‑246 (MQ=255)
                      cgttgTTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCTGTTATATCATCACCAATAAACACGGCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTACCTACCATTCCAATTAATGCACGTcg                                                                                                                                                                                                                                    <  2:88626/247‑1 (MQ=255)
                       gTGGTTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTta                                                                                                                                                                                                                                                                  <  1:794848/224‑1 (MQ=255)
                          gTTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGGTTAACTACCATTCCAATTACTGCAccgcgcgt                                                                                                                                                                                                                                  >  1:572319/1‑243 (MQ=255)
                           tttGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTTCCTCCCATTCCAATTATTGCTCGTcccggac                                                                                                                                                                                                                               >  1:793079/1‑245 (MQ=255)
                               aaCACCAACTCTCCCAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCAACACCAATAAACACGTCATCAACTACGCCTATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGACGCATTGCTGGTGGGATCTGACTCTCTAACCGATCAGGCGGTACTTAAAATGAATGCCACTTGCGGTTTACCTTCCATTCCAATTAATGCACGTCGCGTACCGca                                                                                                                                                                                                                           <  2:898576/251‑1 (MQ=255)
                                acacCAACTCTCACAACCAGCTTGAGTTCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCAGTGCTGGTGGGATCTGTCTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCTCTGGCGGTTTACCTACCATTCCAATTAATGCACGTCGCGTACCGcac                                                                                                                                                                                                                          <  2:886729/251‑1 (MQ=255)
                                            acaACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTACCTACCCATTCCATTAATTCACGTcccgcaacacaccttggcggcctg                                                                                                                                                                                                              >  1:230042/1‑228 (MQ=255)
                                            acaACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTACCTACCATTCCCATTTATGCACGTCCCGGACCCCACAtttgggacgc                                                                                                                                                                                                               >  1:572852/1‑242 (MQ=255)
                                                  aGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTACCTACCATTCCAATTAATGCACGTCGCGTACCGCACATTGGCGACGTAGTACTg                                                                                                                                                                                                        <  1:867664/251‑1 (MQ=255)
                                                       gagaTCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTACCTACCATTCCCATTAATTCACGTCCCGTACCCCACATTTGCGGCGTAGTACTGGCggt                                                                                                                                                                                                   >  1:55661/1‑249 (MQ=255)
                                                        agaTCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCATCATCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTTACTACCCTTCCCATTAATAGTCGTCCCGGACCCCCCATTGGCGACGTAGTACTGGCgctt                                                                                                                                                                                                  >  2:190408/1‑248 (MQ=255)
                                                         gaTCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTACCTACCCATCCAATTAATGCACGTCGCGGACCGCCCATTTGCGACGGAGGACTGGCGATCg                                                                                                                                                                                                >  1:806744/1‑251 (MQ=255)
                                                                aCCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTACCTACCATTCCCATTAATTCAAGTCGCGTACCGCACATTGGCGACGTAGTAATAGCGATCGGGAAccc                                                                                                                                                                                         >  1:180696/1‑251 (MQ=255)
                                                                aCCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTACCTACCATTCCAATTAATGCACCTCGCGTACCGCACCTTGGCGACCGAGTACTGGCGATCGGTaaacc                                                                                                                                                                                         >  1:441886/1‑248 (MQ=255)
                                                                     gggATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTACCTACCATTCCAATTAATGCACGTCGCGTACCGCACATTGGCGACGTAGTACTGGCGATCGGTAACCCGTa                                                                                                                                                                                      <  1:573733/249‑1 (MQ=255)
                                                                      ggATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTACCTACCATTCCCATTAATGTACGTCCCGTCCCGCCCACTGGCGACGTAGTACTGGCGATCGGTAACCCCTACa                                                                                                                                                                                    >  1:297385/1‑250 (MQ=255)
                                                                      ggATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTACCTACCATTCCAATTAATGCACGTCGCGTACCGCTCATTGGCGACGTAGTACTGGCGATCGGTAACCCGTACaa                                                                                                                                                                                   <  1:753725/251‑1 (MQ=255)
                                                                             ttgtAATGATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCATCGACGCCGATAAGATCAGCGTCTCCTTACAGGATGTACGTGTATTTGAAGCATAGCTGGTGGGATCTGACGCTCTAACCGATCTGGTTGTACTTAAAATTAATGCCACTTGCGGTTTACCTACCATTCCAATTAATGCCCGTCGCGGTCCGCACATTGGAGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCg                                                                                                                                                                             <  2:690916/248‑1 (MQ=255)
                                                                              tgtAATCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTTCCTACCATTCCAATTAATTAACGTCGCGTCACGCCCATTGGCGACGTAATACTGGCGGTTGGTaaacccgacacccccaggg                                                                                                                                                                           >  1:377319/1‑234 (MQ=255)
                                                                                 aaTCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTTCCTCCCCATCCAATTAATGCACGTCCCGTACAGCACATTGGCGGCGTAGTACTGGCGATCGGTAACCCGTACAAACTCGGGc                                                                                                                                                                          >  1:538072/1‑248 (MQ=255)
                                                                                 aaTCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTACCTACCATTCCAATTAATGCACGTCGCGTACCGCACATTGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAg                                                                                                                                                                        <  1:779226/250‑1 (MQ=255)
                                                                                  aTCATGGATCAACGCGGTTATATCATCATCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTACCTACCATTCCAATTAATGCACGTCGCGTACCGCACATTGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGAc                                                                                                                                                                     <  1:190408/251‑1 (MQ=255)
                                                                                   tCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTTACTCCCATTCCACTTAATACACGTCGCGTACCCCACATTGGCGGCGGAGTACTGGGGAACGGTAACCCCTACAACCTCGGGCCGAcc                                                                                                                                                                    >  1:889772/1‑251 (MQ=255)
                                                                                                  gTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTACCTACCATTCCCATTAATGCACGGCGCGTACCGCACATTGGAGGCGTAATACTTGCGATCGGGAACCCCTTCAACCCCGGGCAgccacttaaccacggggtt                                                                                                                                                      >  1:711675/1‑233 (MQ=255)
                                                                                                  gTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTACCTACCATTCCAATTAATGCACGTCGCGTACCGCACATTGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGatt                                                                                                                                                   <  1:384643/251‑1 (MQ=255)
                                                                                                      tatCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTTCCTACCCTTCCCCTTAATGTCCCCCGCTTACCGCACCTTTGCGACGGCGTACTGGCGCGCCGGAAACCGCACAACCTCCGGGCGACCAATACCCAgaacgaagagg                                                                                                                                               >  2:283294/1‑241 (MQ=255)
                                                                                                      tatCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTACCTACCATTCCAATTAATGCACGTCGCGTACCGCACATTGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTAcc                                                                                                                                                              <  1:283294/238‑1 (MQ=255)
                                                                                                         catcaCCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTTCCTCCCCTTCCCATTAATGCACGTCCCGTACCGCACATTGGCGACGGAGTACTGGCGTTCGGTAACCCGGACAACCCCCGGCAGGCCATTACCCCGGGGGTTATTAGTg                                                                                                                                            >  1:725386/1‑251 (MQ=255)
                                                                                                         catcaCCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTACCTACCATTCCAATTAATGCACGTCGCGTACCGCACATTGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTg                                                                                                                                            <  2:564836/251‑1 (MQ=255)
                                                                                                         catcaCCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTACCTACCATTCCAATTAATGCACGTCGCGTACCGCACATTGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTg                                                                                                                                            <  1:743831/251‑1 (MQ=255)
                                                                                                            caCCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTACCTACCATTCCAATTAATGCACGTCGCGTACCGCACATTGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCa                                                                                                                                         <  1:33501/251‑1 (MQ=255)
                                                                                                                 aTAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTACCTACCATTCCAATTAATGCACGTCGCGTACCGCACATTGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGt                                                                                                                                    <  1:317796/251‑1 (MQ=255)
                                                                                                                   aaaCACGTCATCAACGACGCCGATCAGATCAGCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTACCTACCATTCCAATTAATGCACGTCGCGTACCGCACATTGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCg                                                                                                                                  <  1:132278/251‑1 (MQ=255)
                                                                                                                       agttcctccaagtcgccgtgcagatcaTCGTCGCCTTTCAGGATGTATGGGGATTTGAGTCATTGCGGGTGGGATATGACACTCTAAACGAGCCGGCGGCACTTAAAATTAATGCGACTGGCGGTGTACCAACCATTCCAATTAATGCATGCCGCGTACCTCACATTGGCGACGTAGTACTGGCGATCGGTAACCCGTGCAACCTCGGTCAGACCATTACCCAGGGGATTATTAGTGCTACTGGTCGAATc                                                                                                                              <  2:145103/231‑1 (MQ=255)
                                                                                                                       acGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTACCTACCATTCCAATTAATGCACGTCGCGTACCGCACATTGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATc                                                                                                                              <  1:96691/251‑1 (MQ=255)
                                                                                                                        cggcatcaACGACGCCGATCAGATCATCGTCGCCTTACAGGAGGGACGTGTATTTGAAGCATGGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTACCTACCATTCCAATTAATGCACGTCGCGTACCGCACATTGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCg                                                                                                                             <  2:216496/248‑1 (MQ=255)
                                                                                                                             tcaACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGACGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGCTTACCTACCCTTTCAATTATTTCCATGCGGGTACCCCCCATTGGCGACGGATTAGTAGCGAGCGGGAAACCGTGAAACCACGGGCCCACCATTACCCAGGGGagtatttgctatacgggtgccgtagggcgt                                                                                                                        >  1:447408/1‑227 (MQ=255)
                                                                                                                              caACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGGTTACCTAACCATCCCATTAATACACGTCGCGTAACGCACAATGGCGACGTAGAAACGGCGATCGGGAAACCGTACCACCCACGGCAGACCAGTACCCAAGGGGATATTAggagcaccgggcggagatttttt                                                                                                                        >  1:291388/1‑228 (MQ=255)
                                                                                                                                   acggcgaacatatcatcTTCGCCGTACAGTATGGCCGTGTATTTGAAGCATTGCGGGTGGGATCTGACGCTCTAACCGATCTGGCGGTCCTTAAAATTAATGCCCCTGGCGGTTTACCTACCATTCCAATTAATGCGCGTCGCGTCCCGCACATTGGCGACGTAGTACTGGCGAGCGGTAACCCGTCCATCCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCg                                                                                                                  <  1:353451/240‑1 (MQ=255)
                                                                                                                                   acgCCGATCAGATCATCGTCGCCTTACAGGATGGACGTTTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTACCTACCATTCCAATTAATGCACGTCGCGTACCGCACATTGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCg                                                                                                                  <  2:180696/251‑1 (MQ=255)
                                                                                                                                      cctatcagatcaACGTCGCCTTTCAGGATGGATGTGTATTTGAAGCATTGCTGGTGGGATCTGAGTCTCTAACCGACCTGGCGGTACTGAAAATTAATGCCACTGGCTGTTTACATACCATGCCAATTAATGCACGTCGCGTACCGCACATTGGCGTCGTAGTACTGGCGAACGGTAATCCGTCCATCCTCGGGCAGACCATTGCCCAGGGGATTATTAGTGCCACGGGTCGAACCGGTCTGAACCCGAcc                                                                                                              <  2:714467/248‑1 (MQ=255)
                                                                                                                                       cGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTACCTACCATTCCAATTAATGCACGTCGCGTACCGCACATTGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCg                                                                                                             <  2:793079/251‑1 (MQ=255)
                                                                                                                                           cagatcaTCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTACCTACCATTCCAATTAATGCACGTCGCGTACCGCACATTGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGggcg                                                                                                         <  1:712123/251‑1 (MQ=255)
                                                                                                                                           cagatcaGCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCGAACCGATCTGGCTGTACTTAAAATTAATGCCACTGGCGGTTTACCTACCAGTCCAATTAATGCACGTCGCGTACCGCACATTGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGggc                                                                                                          <  2:693700/250‑1 (MQ=255)
                                                                                                                                                   cgtcgCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTTCCTCCCATTCCAATTAATACACGTAGCGTACCGCCCATTTGCGACGTAGTACTTGCGGTCGGGAACCCGCACAACCTCGGGCAGACCATTTCCCCGGGGATTTTTATTGCCACGGGTCGAATCGGGCTAAACCCCAACGGGCGGCACAACt                                                                                                 >  1:74362/1‑251 (MQ=255)
                                                                                                                                                          ttACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTACCTACCATTCCAAATAATGCACGTCGCGTAACGCACATTGGCGACGTAGGACTTGCGATCGGTAACCCGGACAAACTCGGGGAGACCATTACCCCGGGGATTATTAATGCCACGGGGCGAAACGGTCTGAACCCGACCCGGCGGGAAAAATtcccccc                                                                                          >  1:833598/1‑247 (MQ=255)
                                                                                                                                                           tACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTACCTCACATTCCAATTAATTCACGTCCCGTAACGCACATTTGCGCCGTAGGACTGGCGAACGGGAAACCGGACAACCACGGGCAGAACCATAACCACGGGATTATTAGAGGCACGGGGCGAAACGGGCTGAccccccacgggggggagaaattccccca                                                                                           >  2:824607/1‑226 (MQ=255)
                                                                                                                                                            aCAGGATGGACTTGTATTTGAAGCATTGCTGGTGGGATCTGACGCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTACCTACCAGTCCAATTAATGCACGTCGCGTACCGCACATTGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCAAAACTTCCTCCaa                                                                                         <  2:306277/250‑1 (MQ=255)
                                                                                                                                                                    gACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTTCCTACCATTCCCATTAATGCACGTCCCGTACCCCACATTGGCGACGGTGGAATGGCGATCGGTAACCCCGACAACCTCGGGCCGACCATTTCCCCGGGGATTATTTTTGCCACAGGGCGAATTGGTATTAAACCGCCCGGGGGGCCAAACTTCCTCCAAAACgctgct                                                                                >  1:235737/1‑251 (MQ=255)
                                                                                                                                                                      cGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGCTTACCCAACATTCCAATTAATGCACGTCGCGGACCGCACATTGGCGACGGAGTAATGGCGATCGGGAACCCGTACAACCTCGGGGAGACCATTACCCCAGGGGTTATTAGTGCCAAGGGGCGAAACGGACTGAAACCGACCGGGCGGGACAACTTTCTCCACACCggtggttc                                                                              >  1:556529/1‑244 (MQ=255)
                                                                                                                                                                           attgtaagaagtGCTGGTCGGATCTGCATCTCTAACCGATCGTGGTGAACGTAAATTTAATGCCACTGTCGGTTTCCCTTACATTGTAATTAAGGCACGCCGCGGTCCGCTCTTTGCCGACGTAGTACTGTCGTCTGTTTACCCGTACAACCTCGTGCAGCCCAATACACAGGGGATTATTAGTGCCACGGGTCGAATCGTTCTGAAACCGACTGGGCTGAAAAAATTCCTACAAACCGATGCTTCCATTa                                                                        <  2:362057/240‑1 (MQ=255)
                                                                                                                                                                               gAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTACCTACCCTTTCCATTAATTCAAGTCGCGTACCGCACCTTGGCGACGTAGTACTGGCGATCGGGAACCCGTACAACCTCGGGCAGACCACTACCCCGGGGGATATTAATGCCCAGGGGCGGATCGGGCTGAACCCGACCGGGCGGCAAAACTTCCCCCAAACCGATGCTTCCCTTAAccc                                                                    >  1:219406/1‑250 (MQ=255)
                                                                                                                                                                                      tGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTTCCTACCATTCCCATTAATGCACGTCGCGGACCGCCCATTGGCGACGTATTACTGGCGATCGGGAACCCCTACAACCTCGGGGCGAACATTACCCAGGGGGTTATTTATGCCCCGGGGCGAAACGGTCCGAACCCGACCGGGCGGCCAAACTTCCTCCCAACCGATGCTTCCATTAACCACGGTAAc                                                             >  1:159452/1‑251 (MQ=255)
                                                                                                                                                                                           gtggGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTACCTACCATTCCAATTAATGCACGTCGCGTACCGCACATTGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCAAAACTTCCTCCAAACCGATGCTTCCATTAACCACGGTAACTCTgg                                                        <  1:50604/251‑1 (MQ=255)
                                                                                                                                                                                             gggATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGGTTTCCCACCCTTCCAATTAATGCACGTCTCGCACCGCCCATTGGAGACGTAGTAATGGCGATCTGTAACCCCTTCAACCCCTGGGCGACCCTTACCCCGGGGGTTCTTAGTGCCCCCGGGCGAATCTGGCTGAACCCCCACGGGGCGCCAACCCTCCTCCAAACCGAAGCTTCCCTTAACCAccggcaatccttcg                                                      >  2:248508/1‑238 (MQ=255)
                                                                                                                                                                                                     aCTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTTCCTACCATTCCAATTAATGCACGGCGCGTACCGCACATTGGGGGCGTAGTAATGGCGGTCGGGAACCCCGACAACCTCGGGGAGACCATTACCCAGGGGGTTATTAGTGCCACGGGGCGGATCGGGCTGGAACCGCCCGGGGGGCCAAACTTCCTCCCAACCGATGCTTTCCTTAAACAACGGGACTCTGGGGGCGCGCTg                                              >  1:217322/1‑251 (MQ=255)
                                                                                                                                                                                                            aaCCGATCTGGCGGTACTTAAAATTAATGCCCCTGGCGGTTTACCTACCATTCCAAATAATGCACGTCGCGTACCGCACAGTGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCAAAACTGCCTCCAAACCGATGCTTCCATTAACCACGGTAACTCTGGCGGCGCGCTGGTGAACt                                       <  2:667814/251‑1 (MQ=255)
                                                                                                                                                                                                            aaCCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTACCTACCATTCCAATGAATGCACGTCGCGTACCGCACATTGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCAAAACTTCCTCCAAACCGATGCTTCCATTAACCACGGTAACTCTGGCGGCGCGCTGGTGAACt                                       <  2:973075/251‑1 (MQ=255)
                                                                                                                                                                                                              ccgtacTGGCTGTACTTAAAATTAAGGCCCCTGGCGGTTTACCTACCATTCCAATTAATGCACGTCGCGTACCGCACATGGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCATGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCAAAACTTCCTCCAAACCGATGCTTCCATTAACCACGGTAACTCTGGCGGCGCGCTGGTGAACTCg                                     <  2:960176/246‑1 (MQ=255)
                                                                                                                                                                                                                gATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTACCTACCATTCCAATTAATACACGTCCCGTACCGCCCATTGGCGACGGTGTACTGGCGGTCGGGAACCCGTACAACCCCGGGCAGACCATTACCCCGGGGATTTTTAGTGGCCCGGGTCGGATCGGGCTTAACCCCACCGGGCGGCCAAACTTCCTCCAAAACGATGCTTCCCTTTACCCCGGTAACTCTCTAGGCGCGCTgggggaatctcc                                   >  1:703513/1‑241 (MQ=255)
                                                                                                                                                                                                                     aggaggTACTGAAAATTACTGCCACTGGCGGTTTACCTACCAGTCCAATTAATGCACGGTGCGTTCCGCACATTGGCGACGTAGTACTGGCGAGAGGTAACCCGTACAACCTCAGACAGACCATTACCCAGGGGATTATTAGCTCCACGGGTCGAATCTGTCTGAACCGGACCGGGCGGCACAACTTCCTCCAAACCGATGCTGCCATTAACCGCGGGAACTCTGGCGGCGCGCTGGTGAACTCGCTgttc                               <  2:9567/247‑4 (MQ=255)
                                                                                                                                                                                                                                aaaaTTAATGCCACTGGCGGTTTACCTACCATTCCAATTAATGCACGTCGCGTACCGCACATTGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCAAAACTTCCTCCAAACCGATGCTTCCATTAACCACGGTAACTCTGGCGGCGCGCTGGTGAACTCGCTGGGCGAACTGATgg                     <  1:296188/250‑1 (MQ=255)
                                                                                                                                                                                                                                   actaaggCCACCGGCGGTTTTCTTACCAGTCCAATTAAAGCTCGGCGCGTACCGCATATGGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCTGGCAGACCATTACCCAGGGGAGTATTAGTGCCGCGGGGCGAATCGGTCTGAACCCGACCGGGCGGCAAAACTTCCACCAAGCCGATGCTTCCATTAACCACGGTAACGCTGGCTGCGCGCTGGTGAACTCGCTTGGTTAGCTGATGGGCAt                 <  2:122182/245‑1 (MQ=255)
                                                                                                                                                                                                                                    gtatatccACCGGCGGTTTTCCTACCATTCCAATTAATGCCCGGCGCGTACCGCCCATTGGCAACTTAGTACAGGCGATCGGTAACCTGTACAACCGCGGGCAGGCCATTACCCAGGGGATTATTAGTGCCACGTGTCGAATCGGTCTGAACCCGACCGGGCGGCAAAACTTCCTCCAAACCGATGCTTCCATTAACCACGGTAACTCTGGCGGCGCGCTGGTGAACTCGCTGGGCGAACTGATGGGAAtt                <  2:513129/245‑1 (MQ=255)
                                                                                                                                                                                                                                     tAATGCCACTGGCGGTTTACCTACCATTCCAATTAATGCACGTCGCGTACCGCACATTGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCAAAACTTCCTCCAAACCGATGCTTCCATTAACCACGGTAACTCTGGCGGCGCGCTGGTGAACTCGCTGGGCGAACTGATGGGCAtt                <  1:891939/250‑1 (MQ=255)
                                                                                                                                                                                                                                         gCCACTGGCGGTTTTCCTACCATTCCAATTAATGCACGTCGTGTGCCGCACTTTGTCGACGTATGACTGGCTATCGGTAACCCGTACAGTCTCGGGCAGACCAGTTCCCAGTGGATGAGGAGTGCCACGTGTCGAATCGGTCTTAACCCGACCGGGCGGTTAAGCTTTGTCCAATCCGATGCTTCCATTAACCTCGGTAGCTCTGGCGGCGCGCTGGTGAACTCGCTGGGCGAACTGATGGGCATTAATAc           <  1:132061/251‑1 (MQ=255)
                                                                                                                                                                                                                                          ccACTGGCGGTTTGCCTACCATTCCAATTAATGCACGTCGCGTACCGCACATTGGCGGCGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGGTTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCAAAACTTCCTCCAAACCGATGCTTCCATTTACCACGGGAACTCTGGCGGCGCGCTGGGGAACTCGCTGGGCGAACTGATGGGCATTAATACg          >  1:665244/1‑251 (MQ=255)
                                                                                                                                                                                                                                           cACTGGCGGTTTACCTACCATTCCAATTAATGCACGGCGCGTACCGCACATTGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCAAAACTTCCTCCAAACCGATGCTTCCATTAACCACGGTAACTCTGGCGGCGCGCTGGTGAACTCGCTGGGCGAACTGATGGGCATTAATACGc         <  1:62618/251‑1 (MQ=255)
                                                                                                                                                                                                                                               ggCGGTTTACCTACCATTCCAATTAATGCACGTCGCGTACCGCACATTGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCAAAACTTCCTCCAAACCGATGCTTCCATTAACCACGGTAACTCTGGCGGCGCGCTGGTGAACTCGCTGGGCGAACTGATGGGCATTAATCCGCTGtg     >  2:176974/1‑250 (MQ=255)
                                                                                                                                                                                                                                               ggCGGTTTACCTACCATTCCAATTAATGCACGTCGCGTACCGCACATTGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCAAAACTTCCTCCAAACCGATGCTTCCATTAACCACGGTAACTCTGGCGGCGCGCTGGTGAACTCGCTGGGCGAACTGATGGGCATTAATACGCTGTc     <  1:257943/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 cGGTTTACCTACCATTCCAATTAATGCACGTCGCGTACCGCACATTGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCAAAACTTCCTCCAAACCGATGCTTCCATTAACCACGGTAACTCTGGCGGCGCGCTGGTGAACTCGCTGGGCGAACTGATGGGCATTAATACGCTGTCGt   <  1:924006/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  ggTTTACCTACCATTCCAATTAATGCACGTCGCGTACCGCACATTGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCAAAACTTCCTCCAAACCGATGCTTCCATTAACCACGGTAACTCTGGCGGCGCGCTGGTGAACTCGCTGGGCGAACTGATGGGCATTAATACGCTGTCGtt  >  1:947807/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                
AGCGGTGGTTAACGTTTACAACCGTGGTTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGTTTACCTACCATTCCAATTAATGCACGTCGCGTACCGCACATTGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCAAAACTTCCTCCAAACCGATGCTTCCATTAACCACGGTAACTCTGGCGGCGCGCTGGTGAACTCGCTGGGCGAACTGATGGGCATTAATACGCTGTCGTT  >  NC_000913/3382367‑3382858

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: