Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,162,1200TG42.6% 108.0 / 16.1 62G219G (GGT→GGG) iceTputative citrate/iron‑citrate/zinc‑citrate efflux transporter
Reads supporting (aligned to +/- strand):  ref base T (6/29);  new base G (26/0);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 1.37e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTA‑TTAACCCTGGCGCTGGACGGCAGTAAAGGT‑ACA‑GGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTCGCTGGGCCTGGCGGG  >  NC_000913/2161970‑2162269
                                                                                                                                                       |                                                                                                                                                       
cAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTA‑TTAACCCTGGCGCTGGACGGCAGTAAAGGt                                                                                                                                                         <  2:707537/151‑1 (MQ=255)
 aaCATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTA‑TTAACCCTGGCGCTGGACGGCAGTAAAGGT‑a                                                                                                                                                       <  1:1176760/151‑1 (MQ=255)
            ggggATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTA‑TTAACCCTGGCGCTGGACGGCAGTAAAGGT‑ACA‑GGTTTATcg                                                                                                                                           <  2:887516/151‑1 (MQ=255)
             gggATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAAGTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTA‑TTAACCCTGGCGCTGGACGGCAGTAAAGGT‑ACA‑GGTTTATc                                                                                                                                            <  1:339972/149‑1 (MQ=255)
                     cGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTA‑TTTAACCTGGGGCCGGAAGGGAGGAAAGGG‑AAA‑GGGTTATtgcc                                                                                                                                         >  2:37074/1‑144 (MQ=255)
                         gCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTA‑TTAACCCTGGCGCTGGACGGCAGTAAAGGT‑ACA‑GGTTTATCGCCGCTGACGATTg                                                                                                                              <  2:147391/151‑1 (MQ=255)
                            aTCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGGA‑TTAAACCTGGGGCTGGGCGGGAGTAAAGGG‑ACA‑AGTTTTTCGCCGCTGGCGGtt                                                                                                                               >  1:102442/1‑147 (MQ=255)
                               gCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGGA‑TTTAACCCCGCGCTGGGGGGCAGTTAAGGG‑TCA‑GGGTTATCGCCGCCGGCCATTGCCGGcc                                                                                                                        >  1:530591/1‑151 (MQ=255)
                                  acaTTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGGA‑TTTAACCTGGGGCTGGGCGGCAGTAAAAGG‑ACA‑GGGTTTTTGCCGCTGGCGAtttccgggccgg                                                                                                                     >  2:337281/1‑141 (MQ=255)
                                  acaTTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGGA‑TTAACCCTGGCGGTGGGAGGCAGGAAAGGG‑ACA‑GGGTTTTCGCCGCCGACGATTGCCGGGCTgg                                                                                                                     >  1:741200/1‑151 (MQ=255)
                                   caTTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCCTGGCGGGA‑TTAAACCTGGGGCTGGGAGGCAGGAAAAGG‑GTA‑GGGTTTTCGCCGCTGACGATTTCGGGCCCggg                                                                                                                    >  1:290224/1‑150 (MQ=255)
                                   caTTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTA‑TTAACCCTGGCGCTGGACGGCAGTAAAGGT‑ACA‑GGTTTATCGCCGCTGACGATTGCAGGCCTgg                                                                                                                     <  2:290241/150‑1 (MQ=255)
                                           ttAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGGA‑TTTAACCCGGCGCTGGGCGGCAGGAAAGGG‑ACA‑GGGTTTTCGGCGCCGAACGTTTCAGGCCCGGGCGGAggt                                                                                                             >  1:864478/1‑148 (MQ=255)
                                                gCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTA‑TTAACCCTGGCGCTGGACGGCAGTAAAGGT‑ACA‑GGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCgtg                                                                                                       <  1:7960/151‑1 (MQ=255)
                                                gCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTA‑TTAACCCTGGCGCTGGACGGCAGTAAAGGT‑ACA‑GGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCgt                                                                                                        <  1:730880/150‑1 (MQ=255)
                                                  cGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTA‑TTAACCCTGGCGCTGGACGGCAGTAAAGGT‑ACA‑GGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCgtggt                                                                                                     <  2:311736/151‑1 (MQ=255)
                                                  cGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTA‑TTAACCCTGGCGCTGGACGGCAGTAAAGGT‑ACA‑GGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCgtggt                                                                                                     <  2:490093/151‑1 (MQ=255)
                                                        acacCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTA‑TTAACCCTGGCGCTGGACGGCAGTAAAGGT‑ACA‑GGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCa                                                                                                 <  2:807110/149‑1 (MQ=255)
                                                          acCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTA‑TTTACCCTGGCGCTCGGCGGCAGTAAAGGG‑ACAGGGTTTTTCG‑CGCTGACGGTTGCAGGGCCGGTCGCAGGTGGCGTGGGGGCCCTg                                                                                              >  1:893814/1‑150 (MQ=255)
                                                             aTGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTA‑TTAACCCTGGCGCTGGACGGCAGTAAAGGT‑ACA‑GGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGc                                                                                          <  2:1153555/151‑1 (MQ=255)
                                                                 aGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGGA‑TTAAACCTGGGGCTGGACGGGAGGAAAGGG‑ACA‑GGGTTATTGCCGCTGGCGATTGCAGGGCTGGGCGCAGTTGGCGGGGGGGCAcctggggttt                                                                                       >  2:1014025/1‑141 (MQ=255)
                                                                    cgcgGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTA‑TTAACCCTGGCGCTGGACGGCAGTAAAGGT‑ACA‑GGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATct                                                                                   <  1:1279361/151‑1 (MQ=255)
                                                                    cgcgGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTA‑TTAACCCTGGCGCTGGACGGCAGTAAAGGT‑ACA‑GGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATc                                                                                    <  1:1130057/150‑1 (MQ=255)
                                                                    cgcgGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTA‑TTAACCCCGGCGCTTGGCGGCAGTAAAGGG‑ACA‑GGGTTTTTGCCGCTGACGATTGCCGGCCTGGGCGCAGGTGGCGGGGGGGCCCTGGGGCttttt                                                                                     >  1:800388/1‑147 (MQ=255)
                                                                    cgcgGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGGA‑TTAACCCTGGGGCTGGGCGGCAGGAAAGGG‑ACA‑GGGTTTTCGCCGCTGACGATTGCAGGCCTGGTCGCCGGTGGGGGGGGGGCACCGGTGCttt                                                                                       >  2:902779/1‑147 (MQ=255)
                                                                    cgcgGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGGA‑TTAAACCTTGCGCCGGGCGGCAGTAAAGGG‑ACA‑GGGTTATCGCCGCTGGCGATTGCAGGCCTGGTCGCAGGTGGCGTGGGGGGACTGGGGCTTTATct                                                                                   >  2:380600/1‑151 (MQ=255)
                                                                    cgcgGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGCTGGCATGGCGGTA‑TTAACCCTGGCGCTGGACGGCAGTAAAGGG‑ACA‑GGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGGTGGCGGGGTGGCACTGGTGCTTTAt                                                                                     >  1:358700/1‑149 (MQ=255)
                                                                     gcggcgCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTA‑TTAACCCTGGCGCTGGACGGCAGTAAAGGT‑ACA‑GGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATctg                                                                                  <  2:1072283/151‑1 (MQ=255)
                                                                       ggcgCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTA‑TTAACCCTGGCGCTGGACGGCAGTAAAGGT‑ACA‑GGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATctgct                                                                                <  1:127353/151‑1 (MQ=255)
                                                                        gcgcTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGGA‑TTAAACCTCGGGCTGGGCCGCAGAAAAAGA‑ACC‑GGGTTATTGCCCCTGACGGTTGCAGGCCTGGTGGGAGGGGGCGGGGGGGCACTgggggcttatttact                                                                                >  1:492196/1‑137 (MQ=255)
                                                                           cTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTA‑TTTAACCTGGGGCCGGGAGGCAGTAAAGGG‑ACA‑GGGTTTTCGCCGCTGACGGATGGAGGGCTGGTGGCAGTTGGCGTGGGGGCACTGGGGCTTTATctgctg                                                                               >  1:223865/1‑148 (MQ=255)
                                                                               gATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTA‑TTAACCCTGGCGCTGGACGGCAGTAAAGGT‑ACA‑GGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGcc                                                                         <  1:983210/150‑1 (MQ=255)
                                                                                  ctctcCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTA‑TTAACCCTGGCGCTGGACGGGAGGAAAGGG‑ACA‑GGGTTTTTGCCCCTGACGATTTCAGGCCTGGGCGCAGGTGGCGGGGGGGCACTGGGGCTTTTTTTGCTGCACGcc                                                                         >  1:411454/1‑147 (MQ=255)
                                                                                   tctcCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTTAACCCT‑GCGCCGGGCGGGAGTAAAAGG‑ACA‑GGGTTTTTGCCGCTGACGGTTTCCGGGCCGGGCGCCGGTTGGGTGGGGGGCCTGGGGCTTTTTCTGCTGCACCCCCGaa                                                                     >  1:533881/1‑150 (MQ=255)
                                                                                   tctcCGGATTTTTATTGCTGGCGGTTGGCATGGCGGGATTTAACCCCGGGCCTGG‑CGGCAGTAAAGGG‑ACC‑GGGTTTTTGGCGCTGGCCATTGCAAGCCCCGTCGCCGGTGGCGTGGGGGGTACGGGGCTTTTTTTGCTTCCCGCCCGaa                                                                     >  1:1235499/1‑150 (MQ=255)
                                                                                    ctcCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTA‑TTAACCCTGGCGCTGGACGGCAGTAAAGGT‑ACA‑GGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGaa                                                                     <  1:909881/149‑1 (MQ=255)
                                                                                    ctcCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTA‑TTAACCCTGGCGCTGGACGGCAGTAAAGGT‑ACA‑GGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGaa                                                                     <  1:1223306/149‑1 (MQ=255)
                                                                                      ccGGATTTTTATTGCTGGCGGTTGGCATGGCGGTA‑TTAACCCTTGCGCTGGGCCGCCGTTAAGGTGACG‑GGTTT‑TATCCGGTGACGGTTGGAGGGCCGGTGGCCGATTGCGTGGGGGGGCTGGGGGTTTTTCTGCTGTTAGCCAGAAATa                                                                  >  2:469375/1‑150 (MQ=255)
                                                                                      ccGGATTTTTATTGCTGGCGGTTGGCATGGCGGTA‑TTAACCCTGGCGCTGGACGGCAGTAAAGGT‑ACA‑GGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAAt                                                                   <  2:635405/149‑1 (MQ=255)
                                                                                      ccGGATTTTTATTGCTGGCGGTTGGCATGGCGGTA‑TTAACCCTGGCGCTGGACGGCAGTAAAGGT‑ACA‑GGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATa                                                                  <  1:441987/150‑1 (MQ=255)
                                                                                       cGGATTTTTATTGCTGGCGGTTGGCATGGCGGTA‑TTTAACCTTGCGGTGGACGGGAGGAAAGGG‑ACA‑AGGTTTTTGCCCCCGACGGTTGCAGGGCTGGGCGCAGTTGGCGTGGGGGGCCTTGGGCTTTTTTTGCTTGCCGCCCCAAATaa                                                                 >  1:644069/1‑150 (MQ=255)
                                                                                          aTTTTTATTGCTGGCGGTTGGCATGGCGGGA‑TTTACCCTTGGGCTGGTAGGCAGGAAAGGG‑ACA‑GGTTTTTTGCCGCCGAAGATTGCAGGCCTGGGCGCAGTTGGCGGGGTGGGCCTGGGGCTTTTTTTTCTGCACCGCCGAAATaaaaa                                                              >  1:42223/1‑150 (MQ=255)
                                                                                            ttttATTGCTGGCGGTTGGCATGGCGGGA‑TTTAACCCGGCGCTGGACGGCAGTAAAGGG‑ACA‑GGGTTTTTGCCGCTGGCGGTTGCAGGCCTGGGCGCAGGTGGCGTGGTGGGACCGGGGCTTTTTCTGCTGCCCCCCCGAAATaaaaaac                                                            >  2:1077054/1‑148 (MQ=255)
                                                                                            ttttATTGCTGGCGGTTGGCATGGCGGGA‑TTAAACCTGGGGCTTGGCGGCAGGAAAAGG‑ACA‑GGGTTATTGCCGCCGAAGATTGGAGGCCTGGTCGCAGGTGGGGGGGGGGCACTGGGGGTTTATTTGCTGCACGGCAGAAAAaaccaac                                                            >  1:137679/1‑146 (MQ=255)
                                                                                            ttttATTGCTGGCGGTTGGCATGGCGGGA‑TTAAACCTGGCGCTGGACGGCAGTAAAGGG‑ACA‑GGGTTTTCGCCGCCGACGATTGCCGGCCTGGTCGCAGTTGGCGGGGTGGCCCTGGTGCTTTTTCTGCTGCCCGCCCGAAATAACAAcc                                                            >  2:1151227/1‑150 (MQ=255)
                                                                                             tttATTGCTGGCGGTTGGCATGGCGGGA‑TTTAACCCGGCGCCGGAAGGCAGGAAAGGG‑ACA‑GGGTTATCGCCGCTGACGATTTCAGGCCTGGGCGCAGTTGGGGGGGTGGGACTGGTGCTTTATCTGCTGCACCCCCGAAATaaaaaac                                                            >  2:374980/1‑147 (MQ=255)
                                                                                                aTTGCTGGCGGTTGGCATGGCGGTA‑TTAACCCCGGGGCTGGGCGGCCGTAAAGGG‑ACA‑GGGTTATCGCCGCTGACGATTGCAGGCCTGGGCGCAGTTGGGGGGGGGGGACTGGTGGTTTTTCTTCTGCCCGCCAGAAAAAAAAAACGGGcc                                                       >  1:809214/1‑151 (MQ=255)
                                                                                                    cTGGCGGTTGGCATGGCGGGA‑TTAACCCTTGCGCTGGACCGCAGTAAAGGG‑ACA‑GGGTTATCGCCGCTGACGGTTGCCGGCCTGGGCGCAGTTTGCGTGGTGGCACTGGGGCTTTATTTGCTGCGCGCCAGAAATAACAACCGTGCCCTGt                                                   >  1:513909/1‑151 (MQ=255)
                                                                                                        cGGTTGGCATGGCGGTA‑TTAACCCTGGCGCTGGACGGCAGTAAAGGT‑ACA‑GGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTTCCCGTCATAACTAACAACCGTGTCCTGTGCAg                                               <  1:1214567/151‑1 (MQ=255)
                                                                                                        cGGTTGGCATGGCGGAA‑TTAACCCTGGGGCTGGACGGCCGTAACGGG‑ACA‑GGTTTATCGCCGCGGACGATTGCAGGCCTGGTCGCAGGTGGGGGGGGGGCACTGGTGCTTTATCTGCTGCACGCCAGGAATAACAACCGTGCCCTGtt                                                  >  1:445666/1‑148 (MQ=255)
                                                                                                             ggCATGGCGGTA‑TTAACCCTGGCGCTGGACGGCAGTAAAGGT‑ACA‑GGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTg                                           <  1:19921/150‑1 (MQ=255)
                                                                                                                   gCGGTA‑TTAACCCTGGCGCTGGACGGCAGTAAAGGT‑ACA‑GGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGt                                    <  2:1292293/151‑1 (MQ=255)
                                                                                                                     ggTA‑TTAACCCTGGCGCTGGACGGCAGTAAAGGT‑ACA‑GGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTc                                  >  1:567498/1‑151 (MQ=255)
                                                                                                                           tAACCCTGGCGCTGGACGGCAGTAAAGGT‑ACA‑GGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTa                              <  1:371874/150‑1 (MQ=255)
                                                                                                                           tAACCCTGGCGCTGGACGGCAGTAAAGGT‑ACA‑GGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTAc                             >  1:314685/1‑151 (MQ=255)
                                                                                                                              cccTGGCGCTGGACGGCAGTAAAGGT‑ACA‑GGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCg                          <  2:15121/151‑1 (MQ=255)
                                                                                                                                       tGGACGGCAGTAAAGGT‑ACA‑GGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTc                 >  2:800609/1‑151 (MQ=255)
                                                                                                                                        ggACGGCAGTAAAGGT‑ACA‑GGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCATGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTc                 <  2:63301/150‑1 (MQ=255)
                                                                                                                                            ggCAGTAAAGGT‑ACA‑GGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTCGCTgg            <  2:270113/151‑1 (MQ=255)
                                                                                                                                            ggCAGTAAAGGT‑ACA‑GGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTCGCTgg            <  2:411262/151‑1 (MQ=255)
                                                                                                                                               aGTAAAGGT‑ACA‑GGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTCGCTGGGcc         >  1:422753/1‑151 (MQ=255)
                                                                                                                                                     ggT‑ACA‑GGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTCGCTGGGCCTGGCgg   >  2:1005554/1‑151 (MQ=255)
                                                                                                                                                      gT‑ACA‑GGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTCGCTGGGCCTGGCggg  <  1:155624/151‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                       
CAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTA‑TTAACCCTGGCGCTGGACGGCAGTAAAGGT‑ACA‑GGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTCGCTGGGCCTGGCGGG  >  NC_000913/2161970‑2162269

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