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→GGGiceTputative 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

CAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAACCCTGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTCGCTGGGCCTGGCGGG  >  NC_000913/2161970‑2162269
                                                                                                                                                       |                                                                                                                                                       
cAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAACCCTGGCGCTGGACGGCAGTAAAGGt                                                                                                                                                         <  2:707537/151‑1 (MQ=255)
 aaCATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAACCCTGGCGCTGGACGGCAGTAAAGGTa                                                                                                                                                       <  1:1176760/151‑1 (MQ=255)
            ggggATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAACCCTGGCGCTGGACGGCAGTAAAGGTACAGGTTTATcg                                                                                                                                           <  2:887516/151‑1 (MQ=255)
             gggATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAAGTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAACCCTGGCGCTGGACGGCAGTAAAGGTACAGGTTTATc                                                                                                                                            <  1:339972/149‑1 (MQ=255)
                     cGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTTAACCTGGGGCCGGAAGGGAGGAAAGGGAAAGGGTTATtgcc                                                                                                                                         >  2:37074/1‑144 (MQ=255)
                         gCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAACCCTGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTg                                                                                                                              <  2:147391/151‑1 (MQ=255)
                            aTCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGGATTAAACCTGGGGCTGGGCGGGAGTAAAGGGACAAGTTTTTCGCCGCTGGCGGtt                                                                                                                               >  1:102442/1‑147 (MQ=255)
                               gCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGGATTTAACCCCGCGCTGGGGGGCAGTTAAGGGTCAGGGTTATCGCCGCCGGCCATTGCCGGcc                                                                                                                        >  1:530591/1‑151 (MQ=255)
                                  acaTTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGGATTTAACCTGGGGCTGGGCGGCAGTAAAAGGACAGGGTTTTTGCCGCTGGCGAtttccgggccgg                                                                                                                     >  2:337281/1‑141 (MQ=255)
                                  acaTTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGGATTAACCCTGGCGGTGGGAGGCAGGAAAGGGACAGGGTTTTCGCCGCCGACGATTGCCGGGCTgg                                                                                                                     >  1:741200/1‑151 (MQ=255)
                                   caTTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCCTGGCGGGATTAAACCTGGGGCTGGGAGGCAGGAAAAGGGTAGGGTTTTCGCCGCTGACGATTTCGGGCCCggg                                                                                                                    >  1:290224/1‑150 (MQ=255)
                                   caTTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAACCCTGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTgg                                                                                                                     <  2:290241/150‑1 (MQ=255)
                                           ttAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGGATTTAACCCGGCGCTGGGCGGCAGGAAAGGGACAGGGTTTTCGGCGCCGAACGTTTCAGGCCCGGGCGGAggt                                                                                                             >  1:864478/1‑148 (MQ=255)
                                                gCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAACCCTGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCgtg                                                                                                       <  1:7960/151‑1 (MQ=255)
                                                gCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAACCCTGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCgt                                                                                                        <  1:730880/150‑1 (MQ=255)
                                                  cGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAACCCTGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCgtggt                                                                                                     <  2:311736/151‑1 (MQ=255)
                                                  cGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAACCCTGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCgtggt                                                                                                     <  2:490093/151‑1 (MQ=255)
                                                        acacCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAACCCTGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCa                                                                                                 <  2:807110/149‑1 (MQ=255)
                                                          acCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTTACCCTGGCGCTCGGCGGCAGTAAAGGGACAGGGTTTTTCGCGCTGACGGTTGCAGGGCCGGTCGCAGGTGGCGTGGGGGCCCTg                                                                                              >  1:893814/1‑150 (MQ=255)
                                                             aTGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAACCCTGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGc                                                                                          <  2:1153555/151‑1 (MQ=255)
                                                                 aGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGGATTAAACCTGGGGCTGGACGGGAGGAAAGGGACAGGGTTATTGCCGCTGGCGATTGCAGGGCTGGGCGCAGTTGGCGGGGGGGCAcctggggttt                                                                                       >  2:1014025/1‑141 (MQ=255)
                                                                    cgcgGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAACCCTGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATct                                                                                   <  1:1279361/151‑1 (MQ=255)
                                                                    cgcgGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAACCCTGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATc                                                                                    <  1:1130057/150‑1 (MQ=255)
                                                                    cgcgGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAACCCCGGCGCTTGGCGGCAGTAAAGGGACAGGGTTTTTGCCGCTGACGATTGCCGGCCTGGGCGCAGGTGGCGGGGGGGCCCTGGGGCttttt                                                                                     >  1:800388/1‑147 (MQ=255)
                                                                    cgcgGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGGATTAACCCTGGGGCTGGGCGGCAGGAAAGGGACAGGGTTTTCGCCGCTGACGATTGCAGGCCTGGTCGCCGGTGGGGGGGGGGCACCGGTGCttt                                                                                       >  2:902779/1‑147 (MQ=255)
                                                                    cgcgGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGGATTAAACCTTGCGCCGGGCGGCAGTAAAGGGACAGGGTTATCGCCGCTGGCGATTGCAGGCCTGGTCGCAGGTGGCGTGGGGGGACTGGGGCTTTATct                                                                                   >  2:380600/1‑151 (MQ=255)
                                                                    cgcgGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGCTGGCATGGCGGTATTAACCCTGGCGCTGGACGGCAGTAAAGGGACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGGTGGCGGGGTGGCACTGGTGCTTTAt                                                                                     >  1:358700/1‑149 (MQ=255)
                                                                     gcggcgCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAACCCTGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATctg                                                                                  <  2:1072283/151‑1 (MQ=255)
                                                                       ggcgCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAACCCTGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATctgct                                                                                <  1:127353/151‑1 (MQ=255)
                                                                        gcgcTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGGATTAAACCTCGGGCTGGGCCGCAGAAAAAGAACCGGGTTATTGCCCCTGACGGTTGCAGGCCTGGTGGGAGGGGGCGGGGGGGCACTgggggcttatttact                                                                                >  1:492196/1‑137 (MQ=255)
                                                                           cTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTTAACCTGGGGCCGGGAGGCAGTAAAGGGACAGGGTTTTCGCCGCTGACGGATGGAGGGCTGGTGGCAGTTGGCGTGGGGGCACTGGGGCTTTATctgctg                                                                               >  1:223865/1‑148 (MQ=255)
                                                                               gATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAACCCTGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGcc                                                                         <  1:983210/150‑1 (MQ=255)
                                                                                  ctctcCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAACCCTGGCGCTGGACGGGAGGAAAGGGACAGGGTTTTTGCCCCTGACGATTTCAGGCCTGGGCGCAGGTGGCGGGGGGGCACTGGGGCTTTTTTTGCTGCACGcc                                                                         >  1:411454/1‑147 (MQ=255)
                                                                                   tctcCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTTAACCCTGCGCCGGGCGGGAGTAAAAGGACAGGGTTTTTGCCGCTGACGGTTTCCGGGCCGGGCGCCGGTTGGGTGGGGGGCCTGGGGCTTTTTCTGCTGCACCCCCGaa                                                                     >  1:533881/1‑150 (MQ=255)
                                                                                   tctcCGGATTTTTATTGCTGGCGGTTGGCATGGCGGGATTTAACCCCGGGCCTGGCGGCAGTAAAGGGACCGGGTTTTTGGCGCTGGCCATTGCAAGCCCCGTCGCCGGTGGCGTGGGGGGTACGGGGCTTTTTTTGCTTCCCGCCCGaa                                                                     >  1:1235499/1‑150 (MQ=255)
                                                                                    ctcCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAACCCTGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGaa                                                                     <  1:909881/149‑1 (MQ=255)
                                                                                    ctcCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAACCCTGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGaa                                                                     <  1:1223306/149‑1 (MQ=255)
                                                                                      ccGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAACCCTTGCGCTGGGCCGCCGTTAAGGTGACGGGTTTTATCCGGTGACGGTTGGAGGGCCGGTGGCCGATTGCGTGGGGGGGCTGGGGGTTTTTCTGCTGTTAGCCAGAAATa                                                                  >  2:469375/1‑150 (MQ=255)
                                                                                      ccGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAACCCTGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAAt                                                                   <  2:635405/149‑1 (MQ=255)
                                                                                      ccGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAACCCTGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATa                                                                  <  1:441987/150‑1 (MQ=255)
                                                                                       cGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTTAACCTTGCGGTGGACGGGAGGAAAGGGACAAGGTTTTTGCCCCCGACGGTTGCAGGGCTGGGCGCAGTTGGCGTGGGGGGCCTTGGGCTTTTTTTGCTTGCCGCCCCAAATaa                                                                 >  1:644069/1‑150 (MQ=255)
                                                                                          aTTTTTATTGCTGGCGGTTGGCATGGCGGGATTTACCCTTGGGCTGGTAGGCAGGAAAGGGACAGGTTTTTTGCCGCCGAAGATTGCAGGCCTGGGCGCAGTTGGCGGGGTGGGCCTGGGGCTTTTTTTTCTGCACCGCCGAAATaaaaa                                                              >  1:42223/1‑150 (MQ=255)
                                                                                            ttttATTGCTGGCGGTTGGCATGGCGGGATTTAACCCGGCGCTGGACGGCAGTAAAGGGACAGGGTTTTTGCCGCTGGCGGTTGCAGGCCTGGGCGCAGGTGGCGTGGTGGGACCGGGGCTTTTTCTGCTGCCCCCCCGAAATaaaaaac                                                            >  2:1077054/1‑148 (MQ=255)
                                                                                            ttttATTGCTGGCGGTTGGCATGGCGGGATTAAACCTGGGGCTTGGCGGCAGGAAAAGGACAGGGTTATTGCCGCCGAAGATTGGAGGCCTGGTCGCAGGTGGGGGGGGGGCACTGGGGGTTTATTTGCTGCACGGCAGAAAAaaccaac                                                            >  1:137679/1‑146 (MQ=255)
                                                                                            ttttATTGCTGGCGGTTGGCATGGCGGGATTAAACCTGGCGCTGGACGGCAGTAAAGGGACAGGGTTTTCGCCGCCGACGATTGCCGGCCTGGTCGCAGTTGGCGGGGTGGCCCTGGTGCTTTTTCTGCTGCCCGCCCGAAATAACAAcc                                                            >  2:1151227/1‑150 (MQ=255)
                                                                                             tttATTGCTGGCGGTTGGCATGGCGGGATTTAACCCGGCGCCGGAAGGCAGGAAAGGGACAGGGTTATCGCCGCTGACGATTTCAGGCCTGGGCGCAGTTGGGGGGGTGGGACTGGTGCTTTATCTGCTGCACCCCCGAAATaaaaaac                                                            >  2:374980/1‑147 (MQ=255)
                                                                                                aTTGCTGGCGGTTGGCATGGCGGTATTAACCCCGGGGCTGGGCGGCCGTAAAGGGACAGGGTTATCGCCGCTGACGATTGCAGGCCTGGGCGCAGTTGGGGGGGGGGGACTGGTGGTTTTTCTTCTGCCCGCCAGAAAAAAAAAACGGGcc                                                       >  1:809214/1‑151 (MQ=255)
                                                                                                    cTGGCGGTTGGCATGGCGGGATTAACCCTTGCGCTGGACCGCAGTAAAGGGACAGGGTTATCGCCGCTGACGGTTGCCGGCCTGGGCGCAGTTTGCGTGGTGGCACTGGGGCTTTATTTGCTGCGCGCCAGAAATAACAACCGTGCCCTGt                                                   >  1:513909/1‑151 (MQ=255)
                                                                                                        cGGTTGGCATGGCGGTATTAACCCTGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTTCCCGTCATAACTAACAACCGTGTCCTGTGCAg                                               <  1:1214567/151‑1 (MQ=255)
                                                                                                        cGGTTGGCATGGCGGAATTAACCCTGGGGCTGGACGGCCGTAACGGGACAGGTTTATCGCCGCGGACGATTGCAGGCCTGGTCGCAGGTGGGGGGGGGGCACTGGTGCTTTATCTGCTGCACGCCAGGAATAACAACCGTGCCCTGtt                                                  >  1:445666/1‑148 (MQ=255)
                                                                                                             ggCATGGCGGTATTAACCCTGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTg                                           <  1:19921/150‑1 (MQ=255)
                                                                                                                   gCGGTATTAACCCTGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGt                                    <  2:1292293/151‑1 (MQ=255)
                                                                                                                     ggTATTAACCCTGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTc                                  >  1:567498/1‑151 (MQ=255)
                                                                                                                           tAACCCTGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTa                              <  1:371874/150‑1 (MQ=255)
                                                                                                                           tAACCCTGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTAc                             >  1:314685/1‑151 (MQ=255)
                                                                                                                              cccTGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCg                          <  2:15121/151‑1 (MQ=255)
                                                                                                                                       tGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTc                 >  2:800609/1‑151 (MQ=255)
                                                                                                                                        ggACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCATGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTc                 <  2:63301/150‑1 (MQ=255)
                                                                                                                                            ggCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTCGCTgg            <  2:270113/151‑1 (MQ=255)
                                                                                                                                            ggCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTCGCTgg            <  2:411262/151‑1 (MQ=255)
                                                                                                                                               aGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTCGCTGGGcc         >  1:422753/1‑151 (MQ=255)
                                                                                                                                                     ggTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTCGCTGGGCCTGGCgg   >  2:1005554/1‑151 (MQ=255)
                                                                                                                                                      gTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTCGCTGGGCCTGGCggg  <  1:155624/151‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                       
CAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGGCGCTTTGATCTCTCCGGATTTTTATTGCTGGCGGTTGGCATGGCGGTATTAACCCTGGCGCTGGACGGCAGTAAAGGTACAGGTTTATCGCCGCTGACGATTGCAGGCCTGGTCGCAGTTGGCGTGGTGGCACTGGTGCTTTATCTGCTGCACGCCAGAAATAACAACCGTGCCCTGTTCAGTCTGAAACTGTTCCGTACTCGTACCTTTTCGCTGGGCCTGGCGGG  >  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: