Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913312,1690TG44.9% 92.3 / 28.0 78L20V (TTG→GTG) ykgLuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (15/28);  new base G (35/0);  total (50/28)
Fisher's exact test for biased strand distribution p-value = 1.77e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGATCAGACTGTTCCTGTGGCGTAAACGTTTTTATGTTAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAAC  >  NC_000913/312026‑312317
                                                                                                                                               |                                                                                                                                                    
ggATCAGACTGTTCCTGTGGCGTAAACGTTTTTATGTTAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATg                                                                                                                                                 >  1:464677/1‑149 (MQ=255)
ggATCAGACTGTTCCTGTGGCGTAAACGTTTTTATGTTAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGcc                                                                                                                                               <  2:601896/151‑1 (MQ=255)
  aTCAGACTGTTCCTGTGGCGTAAACGTTTTTATGTTAGTAACAACAGATTTACCACAGTTATCAACTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGGCAg                                                                                                                                             <  2:906719/151‑1 (MQ=255)
    cAGACTGTTCCTGTGGCGTAAACGTTTTTATGTTAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGc                                                                                                                                            <  1:192108/150‑1 (MQ=255)
    cAGACTGTTCCTGTGGCGTAAACGTTTTTATGTTAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGGTGATGCCAGc                                                                                                                                            >  2:867321/1‑150 (MQ=255)
     aGACTGTTCCTGTGGCGTAAACGTTTTTATGTTAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTCAAAGAGGGTTGATGCCCGCAt                                                                                                                                          >  2:201640/1‑151 (MQ=255)
          gTTCCTGTGGCGTAAACGTTTTTATGTTAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGGTGATGCCAGCATACAgg                                                                                                                                     >  2:305948/1‑151 (MQ=255)
          gTTCCTGTGGCGTAAACGTTTTTATGTTAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGGTGATGCCAGCATACAgg                                                                                                                                     >  1:147865/1‑151 (MQ=255)
           ttCCTGTGGCGTAAACGTTTTTATGTTAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAggg                                                                                                                                    >  1:100847/1‑151 (MQ=255)
             ccTGTGGCGTAAACGTTTTTATGTTAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGc                                                                                                                                  <  2:1086040/151‑1 (MQ=255)
               tgtgGCGTAAACGTTTTTATGTTAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGGTGATGCCAGCATACAgggggg                                                                                                                                 >  2:388644/1‑148 (MQ=255)
                gtgGCGTAAACGTTTTTATGTTAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGa                                                                                                                                <  2:922720/150‑1 (MQ=255)
                  ggCGTAAACGTTTTTATGTTAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGGTGATGCCAGCATACAGGGGCGaa                                                                                                                               >  2:185106/1‑149 (MQ=255)
                    cGTAAACGTTTTTATGTTAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTg                                                                                                                           >  2:487845/1‑151 (MQ=255)
                     gTAAACGTTTTTATGTTAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGGTGATGCCAGCATACAGGGGCGAAGCTgg                                                                                                                          >  2:372455/1‑151 (MQ=255)
                       aaaCGTTTTTATGTTAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGGTGATGCCAGCATACAGGGGCGAAGCTgg                                                                                                                          >  1:700096/1‑149 (MQ=255)
                         aCGTTTTTATGTTAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGGTGATGCCAGCATACAGGGGCGAAGCTGGTca                                                                                                                       >  1:1194775/1‑150 (MQ=255)
                           gTTTTTATGTTAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGGTGATGGCAGCCTACAAGGGCGAAGCTGGTcaa                                                                                                                      >  1:1190154/1‑149 (MQ=255)
                             ttttATGTTAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCCGTCTAAAAGAGGGGTGATGCCGGCATACAGGGGCGAAGCTGGTcaaca                                                                                                                    >  1:313537/1‑149 (MQ=255)
                             ttttATGTTAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTcaacaa                                                                                                                   <  1:159318/150‑1 (MQ=255)
                             ttttATGTTAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTcaaca                                                                                                                    >  1:1004893/1‑149 (MQ=255)
                             ttttATGTTAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGGTGATGCCAGGATACAGGGGCGAAGCTGGTcaaca                                                                                                                    >  2:983514/1‑149 (MQ=255)
                                 aTGTTAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCa                                                                                                               <  1:691495/150‑1 (MQ=255)
                                   gTTAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGGTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAAc                                                                                                             >  1:1149632/1‑150 (MQ=255)
                                    ttAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCCGCATACAGGGGCGAAGCTGGTCAACAAGTCAAc                                                                                                             >  1:1051176/1‑149 (MQ=255)
                                           cAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATa                                                                                                    <  2:1109048/151‑1 (MQ=255)
                                              cAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATg                                                                                                 <  2:70181/151‑1 (MQ=255)
                                              cAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGGTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATaa                                                                                                   >  2:79549/1‑149 (MQ=255)
                                               aGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTTAGTCTAAAAGAGGGGTGATGCCAGCATACAGGGGCGGAGCTGGTCAACAAGACAACATTAAGATAATgg                                                                                                >  2:915107/1‑151 (MQ=255)
                                               aGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGGTCAGGCTAACAAAGGGTTTGTTCCCGCCTTCAGGGGCGAAGGTGGCCAAAAAGGCAACATTAAGATaaag                                                                                                 >  2:586196/1‑148 (MQ=255)
                                                 aTTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGa                                                                                               <  2:421954/150‑1 (MQ=255)
                                                    tACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGGCTAAAAGAGGGGTGATGCCAGCATACAGGGGCGAAGCTGGGCAACAAGTCAACATTAAGATAATGGAatat                                                                                           >  1:1170504/1‑151 (MQ=255)
                                                      ccACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCATCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATAtt                                                                                          <  1:911990/150‑1 (MQ=255)
                                                      ccACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGATCAAAAGAGGGGTGATGCCAGGATACAGGGGGGAAGCTGGTCCACAAATCAAAATTAAAATAATGGAATAttt                                                                                         >  2:1081139/1‑150 (MQ=255)
                                                       cacaTTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTTAGTCTAAAAAAGGGGTGATGCCAGCATACAGGGGCGGAGCTGGTCAACAAGTCAACATTAAGATAATGGAATAtt                                                                                          >  2:678743/1‑149 (MQ=255)
                                                          aTTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGGTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCgg                                                                                       >  1:367659/1‑149 (MQ=255)
                                                            ttATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTTAGTCTAAAAGAGGGGTGATGGCCAAAAAAAGGGGGGAAGCTGGTCAACAAATCCACCTTAAGAAAATGGAATATTCGgag                                                                                     >  2:290252/1‑149 (MQ=255)
                                                                catcTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGGTGATGCCAGCATACAGGGGGGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGaaa                                                                                >  2:612166/1‑150 (MQ=255)
                                                                    ttAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTTAGTCTAAAAGAGGGGTGATGCCAGCATACAGGGGGGGAGCTGGTCAACAAATCAACATTAAGATAATGGAAAATTCGGGGAGAAAAGTc                                                                            >  2:437140/1‑150 (MQ=255)
                                                                       aCAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGGTGGTGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAgg                                                                         >  2:462943/1‑150 (MQ=255)
                                                                       aCAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAAAGGGGTGATGCCAGCCTACAGGGGCGGAGCTTGTCAACAAGTCAACATTTAGATAATGGAATATTCGGAGAGAAACGTCAgg                                                                         >  1:349586/1‑150 (MQ=255)
                                                                        cAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTTAAAGAGGGGTGATGCCAGCATACAGGGGCGAAGGTGGGCAACAAGTCAACATTAAGATAAAGGAATATTCGGAGAGGAACGTCa                                                                           >  2:888686/1‑147 (MQ=255)
                                                                          gTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTTAATCTAAAAGAGGGGTGATGCCAGCATAAAGGGGCGAAGCTGGTCAAAAAATCAACATTAAGATAATGGAATATTCGGGGAGAAAAGACAggta                                                                       >  2:1016753/1‑147 (MQ=255)
                                                                            ttGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGGTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGtt                                                                    >  2:105174/1‑150 (MQ=255)
                                                                              gAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTgg                                                                  <  2:189968/150‑1 (MQ=255)
                                                                                gTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGcc                                                                <  2:711244/150‑1 (MQ=255)
                                                                                gTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTAAGTCTAAAAAAGGGGTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGcc                                                                >  2:519519/1‑150 (MQ=255)
                                                                                  gAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCa                                                             <  1:34955/151‑1 (MQ=255)
                                                                                   aaaaTGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGGTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCa                                                             >  1:1064300/1‑150 (MQ=255)
                                                                                        gCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGaa                                                       <  1:726496/151‑1 (MQ=255)
                                                                                           ggCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAAAGGGGTGATGCCAGCATACAGGGGCGAAGCTGGGCAAAAAGTCAACAATAAGATAATGGAAAATTCGGAGAGAAACGTCAGGCAGGTGGCCTCAAAAGaa                                                       >  1:436350/1‑148 (MQ=255)
                                                                                              tGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGGTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGga                                                  >  2:837604/1‑150 (MQ=255)
                                                                                                 gccgcAAAGGGAGGAGACTGAGTGGCGTGTTTAGTCTAAAAGAGGGGTGATGCCAGCATACAGGGGCGGAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAAAAGgaaga                                               >  2:476091/1‑150 (MQ=255)
                                                                                                   cgcAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAata                                            <  1:186302/151‑1 (MQ=255)
                                                                                                     cAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAatata                                          <  2:912320/151‑1 (MQ=255)
                                                                                                      aaaGGGAGGAGACTGAGTGGCGTGTTCAGTCTCAAAGAGGGGTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAAAAGGAAGAat                                             >  2:996022/1‑147 (MQ=255)
                                                                                                      aaaGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAatata                                          >  2:225356/1‑150 (MQ=255)
                                                                                                      aaaGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACCTCAGGCAGTTGGCCTCAAATGAACAGgaaga                                               >  1:415289/1‑145 (MQ=255)
                                                                                                      aaaGGGAGGAGACTGAGTGGCGTGTTCAGGCTAAAAGAGGGGTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAatata                                          >  1:623145/1‑150 (MQ=255)
                                                                                                       aaGGGAGGAGACTGAGTGGCGTGCTCAGTCTAAAAGAAGGGTGATGCCAGCATACAGGGGGGGAGCTGGTCAACAAAGCAACCTTAAAATAATGGGAAATTCGGAGAGAAACGTCAGGGAGTTGGCCCCCAATGAAAAGGAAGAatatat                                         >  2:669656/1‑150 (MQ=255)
                                                                                                         gggAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATccc                                      <  1:289313/151‑1 (MQ=255)
                                                                                                          ggaggaGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGGTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATcccc                                     >  2:354455/1‑151 (MQ=255)
                                                                                                          ggagcaGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGAATACAGGTGGGAAGCAGGTCAACAACTTAACATTAAAATAATGGAATATTCGGAGACAAACGTCAGGTAGTTGGCCTCAAATGATCTGGATTAATATATc                                        <  2:2275/148‑1 (MQ=255)
                                                                                                            aggagACTGAGTGGCGTGTTCAGGCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATcccc                                     >  1:154025/1‑149 (MQ=255)
                                                                                                             ggagACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAgg                                  <  1:932568/151‑1 (MQ=255)
                                                                                                               agaCTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGaa                                <  1:568004/151‑1 (MQ=255)
                                                                                                                gaCTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGaaa                               <  2:1051470/151‑1 (MQ=255)
                                                                                                                    gAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAAtc                            <  2:510860/150‑2 (MQ=255)
                                                                                                                      gTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAAt                         <  1:289150/151‑1 (MQ=255)
                                                                                                                       tGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATg                        >  2:360278/1‑151 (MQ=255)
                                                                                                                          cGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTg                     <  2:1190095/151‑1 (MQ=255)
                                                                                                                           gtgtTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTgg                    >  2:613121/1‑151 (MQ=255)
                                                                                                                                 aGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTAAACAAGACAACATTTAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCGCCAGGAAAATTAATGTTGGTGTtat              >  2:1168846/1‑151 (MQ=255)
                                                                                                                                      aaaaGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATaa            <  1:201645/148‑1 (MQ=255)
                                                                                                                                       aaaGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAacac        >  1:1171789/1‑151 (MQ=255)
                                                                                                                                       aaaGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAacac        >  2:295536/1‑151 (MQ=255)
                                                                                                                                             ggTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAAc  <  1:185108/151‑1 (MQ=255)
                                                                                                                                             ggTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAAc  <  2:747617/151‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GGATCAGACTGTTCCTGTGGCGTAAACGTTTTTATGTTAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAAC  >  NC_000913/312026‑312317

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