Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913159,1010CA27.1% 220.9 / 28.7 96C9F (TGC→TTC) pcnBpoly(A) polymerase I
Reads supporting (aligned to +/- strand):  ref base C (53/17);  new base A (0/26);  total (53/43)
Fisher's exact test for biased strand distribution p-value = 2.05e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAGCATGCTGCTCACGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGGTAAGCCGCC  >  NC_000913/158952‑159250
                                                                                                                                                     |                                                                                                                                                     
ccGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAGCATGCTGCTCACGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGc                                                                                                                                                       <  2:1241984/150‑1 (MQ=255)
     ttATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAGCATGCTGCTCACGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCaaaaa                                                                                                                                                  >  1:1751194/1‑150 (MQ=255)
        ttGAGCCTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAGCATGCTGCTCACGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTa                                                                                                                                               <  2:461968/150‑1 (MQ=255)
        ttGAGCCTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAGCATGCTGCTCACGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTa                                                                                                                                               <  2:21491/150‑1 (MQ=255)
        ttGAGCCTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAGCATGCTGCTCACGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTa                                                                                                                                               <  1:2015928/150‑1 (MQ=255)
         tGAGCCTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAGCATGCTGCTCACGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAg                                                                                                                                              <  1:896903/150‑1 (MQ=255)
         tGAGCCTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAGCATGCTGCTCACGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAg                                                                                                                                              <  1:897129/150‑1 (MQ=255)
            gCCTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAGCATGCTGCTCACGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGa                                                                                                                                           >  1:133024/1‑150 (MQ=255)
             ccTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAGCATGCTGCTCACGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGAc                                                                                                                                          <  1:1891770/150‑1 (MQ=255)
                   cATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAGCATGCTGCTCACGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGt                                                                                                                                    >  1:513309/1‑150 (MQ=255)
                      tACCTTCCGGGAATTTTCACTGATATCTTTGCGGTAAATAGCATGCTGCTCACGGGGTATCACCGTCACCTGTGGACGGGCTACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGAAAAAATTAGCGACTCGGGTaaa                                                                                                                                 <  2:614272/150‑1 (MQ=255)
                      tACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAGCATGCTGCTCACGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTaaa                                                                                                                                 >  1:1840209/1‑150 (MQ=255)
                       aCCTTCTGGGCTTTTTCACTTATATCTTTGCGGGAAATAGCTTGCGGCTCACGCGGGATCACCGTCACCTGTGGAGGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGAAAAAATAAGCGACTCGGGTaaaa                                                                                                                                <  2:1030014/150‑1 (MQ=255)
                           tCAGGGCATTTTCACTGATATCTTTGCGGGAAATAGCATGCTGCTCACGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAg                                                                                                                            >  1:2099327/1‑150 (MQ=255)
                                   ttttCACTGATATCTTTGCGGGAAATAGCATGCTGCTCACGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTc                                                                                                                    >  1:2054686/1‑150 (MQ=255)
                                        aCTGATATCTTTGCGGGAAATAGCATGCTGCTCACGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAg                                                                                                               >  2:530664/1‑150 (MQ=255)
                                          tGATATCTTTGCGGGAAATAGCATGCTGCTCACGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAgtg                                                                                                             >  2:1477088/1‑150 (MQ=255)
                                            atatCTTTGCGGGAAATAGCATGCTGCTCACGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTc                                                                                                           >  1:2068512/1‑150 (MQ=255)
                                            atatCTTTGCGGGAAATAGCATGCTGCTCACGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTc                                                                                                           >  1:191625/1‑150 (MQ=255)
                                            atatCTTAGCGGGAAATAGCATGCTGCTCACGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTACTGTc                                                                                                           >  1:1282120/1‑150 (MQ=255)
                                                 tttGCGGGAAATAGCATGCTGCTCACGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAAca                                                                                                      >  1:1991391/1‑150 (MQ=255)
                                                 tttGCGGGAAATAGCATGCTGCTCACGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAAca                                                                                                      >  1:1052752/1‑150 (MQ=255)
                                                  ttGCGGGAAATAGCATGCTGCTCACGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAAcat                                                                                                     >  2:1328500/1‑150 (MQ=255)
                                                      gggAAATAGCATGCTGCTCACGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCAtt                                                                                                 <  2:1770070/150‑1 (MQ=255)
                                                      gggAAATAGCATGCTGCTCACGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCAtt                                                                                                 >  2:570663/1‑150 (MQ=255)
                                                      gggAAATAGCATGCTGCTCACGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCAtt                                                                                                 >  2:65912/1‑150 (MQ=255)
                                                         aaaTAGCTTGCTGCACACGGGGTAACACCGTCACCTGTGGAGGGGCGACTGCCTTTTCAGCCTCGATTTCCGCGGGGTTTACCACCTGGCGTAAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAg                                                                                              <  1:1075229/150‑1 (MQ=255)
                                                         aaaTAGCATGCTGCTCACGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAg                                                                                              <  2:1089903/150‑1 (MQ=255)
                                                         aaaTAGCATGCTGCTCACGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAg                                                                                              >  1:1289584/1‑150 (MQ=255)
                                                          aaTAGCATGCTGCTCACGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAgg                                                                                             >  1:1731670/1‑150 (MQ=255)
                                                            tAGCATGCTGCTCACGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGAc                                                                                           >  1:278330/1‑150 (MQ=255)
                                                             aGCATGCTGCTCACGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACa                                                                                          >  1:1199059/1‑150 (MQ=255)
                                                                 tgctgcTCAGGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGATTAGCACCTTGCGGAAAAATTTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACaaaaa                                                                                      <  2:336896/150‑1 (MQ=255)
                                                                  gctgctCCCGCGGGATCACCGTCACCTGTGTAGGGGCGACTGCGTGTTCAGCCTCGCTTTCCTCGCGGATTAGCACCTTGCGGAAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACaaaaaa                                                                                     <  1:1184699/150‑1 (MQ=255)
                                                                  gctgctCACGCGGGATCACGTTCACCTGTGGCGGGGCGACGGCGTTTTCAGCCTCGGTTTCCTCGGGGTTTACAACATGGCGGAAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACaaaaaa                                                                                     <  1:1640333/150‑1 (MQ=255)
                                                                  gctgctCACGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACaaaaaa                                                                                     >  1:1327298/1‑150 (MQ=255)
                                                                     gctCACGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATa                                                                                  >  1:2134541/1‑150 (MQ=255)
                                                                      ctCACGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAg                                                                                 >  2:561927/1‑150 (MQ=255)
                                                                       acccgggggATCACCGTCACCTGTGGAGGGGCGACTGCCTTTTCTGCCTCGGTTTCCGCGCGGATTTGCCCCTTGCGGAAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGc                                                                                <  1:893451/144‑1 (MQ=255)
                                                                         ccgcgGGATCCCCGTCACCTGTGGAGGGGGGACTGCCTTTTCAGCCTCGCTTTCCGCGCGGCTTACCACCTTGCGGAAAAAATCAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCgg                                                                              <  2:1155/149‑1 (MQ=255)
                                                                         aCGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCgg                                                                              >  2:1473204/1‑150 (MQ=255)
                                                                         aCGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCgg                                                                              >  1:739659/1‑150 (MQ=255)
                                                                             ggtaTCACCTTCACCTGTGGAGGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGATTAGCACCTTGCGGAAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTaa                                                                          <  2:1819852/147‑1 (MQ=255)
                                                                             ggtaTCACCGTCACATGTGGACGGGCGACTGCCTGTACAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGAAAAAATTCGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTaa                                                                          <  2:1819647/147‑1 (MQ=255)
                                                                             gggATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTaa                                                                          >  1:582157/1‑150 (MQ=255)
                                                                                 tCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCAt                                                                      >  2:594337/1‑150 (MQ=255)
                                                                                 tCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCAt                                                                      >  2:1042447/1‑150 (MQ=255)
                                                                                  cACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATa                                                                     >  2:2074539/1‑150 (MQ=255)
                                                                                   cccttCACGTGTGGTGGGGCGACTGACTGTTCAGACTCGCTTTCCTCGCGGTTTAGCACCTTGCGTAAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAg                                                                    <  2:1899515/146‑1 (MQ=255)
                                                                                   aCCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGATTAGCACCTTGCGGAAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAg                                                                    <  2:1899543/150‑1 (MQ=255)
                                                                                    ccGTCACCTGTGGAGGGGCTAGTGCCTTTACAGCCTCGTTTTCCTCGGGTGTTAGCACTTTGCGCAAAAAATACGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGc                                                                   <  1:498544/150‑1 (MQ=255)
                                                                                    ccGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGc                                                                   >  1:434464/1‑150 (MQ=255)
                                                                                            tgtgGCCGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGTTTCGCACCTTGCGGAAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATg                                                           <  2:202936/150‑1 (MQ=255)
                                                                                             gtgGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGa                                                          >  1:141125/1‑150 (MQ=255)
                                                                                              tgGACGGGCGACTGCCTTTTCAGCCGCGCTTTCCTCGCGGCTTAGCACCTTGCGGAAAAATTTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGAc                                                         <  2:254387/150‑1 (MQ=255)
                                                                                              tgGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGAAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGAc                                                         <  2:94686/150‑1 (MQ=255)
                                                                                              tgGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGATTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGAc                                                         <  2:94696/150‑1 (MQ=255)
                                                                                                   gggCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCAtt                                                    >  2:1996765/1‑150 (MQ=255)
                                                                                                    ggCGACTGCCTGTTCAGCCTCGATTTCCTCGCGGCTTAGCACCTTGCGGAAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCAttt                                                   <  1:1145306/150‑1 (MQ=255)
                                                                                                    ggCGACTGCCTGTCCAGCCTCGCTTACCTCGCGGATTAGCACCTTGCGGAAAAAATCCGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCAttt                                                   <  1:1145264/150‑1 (MQ=255)
                                                                                                    ggCGACTGCCTGTACCGCCTCGCTTTCCTCGCGGATTACAACCTGGCGGAAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCAttt                                                   <  1:452583/150‑1 (MQ=255)
                                                                                                     gCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTg                                                  <  2:2068512/150‑1 (MQ=255)
                                                                                                          tGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAAt                                             >  2:723384/1‑150 (MQ=255)
                                                                                                          tGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGAAAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAAt                                             >  1:1930053/1‑150 (MQ=255)
                                                                                                          gggttGTTCAGCCTCGCTTTCTAGGCTGCTTTGAAACTGGCGTCAAAAATTAGCGAATCGGGTAAAAATAGTACCCCTCGGTAGTGTAAAACATCATTCAGGACAAAAAAATAGCGGCTAATCCTTGCTCAGCATGACGCATTTGAGAAt                                             <  2:1598251/146‑1 (MQ=255)
                                                                                                           gCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATg                                            >  2:229550/1‑150 (MQ=255)
                                                                                                           gCCAGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCCAAGCTCAGCATGACGCAGGTGAGaaag                                            >  2:230167/1‑148 (MQ=255)
                                                                                                             cTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGtt                                          >  1:1595900/1‑150 (MQ=255)
                                                                                                               gTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGa                                        >  1:129248/1‑150 (MQ=255)
                                                                                                               gTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGa                                        >  2:1999720/1‑150 (MQ=255)
                                                                                                                ttCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGaa                                       >  2:505909/1‑150 (MQ=255)
                                                                                                                 tCAGCCTCGCTTTCCTAGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAAt                                      <  2:1130339/150‑1 (MQ=255)
                                                                                                                   aGCCTCGCTTTCCTCGCGGCTTAGAACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAAttt                                    >  1:1652441/1‑150 (MQ=255)
                                                                                                                    gCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGAAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTa                                   <  2:1991391/150‑1 (MQ=255)
                                                                                                                        cGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAAt                               <  1:65912/150‑1 (MQ=255)
                                                                                                                          cTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTg                             >  1:310437/1‑150 (MQ=255)
                                                                                                                           tttCCTCGCGGATTAGCACCTTGCGCAAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGc                            <  1:445372/150‑1 (MQ=255)
                                                                                                                             tCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCg                          >  2:1529102/1‑150 (MQ=255)
                                                                                                                             tCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCg                          >  2:1529315/1‑150 (MQ=255)
                                                                                                                               cTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGAc                        >  2:331219/1‑150 (MQ=255)
                                                                                                                                tCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACt                       >  1:1620950/1‑150 (MQ=255)
                                                                                                                                tCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACt                       >  1:439112/1‑150 (MQ=255)
                                                                                                                                  gcgGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCg                     >  1:496363/1‑150 (MQ=255)
                                                                                                                                   gggattAGCACCTTGCGGAAAAAATTCGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCgg                    <  1:70591/146‑1 (MQ=255)
                                                                                                                                   cgGCTTAGCACCTTGCGGAAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCgg                    <  2:1166667/150‑1 (MQ=255)
                                                                                                                                    ggCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCggg                   <  2:1873226/150‑1 (MQ=255)
                                                                                                                                      ctaaccACCTTGCGGAAAAAATACGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCa                 <  1:552178/145‑1 (MQ=255)
                                                                                                                                         agcagcTGGCGTAAAAAATTCGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATCGCTCATCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCAcgg              <  1:122686/150‑1 (MQ=255)
                                                                                                                                         aGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCAcgg              >  2:1831671/1‑150 (MQ=255)
                                                                                                                                            aCCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCAcggcgg           >  2:963084/1‑150 (MQ=255)
                                                                                                                                               ttGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGGTaa        <  1:1529102/150‑1 (MQ=255)
                                                                                                                                               ttGCGGAAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGGTaa        <  1:1529315/150‑1 (MQ=255)
                                                                                                                                                   ggCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGGTAAgccg    >  1:1509316/1‑150 (MQ=255)
                                                                                                                                                    gCAAAAATTCGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGGTAAgccgc   <  2:1561070/150‑1 (MQ=255)
                                                                                                                                                    gCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGGTAAgccgc   <  2:1561056/150‑1 (MQ=255)
                                                                                                                                                    gCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGGTAAgccgc   >  1:1634561/1‑150 (MQ=255)
                                                                                                                                                    gCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGGTAAgccgc   >  1:1264809/1‑150 (MQ=255)
                                                                                                                                                     cAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGGTAAgccgcc  >  2:1000660/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAGCATGCTGCTCACGCGGGATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTCGGTAGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGGTAAGCCGCC  >  NC_000913/158952‑159250

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