Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,042,6040AC39.4% 136.6 / 16.8 71V182G (GTG→GGG) yqfAhemolysin III family HylIII inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (33/10);  new base C (1/27);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 2.64e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGC  >  NC_000913/3042456‑3042749
                                                                                                                                                       |                                                                                                                                                 
gcgcGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCAcc                                                                                                                                                 >  2:276518/1‑151 (MQ=255)
 cgcgTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCAccc                                                                                                                                                >  1:246747/1‑151 (MQ=255)
 cgcgTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCAccc                                                                                                                                                >  1:40848/1‑151 (MQ=255)
   cgTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCg                                                                                                                                              <  1:535339/151‑1 (MQ=255)
    gTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGa                                                                                                                                             <  1:734388/151‑1 (MQ=255)
     tCTTATCAGGCCTACTAAGGCATACCCATTCCGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCCCCCCg                                                                                                                                              <  2:795230/149‑1 (MQ=255)
      cTTATCCGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGc                                                                                                                                           >  2:902950/1‑151 (MQ=255)
            aGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATaa                                                                                                                                     <  1:947299/151‑1 (MQ=255)
            aGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGAACGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATa                                                                                                                                      <  2:1006137/150‑1 (MQ=255)
             ggCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATaaa                                                                                                                                    <  1:305014/151‑1 (MQ=255)
             ggCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATaa                                                                                                                                     >  2:479837/1‑150 (MQ=255)
               ccTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAAcc                                                                                                                                  <  2:347387/151‑1 (MQ=255)
                 tACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCAc                                                                                                                                >  2:989546/1‑151 (MQ=255)
                  aCAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACGACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCAc                                                                                                                                <  2:1005224/150‑1 (MQ=255)
                   cAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCAcgc                                                                                                                              <  1:349711/151‑1 (MQ=255)
                     aaGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCAcgccg                                                                                                                            >  1:2474/1‑151 (MQ=255)
                               cATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCa                                                                                                                  <  1:502771/151‑1 (MQ=255)
                                aTTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGcc                                                                                                                   >  1:933535/1‑149 (MQ=255)
                                 ttACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGGTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGt                                                                                                                >  2:330253/1‑151 (MQ=255)
                                   aCGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTa                                                                                                               >  2:294025/1‑150 (MQ=255)
                                         gCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAAc                                                                                                        >  1:591146/1‑151 (MQ=255)
                                          cccAATATACAAATAGATCGCCAGAAAGGGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTaa                                                                                                         <  1:324050/149‑1 (MQ=255)
                                                ataCAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGccc                                                                                                 >  1:1051787/1‑151 (MQ=255)
                                                    aaaTAGACCGCAAAAAAGTGGCACCCACTCCCGCCGAGCACGAACCCGTGCCAAATGGCTTGTTTTTATGGAATGCGTTTGCAGCCGTAAAAAACCCCCCCGAGCGAATAACCCCCGCCCCCTACCGCCAGTAAGGTAACGCTGCccgccg                                                                                             <  2:362658/151‑1 (MQ=255)
                                                        aGACCGCAAAAAAGGGCCACCCCCTCCCCCCGGGCAGGAGGCGGTGCCAGAGGCCATGGTTGTAGGAAATGGGTTTGAAGAGGAAAAAAACCCCCCCGGGCAAATAACCCACGCCCCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAg                                                                                         <  1:56335/151‑1 (MQ=255)
                                                         gATCGCCAGAAAGCGGCACACACTACCGCCGAGCACGACGCGGTGCCAGATGGCATGGATGTATGGAATGCGTTTGCAGACGTAGAAAATCACCGCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCAGCCcgccgc                                                                                            >  2:45325/1‑147 (MQ=255)
                                                            cccaaaaaaGTGCCCCACCCTCCCGCCGAGCACGAAGCCGTGCCAGATGGCATGTTTGTATGGAATGCGTTTGCAGACGTAAAAAACCCCCCCGAGCAAATAAACCCCCCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCtt                                                                                      <  1:758247/144‑1 (MQ=255)
                                                              ccAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAAc                                                                                   >  2:1090303/1‑151 (MQ=255)
                                                              ccAAAAAGGGGCACACCCTCCCGCCGAGCACGAAGCCGTGCCAGATGGCATGTTTGTATGGAATGCGTTTGCAGACGTAAAAAATCCCCCCGAGCGAAAAAACCACGCCCCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTaa                                                                                    <  1:270009/150‑1 (MQ=255)
                                                                aGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACt                                                                                  >  2:586358/1‑150 (MQ=255)
                                                                aGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTg                                                                                 >  1:295172/1‑151 (MQ=255)
                                                                aGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTg                                                                                 >  2:849629/1‑151 (MQ=255)
                                                                  aaaGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGcc                                                                               >  2:849066/1‑151 (MQ=255)
                                                                  aaaGGGGCACACACTCCCGCGGAGCACGAAGCCGTGCCAGATGGCATGGTTTTATGAAATGCTTTTGCAGCCGTAAAAAACCCCCCCGGGCGAAAAACCCCCGCCGCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGc                                                                                <  1:877727/150‑1 (MQ=255)
                                                                  aaaGGGCCACACACTCCCGCCGAGCACGAAGCCGTGCCAGAGGGCATGGTTGTATGAAATGCTTTTGCAGACGTAGAAAACCCCCCCGAGCAAATAAACCCCGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAAc                                                                                   <  1:795633/147‑1 (MQ=255)
                                                                             cccccccgccgAGCAGGAACCCGTGCAAGATGGCATGGTTGTAGGGAATGCTTTTGCAGACGTAAAAAACCCCCCGGAGCAAAAAAACCCCGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATaaa                                                                     <  1:49399/145‑1 (MQ=255)
                                                                                  ccgccgAGCCAGAAGCCGTGCCAGATGGGATGGTTGAAGGGAAGGCGTTTTCAGGCGAAGAAAATCCCCCAGGGCGAAAAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCACCGCGAGCTTAACTGCCATTTCATAAATTacca                                                               <  2:546669/151‑1 (MQ=255)
                                                                                   ccgcaaGCACCAACCCGGGCCAGAATGGCATGTTTTTAAGGAAGGGTTTGACAATGTTAAAAAACCCCCCCGGGGAAATAACCACCCCCTCTTACCGCCAGTAAGTTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTacca                                                               <  2:556041/146‑1 (MQ=255)
                                                                                       gAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAg                                                           >  1:538990/1‑150 (MQ=255)
                                                                                             gAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGAcagc                                                    >  1:555074/1‑151 (MQ=255)
                                                                                                 ccGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGAcagccagc                                                >  1:106286/1‑151 (MQ=255)
                                                                                                   gTGCCGGATGGCAGGTTTTTAGGAAATGGTTTGGCAGCCGTAAAAAACCCCCCCGAGCAAAAAACCCCCGCCGCTTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGccc                                              <  1:989827/151‑1 (MQ=255)
                                                                                                   gTGCAAGAGGGCATGGTTGTAGGGAATGCTTTTGCAGACGAAAAAAACCCCCCCGAGCAAAAAACCCCCGCCCCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGccc                                              <  1:761771/151‑1 (MQ=255)
                                                                                                       caggagggCATGTTTGTATGGAATGCTTTTGCAGACGTAAAAAACCCCCCCGAGCAAATAACCCCCGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCa                                             <  1:368708/143‑1 (MQ=255)
                                                                                                       cAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCg                                          >  1:291732/1‑151 (MQ=255)
                                                                                                       cAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCg                                          >  2:919803/1‑151 (MQ=255)
                                                                                                        aGATGGCTTGGTTGTATGAAATGCGTTTGCAGACGTAAAAAACCCCCCCGAGCGAATAACCCCCGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCg                                          <  1:127306/150‑1 (MQ=255)
                                                                                                             ggagtGGTTTTGTGAAGTGCTTTTGCAGACGAAAAAAACCCCCCCGAGAAAAAAACCCCCCCCGCCGCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCg                                          <  2:772147/142‑1 (MQ=255)
                                                                                                             ggCATGTTTTTATGAAATGCTTTTGCAGACGTAGAAAATCCCCCCGAGCGAATAACCCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAg                                      <  1:765985/150‑1 (MQ=255)
                                                                                                                aTGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGAt                                    >  2:179396/1‑149 (MQ=255)
                                                                                                                  gtttttATGGAATGGTTTTGCAGACGTAAAAAATCCCCCCGAGCAAATAACCCCCGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGt                                <  1:912897/146‑1 (MQ=255)
                                                                                                                  gtttGTAGGGAATGCGTTTGCAGACGAAAAAAACCCCCCCGAGCGAAAAACCCCCCCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATa                                   <  1:809606/146‑1 (MQ=255)
                                                                                                                  ggTTGTATGGAATGGTTTTGCAGCCGAAAAAAACCCCCCCGAGCGAATAACCCCCCCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAgg                                 <  1:170044/150‑1 (MQ=255)
                                                                                                                    tgttATGGAATCGTTTGGCAGACGAAAAAATCCCCCCCGGGCAAAAAACCCCCCCCCCCACCCGCCATTAAGGTACCCCGGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATa                                   <  2:601894/143‑1 (MQ=255)
                                                                                                                       ttgggAAGGCTTTTGCCGACGTAAAAACCCCCCCCGAGCAAAAAACCCCCCCCCCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCAc                             <  1:35354/146‑1 (MQ=255)
                                                                                                                       taggaaaTGCGTTTGCAGACGTAAAAAACCCCCCCGAGCAAAAAACCCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCAcc                            <  1:453522/145‑1 (MQ=255)
                                                                                                                          ggAATGCGTTTGCAGACGTAAAAAACCCCCCCGAGCAAAAAACCCCCGCCGCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCaga                         <  1:668207/150‑1 (MQ=255)
                                                                                                                                cttttGCAGCGCTAAAAAATCCCCCCGGGCGAAGAAACCCCCCCCCCTACCGCCAGTAAGGTAACGCTCCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATaa                    <  2:957162/147‑1 (MQ=255)
                                                                                                                                  tttGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAAtt                >  1:765325/1‑151 (MQ=255)
                                                                                                                                   ttGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAAttt               >  1:345395/1‑151 (MQ=255)
                                                                                                                                    ttggcggACGTAAAAAACCCCCCGGAGGAAATAACCCCCCCCCCCGCCCGCCATTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAAt                 <  2:1050018/144‑1 (MQ=255)
                                                                                                                                         aGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAACTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTg             >  2:208412/1‑149 (MQ=255)
                                                                                                                                         aGACGTAGAAAACCCCCCCGAGCGAATAAACCACGCCGCCTACCCCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCCGATAGCTCACCAGGGATACAATTTTGa            >  2:405556/1‑150 (MQ=255)
                                                                                                                                           aggTAAAAAATCCCCCGGAGCGAAAAAACCACCCCGCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAAt          <  2:533929/148‑1 (MQ=255)
                                                                                                                                           aCGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATc         >  1:913608/1‑151 (MQ=255)
                                                                                                                                             gTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCgg       >  1:1048608/1‑151 (MQ=255)
                                                                                                                                             gTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCgg       >  2:306491/1‑151 (MQ=255)
                                                                                                                                                 aaaaTCCCCCCGAGCAAATAACCCCCCCCCCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAAtttt              <  2:337047/140‑1 (MQ=255)
                                                                                                                                                 aaaaTCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTgcg   >  2:3990/1‑151 (MQ=255)
                                                                                                                                                  aaaTCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTgcgc  >  1:849471/1‑151 (MQ=255)
                                                                                                                                                  aaaTCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTgcg   >  1:272445/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                 
GCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGC  >  NC_000913/3042456‑3042749

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