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

GCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAG‑ATGGCATGGTTGTATGGAATGCGTTT‑G‑CAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGC  >  NC_000913/3042456‑3042749
                                                                                                                                                       |                                                                                                                                                 
gcgcGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAG‑ATGGCATGGTTGTATGGAATGCGTTT‑G‑CAGACGTAGAAAATCAcc                                                                                                                                                 >  2:276518/1‑151 (MQ=255)
 cgcgTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAG‑ATGGCATGGTTGTATGGAATGCGTTT‑G‑CAGACGTAGAAAATCAccc                                                                                                                                                >  1:246747/1‑151 (MQ=255)
 cgcgTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAG‑ATGGCATGGTTGTATGGAATGCGTTT‑G‑CAGACGTAGAAAATCAccc                                                                                                                                                >  1:40848/1‑151 (MQ=255)
   cgTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAG‑ATGGCATGGTTGTATGGAATGCGTTT‑G‑CAGACGTAGAAAATCACCCCg                                                                                                                                              <  1:535339/151‑1 (MQ=255)
    gTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAG‑ATGGCATGGTTGTATGGAATGCGTTT‑G‑CAGACGTAGAAAATCACCCCGa                                                                                                                                             <  1:734388/151‑1 (MQ=255)
     tCTTATCAGGCCTACTAAGGCATACCCATTCCGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAG‑ATGGCATGGTTGTATGGAATGCGTTT‑G‑CAGACGTAGAAAATCCCCCCg                                                                                                                                              <  2:795230/149‑1 (MQ=255)
      cTTATCCGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAG‑ATGGCATGGTTGTATGGAATGCGTTT‑G‑CAGACGTAGAAAATCACCCCGAGc                                                                                                                                           >  2:902950/1‑151 (MQ=255)
            aGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAG‑ATGGCATGGTTGTATGGAATGCGTTT‑G‑CAGACGTAGAAAATCACCCCGAGCGAATaa                                                                                                                                     <  1:947299/151‑1 (MQ=255)
            aGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGAACGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAG‑ATGGCATGGTTGTATGGAATGCGTTT‑G‑CAGACGTAGAAAATCACCCCGAGCGAATa                                                                                                                                      <  2:1006137/150‑1 (MQ=255)
             ggCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAG‑ATGGCATGGTTGTATGGAATGCGTTT‑G‑CAGACGTAGAAAATCACCCCGAGCGAATaaa                                                                                                                                    <  1:305014/151‑1 (MQ=255)
             ggCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAG‑ATGGCATGGTTGTATGGAATGCGTTT‑G‑CAGACGTAGAAAATCACCCCGAGCGAATaa                                                                                                                                     >  2:479837/1‑150 (MQ=255)
               ccTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAG‑ATGGCATGGTTGTATGGAATGCGTTT‑G‑CAGACGTAGAAAATCACCCCGAGCGAATAAAcc                                                                                                                                  <  2:347387/151‑1 (MQ=255)
                 tACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAG‑ATGGCATGGTTGTATGGAATGCGTTT‑G‑CAGACGTAGAAAATCACCCCGAGCGAATAAACCAc                                                                                                                                >  2:989546/1‑151 (MQ=255)
                  aCAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACGACCGCCGAGCACGAAGCCGTGCCAG‑ATGGCATGGTTGTATGGAATGCGTTT‑G‑CAGACGTAGAAAATCACCCCGAGCGAATAAACCAc                                                                                                                                <  2:1005224/150‑1 (MQ=255)
                   cAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAG‑ATGGCATGGTTGTATGGAATGCGTTT‑G‑CAGACGTAGAAAATCACCCCGAGCGAATAAACCAcgc                                                                                                                              <  1:349711/151‑1 (MQ=255)
                     aaGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAG‑ATGGCATGGTTGTATGGAATGCGTTT‑G‑CAGACGTAGAAAATCACCCCGAGCGAATAAACCAcgccg                                                                                                                            >  1:2474/1‑151 (MQ=255)
                               cATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAG‑ATGGCATGGTTGTATGGAATGCGTTT‑G‑CAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCa                                                                                                                  <  1:502771/151‑1 (MQ=255)
                                aTTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAG‑ATGGCATGGTTGTATGGAATGCGTTT‑G‑CAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGcc                                                                                                                   >  1:933535/1‑149 (MQ=255)
                                 ttACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAG‑ATGGCATGGTTGTATGGAATGCGGTT‑G‑CAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGt                                                                                                                >  2:330253/1‑151 (MQ=255)
                                   aCGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAG‑ATGGCATGGTTGTATGGAATGCGTTT‑G‑CAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTa                                                                                                               >  2:294025/1‑150 (MQ=255)
                                         gCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAG‑ATGGCATGGTTGTATGGAATGCGTTT‑G‑CAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAAc                                                                                                        >  1:591146/1‑151 (MQ=255)
                                          cccAATATACAAATAGATCGCCAGAAAGGGGCACACACTACCGCCGAGCACGAAGCCGTGCCAG‑ATGGCATGGTTGTATGGAATGCGTTT‑G‑CAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTaa                                                                                                         <  1:324050/149‑1 (MQ=255)
                                                ataCAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAG‑ATGGCATGGTTGTATGGAATGCGTTT‑G‑CAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGccc                                                                                                 >  1:1051787/1‑151 (MQ=255)
                                                    aaaTAGACCGCAAAAAAGTGGCACCCACTCCCGCCGAGCACGAACCCGTGCCAA‑ATGGCTTGTTTTTATGGAATGCGTTT‑G‑CAGCCGTAAAAAACCCCCCCGAGCGAATAACCCCCGCCCCCTACCGCCAGTAAGGTAACGCTGCccgccg                                                                                             <  2:362658/151‑1 (MQ=255)
                                                        aGACCGCAAAAAAGGGCCACCCCCTCCCCCCGGGCAGGAGGCGGTGCCAG‑AGGCCATGGTTGTAGGAAATGGGTTT‑G‑AAGAGGAAAAAAACCCCCCCGGGCAAATAACCCACGCCCCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAg                                                                                         <  1:56335/151‑1 (MQ=255)
                                                         gATCGCCAGAAAGCGGCACACACTACCGCCGAGCACGACGCGGTGCCAG‑ATGGCATGGATGTATGGAATGCGTTT‑G‑CAGACGTAGAAAATCACCGCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCAGCCcgccgc                                                                                            >  2:45325/1‑147 (MQ=255)
                                                            cccaaaaaaGTGCCCCACCCTCCCGCCGAGCACGAAGCCGTGCCAG‑ATGGCATGTTTGTATGGAATGCGTTT‑G‑CAGACGTAAAAAACCCCCCCGAGCAAATAAACCCCCCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCtt                                                                                      <  1:758247/144‑1 (MQ=255)
                                                              ccAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAG‑ATGGCATGGTTGTATGGAATGCGTTT‑G‑CAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAAc                                                                                   >  2:1090303/1‑151 (MQ=255)
                                                              ccAAAAAGGGGCACACCCTCCCGCCGAGCACGAAGCCGTGCCAG‑ATGGCATGTTTGTATGGAATGCGTTT‑G‑CAGACGTAAAAAATCCCCCCGAGCGAAAAAACCACGCCCCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTaa                                                                                    <  1:270009/150‑1 (MQ=255)
                                                                aGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAG‑ATGGCATGGTTGTATGGAATGCGTTT‑G‑CAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACt                                                                                  >  2:586358/1‑150 (MQ=255)
                                                                aGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAG‑ATGGCATGGTTGTATGGAATGCGTTT‑G‑CAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTg                                                                                 >  1:295172/1‑151 (MQ=255)
                                                                aGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAG‑ATGGCATGGTTGTATGGAATGCGTTT‑G‑CAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTg                                                                                 >  2:849629/1‑151 (MQ=255)
                                                                  aaaGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAG‑ATGGCATGGTTGTATGGAATGCGTTT‑G‑CAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGcc                                                                               >  2:849066/1‑151 (MQ=255)
                                                                  aaaGGGGCACACACTCCCGCGGAGCACGAAGCCGTGCCAG‑ATGGCATGGTTTTATGAAATGCTTTT‑G‑CAGCCGTAAAAAACCCCCCCGGGCGAAAAACCCCCGCCGCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGc                                                                                <  1:877727/150‑1 (MQ=255)
                                                                  aaaGGGCCACACACTCCCGCCGAGCACGAAGCCGTGCCAG‑AGGGCATGGTTGTATGAAATGCTTTT‑G‑CAGACGTAGAAAACCCCCCCGAGCAAATAAACCCCGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAAc                                                                                   <  1:795633/147‑1 (MQ=255)
                                                                             cccccccgccgAGCAGGAACCCGTGCAAG‑ATGGCATGGTTGTAGGGAATGCTTTT‑G‑CAGACGTAAAAAACCCCCCGGAGCAAAAAAACCCCGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATaaa                                                                     <  1:49399/145‑1 (MQ=255)
                                                                                  ccgccgAGCCAGAAGCCGTGCCAG‑ATGGGATGGTTGAAGGGAAGGCGTTT‑T‑CAGGCGAAGAAAATCCCCCAGGGCGAAAAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCACCGCGAGCTTAACTGCCATTTCATAAATTacca                                                               <  2:546669/151‑1 (MQ=255)
                                                                                   ccgcaaGCACCAACCCGGGCCAGAATGGCATGTTTTTAAGGAA‑GGGTTT‑GACAATGTTAAAAAACCCCCCCGGGGAAATAACCACCCCCTCTTACCGCCAGTAAGTTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTacca                                                               <  2:556041/146‑1 (MQ=255)
                                                                                       gAGCACGAAGCCGTGCCAG‑ATGGCATGGTTGTATGGAATGCGTTT‑G‑CAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAg                                                           >  1:538990/1‑150 (MQ=255)
                                                                                             gAAGCCGTGCCAG‑ATGGCATGGTTGTATGGAATGCGTTT‑G‑CAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGAcagc                                                    >  1:555074/1‑151 (MQ=255)
                                                                                                 ccGTGCCAG‑ATGGCATGGTTGTATGGAATGCGTTT‑G‑CAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGAcagccagc                                                >  1:106286/1‑151 (MQ=255)
                                                                                                   gTGCCGG‑ATGGCAGGTTTTTAGGAAATGGTTTG‑G‑CAGCCGTAAAAAACCCCCCCGAGCAAAAAACCCCCGCCGCTTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGccc                                              <  1:989827/151‑1 (MQ=255)
                                                                                                   gTGCAAG‑AGGGCATGGTTGTAGGGAATGCTTTT‑G‑CAGACGAAAAAAACCCCCCCGAGCAAAAAACCCCCGCCCCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGccc                                              <  1:761771/151‑1 (MQ=255)
                                                                                                       caggagggCATGTTTGTATGGAATGCTTTT‑G‑CAGACGTAAAAAACCCCCCCGAGCAAATAACCCCCGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCa                                             <  1:368708/143‑1 (MQ=255)
                                                                                                       cAG‑ATGGCATGGTTGTATGGAATGCGTTT‑G‑CAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCg                                          >  1:291732/1‑151 (MQ=255)
                                                                                                       cAG‑ATGGCATGGTTGTATGGAATGCGTTT‑G‑CAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCg                                          >  2:919803/1‑151 (MQ=255)
                                                                                                        aG‑ATGGCTTGGTTGTATGAAATGCGTTT‑G‑CAGACGTAAAAAACCCCCCCGAGCGAATAACCCCCGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCg                                          <  1:127306/150‑1 (MQ=255)
                                                                                                             ggagtGGTTTTGTGAAGTGCTTTT‑G‑CAGACGAAAAAAACCCCCCCGAGAAAAAAACCCCCCCCGCCGCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCg                                          <  2:772147/142‑1 (MQ=255)
                                                                                                             ggCATGTTTTTATGAAATGCTTTT‑G‑CAGACGTAGAAAATCCCCCCGAGCGAATAACCCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAg                                      <  1:765985/150‑1 (MQ=255)
                                                                                                                aTGGTTGTATGGAATGCGTTT‑G‑CAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGAt                                    >  2:179396/1‑149 (MQ=255)
                                                                                                                  gtttttATGGAATGGTTTT‑G‑CAGACGTAAAAAATCCCCCCGAGCAAATAACCCCCGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGt                                <  1:912897/146‑1 (MQ=255)
                                                                                                                  gtttGTAGGGAATGCGTTT‑G‑CAGACGAAAAAAACCCCCCCGAGCGAAAAACCCCCCCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATa                                   <  1:809606/146‑1 (MQ=255)
                                                                                                                  ggTTGTATGGAATGGTTTT‑G‑CAGCCGAAAAAAACCCCCCCGAGCGAATAACCCCCCCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAgg                                 <  1:170044/150‑1 (MQ=255)
                                                                                                                    tgttATGGAAT‑CGTTTGG‑CAGACGAAAAAATCCCCCCCGGGCAAAAAACCCCCCCCCCCACCCGCCATTAAGGTACCCCGGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATa                                   <  2:601894/143‑1 (MQ=255)
                                                                                                                       ttgggAAGGCTTTT‑G‑CCGACGTAAAAACCCCCCCCGAGCAAAAAACCCCCCCCCCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCAc                             <  1:35354/146‑1 (MQ=255)
                                                                                                                       taggaaaTGCGTTT‑G‑CAGACGTAAAAAACCCCCCCGAGCAAAAAACCCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCAcc                            <  1:453522/145‑1 (MQ=255)
                                                                                                                          ggAATGCGTTT‑G‑CAGACGTAAAAAACCCCCCCGAGCAAAAAACCCCCGCCGCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCaga                         <  1:668207/150‑1 (MQ=255)
                                                                                                                                ctttt‑G‑CAGCGCTAAAAAATCCCCCCGGGCGAAGAAACCCCCCCCCCTACCGCCAGTAAGGTAACGCTCCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATaa                    <  2:957162/147‑1 (MQ=255)
                                                                                                                                  ttt‑G‑CAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAAtt                >  1:765325/1‑151 (MQ=255)
                                                                                                                                   tt‑G‑CAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAAttt               >  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)
                                                                                                                                                       |                                                                                                                                                 
GCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAG‑ATGGCATGGTTGTATGGAATGCGTTT‑G‑CAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGC  >  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: ‑