Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913948,5280AG34.2% 52.2 / 32.5 44I434V (ATC→GTC) ycaMputative transporter YcaM
Reads supporting (aligned to +/- strand):  ref base A (13/14);  new base G (7/7);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGTGGCGATTCTGATTATGCTGCCGACTCTCGGCATTGGTGATATGAACAATCTCTATAAATGGCTGTTGAACCTTAATTCGGTAGTGATGCCGCTGCGTTATTTATGGGTATTTGTTGCATTTATTGCAGTCGTTCGCTTGGCGCAGAAATATAAACCAGAATATGTCTTTATTCGTAATAAGCCTCTGGCAATGACCGTCGGGATTTGGTGTTTTGCCTTTACCGCCTTTGCCTGTTTGACGGGGATCTTCCCGAAAATGGAAGCCTTCACTGCAGAGTGGACCTTCCAGTTGGCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGACTCTGGCAGATCTGACCCGTTTTTACATAATCGAGAAATGCGCGCAAGCCGCAGGACATATACTGGCGGTTTGGATAATAGATCTGGAAGCCT  >  NC_000913/948281‑948774
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
ggtggCGATTCTGATTATGCTGCCGACTCTCGGCATTGGTGATATGAACAATCTCTATAAATGGCTGTTGAACCTTAATTCGGTAGTGATGCCGCTGCGTTATTTATGGGTATTTGTTGCATTTATTGCAGTCGTTCGCTTGGCGCAGAAATATAAACCAGAATATGTCTTTATTCGTAATAAGCCTCTGGCAATGACCGTCGGGATTTGGTGTTTTGCCTTTACCGCCTTTGCCTGTTTGACGGGGAt                                                                                                                                                                                                                                                        <  1:145987/249‑1 (MQ=255)
   ggCGATTCTGATTATGCTGCCGACTCTCGGCATTGGTGATATGAACAATCTCTATAAATGGCTGTTGAACCTTAATTCGGTAGTGAGGCCGCTGCGTTATTTATGGGTATTTGTTGCATTTATTGCAGTCGTTCGCTTGGCGCATAAATATAAACCAGAATATGTCTTTATTCGTAATAAGCCTCTGGCAATGACCGTCGGGATTTGGTGTTTTGCCTTTACCGCCGTTGCCGGTTTGACGGGGATCTTcc                                                                                                                                                                                                                                                   <  2:391915/251‑1 (MQ=255)
         tCTGATTATGCTGCCGACTCTCGGCATTGGTGATATGAACAATCTCTATAAATGGCTGTTGAACATTAATTCGGTAGTGATGCCGATGCGTTATTTATGGGTATTTGTTGCATTTATTGCAGTCGTTCGCTTGGCGCAGACATATAAACAAGAATATGTCTTTATTCGTAATAAGCCTCTGGCAATGACCGTCGGGATTTGGGGTTTTGCCTTTACCGCCTTTGCCTGTTTGACGGGGGTCTTCCCGaaa                                                                                                                                                                                                                                              >  1:604841/1‑250 (MQ=255)
                aggctgcCGACTCTCGGCATTGGTGAGATGAACAATCTCTATAAATGGCTGTTGAACCTTAATTCGGTAGTGATGCCTCTGCGCTATTTATGGGCATTTGTTGCATTTATTGCAGTCGTTCGCTTGGCGCAGAAATATAAACCAGAATCTGTCTTTATTCGTAATAAGCCTCTGGCAATGAACGACGGGATTTGGTGTTTTGCCTTTACCGCCTTTGCCTGTTTGACGGGGCTCTTCCAGGAAATGGAAGc                                                                                                                                                                                                                                      >  2:58404/3‑251 (MQ=255)
                   ctgcCGACTCTCGGCATTGGTGATATGAACAATCTCTATAAATGGCTGTTGAACCTTAATTCGGTAGTGATGCCGCTGCGTTATTTATGGGTATTTGTTGCATTTATTGCAGTCGTTCGCTTGGCGCAGAAATATAAACCAGAATATGTCTTTATTCGTAATAAGCCTCTGGCAATGACCGTCGGGATTTGGTGTTTTGCCTTTACCGCCTTTGCCTGTTTGACGGGGGTCTTCCCGAAAATGGAAGCCtt                                                                                                                                                                                                                                   >  2:230787/1‑251 (MQ=255)
                        gACTCTCGGCAGGGGTGATATGAACAATCTCTTTAAAGGGCTGTTGAACCTTAATTGGGTAGTGATGCCGCTGCGTGATGTATGGGTATTTGTGGAATTTATTGCAGTCGTTCGCTTGGCGCAGAAATATAAACCAGAATATGTCTTTATTCATAATAAGCCTCTGGCAATGACCGTCGGGATTTGGTGTTTTGCCTTTACAGCCTGTGCCTGTTTGTCGGGGATCTTCCCGAAAATGGAAGCCTTCACtt                                                                                                                                                                                                                              <  2:321622/251‑2 (MQ=255)
                         aCTCTCGGCATTGGTAATATGAACAATCTCTATAAATGGCTGTTGAACCTTAATTCGGTAGTGATGCCGCTGCGTTATTTATGGGTATTTGTTGCATTTATTGCAGTCGTTCGCTTGGCGCAGAAATATAAACCAGAATATGTCTTTATTCGTAATAAGCCTCTGGCAATGACCGTCGGGATTTGGTGTTTTGCCTTTACCGCCTTTGCCTGTTTGACGGGGATCTTCCCGAAAATGGAAGCCTTCACTGc                                                                                                                                                                                                                             <  1:156424/251‑1 (MQ=255)
                            ctcGGCATTGGTGATATGAACAATCTCTATAAATGGCTGTTGAACCTTAATTCGGTAGTGATGCCGCTGCGTTATTTATGGGTATTTGTTGCATTTATTGCAGTCGTTCGCTTGGCGCAGAAATATAAACCAGAATATGTCTTTATTCGTAATAAGCCTCTGGCAATGACCGTCGGGATTTGGTGTTTTGCCTTTACCGCCTTTGCCTGTTTGACGGGGATCTTCCCGAAAATGGAAGCCTTCACTGc                                                                                                                                                                                                                             >  2:206281/1‑248 (MQ=255)
                                       tGATATGAACAATCTCTATAAATGGCTGTTGAACCTTAATTCGGTAGTGATGCCGCTGCGTTATTTATGGGTATTTGTTGCATTTATTGCAGTCGTTCGCTTGGCGCAGAAATATAAACCAGAATATGTCTTTATTCGTAATAAGCCTCTGGCAATGACCGTCGGGATTTGGTGTTTTGCCTTTACCGCCTTTGCCTGTTTGACGGGGATCTTCCCGAAAATGGAAGCCTTCACTGCAGAGTGGACCTTcc                                                                                                                                                                                                               >  1:6228/1‑251 (MQ=255)
                                        gATATGAACAATCTCTATAAATGGCTGTTGAACCTTAATTCGGTAGTGATGCTGCTGCGTTATTTATGGGTATTTGTTGCATTTATTGCAGTCGTTCGCTTGGCGCAGAAATATAAACCAGAATATGTCTTTATTCGTAATAAGCCTCTGGCAATGACCGTCGGGATTTGGTGTTTTGCCTTTACCGCCTTTGCCTTTTTGACGGGGATCTTCCCGAACATGGAAGCCTTCACTGCAGAGTGGACCTTccc                                                                                                                                                                                                              >  1:196557/1‑250 (MQ=255)
                                           atGAACAATCTCTATAAATGGCTGTTGAACCTTAATTCGGTAGTGATGCCGCTGCGTTATTTATGGGTATTTGTTGCATTTATTGCAGTCGTTCGCTTGGCGCAGAAATATAAACCAGAATATGTCTTTATTCGTAATAAGCCTCTGGCAATGACCGTCGGGATTTGGTGTTTTGCCTTTACCGCCTTTGCCTGTTTGACGGGGGTCTTCCCGAAAATGGAAGCCTTCACTGCAGAGTGGACCTTCCAGtt                                                                                                                                                                                                           >  2:386317/1‑251 (MQ=255)
                                             gAACACTCTCTATAACTGGCTGTTGAACCTTAATTCGGTAGTGATGCCGCTGCGTTATTTATGGGTATTTGTTGCATTTATTGCAGTCGTTCGCTTGGCGCAGAAATATAAACCAGAATATGTCTTTATTCGTAATAAGCCTCTGGCAATGACAGTCGGGATTTGGTGTTTTGCCTTTACCGCCTTTGCCTGTTTGACGGGGATCTTCCCGAAAATGGAAGCCTTCACTGCCGAGTGGACCTTCCAGTTgg                                                                                                                                                                                                         >  2:237616/1‑251 (MQ=255)
                                                             tGGCTGTTGAACCTTAATTCGGTAGTGATGCCGCTGCGTTATTTATGGGTATTTGTTGCATTTATTGCAGTCGTTCGCTTGGAGAAGAAATATAAACCAGAATATGACTTTATTCGTAATAAGCCTCTGGCAATGACCGTCGGGATTTGGTGTTTTGCCTTTACCGCCTTTGCCTGTTTGACGGGGATCTTCCCGAAAATGGACGCCTTCACTGCAGAGTGGACCTTCCAGTTGGCGCTGAATGTTGCacc                                                                                                                                                                                         >  1:512570/1‑249 (MQ=255)
                                                                  gTTGAACCTTAATTCGGTAGTGATGCCGCTGCGTTATTTATGGGTATTTGTTGCATTTATTGCAGTCGTTCGCTTGGCGCAGAAATATAAACCAGAATATGTCTTTATTCGTAATAAGCCTCTGGCAATGACCGTCGGGATTTGGTGTTTTGCCTTTACCGCCTTTGCCTGTTTGACGGGGATCTTCCCGAAAATGGAAGCCTTCACTGCAGAGTGGACCTTCCAGTTGGCGCTGAATGTTGCAACGCCg                                                                                                                                                                                     <  2:464076/250‑1 (MQ=255)
                                                                       aCCTTAATTCGGTAGTGATGCCGCTGCGTTATTTATGGGTATTTGTTGCATTTATTGCAGTCGTTCGCTTGGCGCAGAAATATAAACCAGAATATGTCTTTATTCGTAATAAGCCTCTGGCAATGACCGTCGGGATTTGGTGTTTTGCCTTTACCGCCTTTGCCTGTTTGACGGGGGTCTTCCCGAAAATGGAAGCCTTCACTGCAGAGTGGACCTTCCAGTTGGCGCTGAATGTTGCAACGCCGTTtgtg                                                                                                                                                                               >  1:127255/1‑251 (MQ=255)
                                                                                            cGCTGCGTTATTTATGGGTATGTGTTGCATTTATTGCAGTCGTTCGCTTGGCGCAGAAATATAAACCAGAATATGTCTTTATTCGTAATAAGCCTCTGGCAATGACCGTCGGGATTTGGTGTTTTGCCTTTACCGCCTTTGCCTGTTTGACGGGGATCTTCCCGAAAATGGAAGCCTGCACTGCAGAGTGGACCTTCCAGTTGGCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGata                                                                                                                                                          <  2:450824/251‑1 (MQ=255)
                                                                                                gCGTTATTTATGGGTATTTGTTGCATTTATTGCAGTCGTTCGCTTGGCGCAGAAATATAAACCAGAATATGTCTTTATTCGTAATAAGCCTCTGGCAATGACCGTCGGGATTTGGTGTTTTGCCTTTACCGCCTTTGCCTGTTTGACGGGGGTCTTCCCGAAAATGGAAGCCTTCACTGCAGAGTGGACCTTCCAGTTGGCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTc                                                                                                                                                       <  2:143163/250‑1 (MQ=255)
                                                                                                             gTATTTGTTGCATTTATTGCAGTCGTTCGCTTGGCGCAGAAATATAAACCAGAATATGTCTTTATTCGTAATAAGCCTCTGGCAATGACCGTCGGGATTTGGTGTTTTGCCTTTACCGCCTTTGCCTGTTTGACGGGGATCTTCCCGAAAATGGAAGCCTTCACTGCAGAGTGGACCTTCCAGTTGGCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGcgcg                                                                                                                                         >  1:588873/1‑251 (MQ=255)
                                                                                                                     ggCATTTATTGCTGTCGTTCGCTTGGCGCAGAAATATAAACCAGAAGATGTCTTTATTCGTAATAAGCCTCTGGCAATGACCGTCTGGATTTGGTGTTGTGCCTTTACCGCCTTTGCCTGTTTGACGGGGATCTTCCCGAAAATGGAAGCCTTCACTGCAGAGTGGACCTTCCAGTTGGCGCTGAATGTTGCAACGCCGTTTTTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCg                                                                                                                                  <  2:588873/249‑1 (MQ=255)
                                                                                                                        aTTTATTGCAGTCGTTCGCTTGGCGCAGAAATATAAACCAGAATATGTCTTTATTCGTAATAAGCCTCTGGCAATGACCGTCGGGATTTGGTGTTTTGCCTTTACCGCCTTTGCCTGTTTGACGGGGGTCTTCCCGAAAATGGAAGCCTTCACTGCAGAGTGGACCTTCCAGTTGGCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGAATa                                                                                                                              <  1:320840/251‑1 (MQ=255)
                                                                                                                           tATTGCAGTCGTTCGCTTGGCGCAGAAATATAAACCAGAATATGTCTTTATTCGTAATAAGCCTCTGGCAATGACCGTCGGGATTTGGTGTTTTGCCTTTACCGCCTTTGCCTGTTTGACGGGGGTCTTCCCGAAAATGGAAGCCTTCACTGCAGAGTGGACCTTCCAGTTGGCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGAATAGt                                                                                                                            <  2:401222/250‑1 (MQ=255)
                                                                                                                                            tGGCGCATAAATATAAACCAGAATATGTCTTTATTCGTAATAAGCCTCTGGCAATGACCGTCGGGATTTGGTGTTTTGCCTTTGCCGCCTTTGCCTGTTTGACGGGGATCTTCCCGAAAATGGAAGCCTTCACTGCAGAGTGGACCTTCGAGTTGGCGCGGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCTGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTc                                                                                                           <  1:568960/250‑1 (MQ=255)
                                                                                                                                            tGGCGCAGAAATAGAAACCAGAATATGTCTTTATTCGTAATAAGCCTCTGGCAATGACCGTCGGGATTTGGTGTTTTGCCTTTACCGCCTTTGCCTGGTTGACGGGGATCTTCCCGAAAATGGAAGCCTTCACTGCAGAGTGGACCTTCCAGTTGGCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTc                                                                                                           <  2:298733/250‑1 (MQ=255)
                                                                                                                                            gggCGCAGAAAGATAAACCAGAATATGTCTTTATTCGTAATAAGCCTCTGGCAATGACCGTCGGGATTTGGTGTTTTGCCGTTACCGCCTTTGCCTGTTTGACGGGGATCTTCCCGAAAATGGAAGCCTTCACTGCAGAGTGGACCTTCCAGTTGGCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGt                                                                                                            <  2:198975/248‑1 (MQ=255)
                                                                                                                                                       tataAACCAGAATATGTCTTTATTCGTAATAAGCCTCTGGCAAAGACAGTCGCGATTTGGTGTTTTGCCTTTACCGCCTTTCCCTGTTCGACGGGACTCTTCCCGAAAATGGACGCCTTCACTGCAGAGTGGACCTTCCAGCTTGCGCTGAATGTTGCAACGCCGTTTGTGCTGGCAGGATTAGGACTGCTATTCCCGCTGCTGGCGAGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGAc                                                                                               >  1:312213/1‑251 (MQ=255)
                                                                                                                                                          aaaCCAGAATATGTCTTTATTCGTAATAAGCCTCTGGCAATGACCGTCGGTATTTGGTGTTTTGCCTTTACCGCCTTTGCCTGTTTGACGGGGGTCTTCCCGAAAATGGAAGCCTTCACTGCAGAGTGGACCTTCCAGTTGGCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGActct                                                                                            <  1:386317/251‑1 (MQ=255)
                                                                                                                                                                  atatGTCTTTATTCGTAATAAGCCTCTGGCAATGACCGTCGGGATTTGGTGTTTTGCCTTTACCGCCTTTGCCTGTTTGACGGGGGTCTTCCCGAAAATGGAAGCCTTCACTGCAGATTGGACCTTCCAGTTGGCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTCTTGTGGTCGGTGACGCTGCCTCGGGc                                                                                         >  2:121494/1‑246 (MQ=255)
                                                                                                                                                                       tCTTTATTCGTAATAAGCCTCTGGCAATGACCGTCGGGCTTTGGTGTTTTGCCTTTACCGCCTTTGCCTGTTTGACGGGGGTCTTCCCGAAAATGGAAGCCTTCACTGCAGAGTGGACCTTCCAGTTGGCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGACTAGTAAATAAGTATTGTGCTCGGTGACGCGGACTCTGGCAGATCTGAcc                                                                               >  2:420408/1‑251 (MQ=255)
                                                                                                                                                                                  aataaGCCTCTGGCAATGACCGTCGGGATTTGGTGTTTTGCCTTTACCGCCTTTGCCTGTTTGACGGGGATCTTCCCGAAAATGGAAGCCTTCACTGCAGAGTGGACCTTCCAGTTGGCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGACTCTGGCAGATCTGACCCGTTTTTACAt                                                                    >  2:139635/1‑251 (MQ=255)
                                                                                                                                                                                  aataaGCCTATAGAAATGACCGTCGGGATATGGTGTTTTGACTTTACACACTCAGCTTGTTTGACGGGGATATTCCCGCAAATGGAAGCCTTCACTGCAGAGTGGACCTTCCAGTTGGCGATGAATGTTGCAACGCCGTTTGTGCTGGTCGGATTCGGACTGATATTCACGTTGCTGGCGCGTAAAGCGAATGGTAAATAATTATTGTGGTCGGGGACGCCGACCATGGCAGCTCTGACCCGTTTTTACAt                                                                    >  2:296280/1‑251 (MQ=255)
                                                                                                                                                                                      aGCCTCTGGCAATGACCGTCGGGATTTGGTGTTTTGCCTTTACCGCCTTTGCCTGTTTGACGGGGGTCTTCCCGAAAATGGAAGCCTTCACTGCAGAGTGGACCTTCCAGTTGGCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGACTCTGGCAGATCTGACCCGTTTTTACATAATc                                                                <  1:199875/251‑1 (MQ=255)
                                                                                                                                                                                             ggCAATGACCGTCGGGATTTGGTGTTTTGCCTTTACCGCCTTTGCCTGTTTGACGGGGATCTTCCCGAAAATGGAAGCCTTCACTGCAGAGTGGACCTTCCAGTTGGCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGACTCTGGCAGATCTGACCCGTTTTTACATAATCGAGaa                                                           <  1:139635/249‑1 (MQ=255)
                                                                                                                                                                                               cAATGACCGTCGGGATTTGGTGTTTTGCCTTTACCGCCTTTGCCTGTTTGACGGGGGTCTTCCCGAAAATGGAAGCCTTCACTGCAGAGTGGACCTTCCAGTTGGCGCTGAATGTTGCAACGCCGGTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGAAGAGTAAATAATTATTGTGGTCGGTGACGCTGACTCTGGCAGATCTGACCCGTTTTTACATAATCGAGAAATgc                                                       <  2:332423/251‑1 (MQ=255)
                                                                                                                                                                                               cAATGACCGTCGGGATTTGGTGTTTTGCCTTTACCGCCTTTGCCTGTTTGACGGGGATCTTCCCGAAAATGGAAGCCTTCACTGCAGAGTGGACCTTCCAGTTGGCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGACTCTGGCAGATCTGACCCGTTTTTACATAATCGAGAAATgc                                                       >  1:591064/1‑251 (MQ=255)
                                                                                                                                                                                                aaTGACCGTCGGGATTTGGTGTTTTGCCTTTACCGCCTTTGCCTGTTTGACGGGGATCTTCCCGAAAATGGAAGCCTTCACTGCAGAGTGGACCTTCCAGTTGGCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAACGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGACTCTGGCAGATCTGACCCGTTTTTACATAATCGAGAAATgcg                                                      >  2:208433/1‑251 (MQ=255)
                                                                                                                                                                                                   gACCGTCGGGATTTGGTGTTTTGCCTTTACCGCCTTTGCCTGTTTGACGGGGATCTTCCCGAAAATGGAAGCCTTCACTGCAGAGTGGACCTTCCAGTTGGCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGAATAGGAAATAATTATTGTGGTCGGTGACGCTGACTCTGGCAGATCTGACCCGTTTTTACATAATCGAGAAATgcgcgc                                                   >  1:25355/1‑251 (MQ=255)
                                                                                                                                                                                                     ccGTCGGGATTTGGTGTTTTGCCTTTACCGCCTTTGCCTGTTTGACGGGGATCTTCCCGAAAATGGAAGCCTTCACTGCAGAGTGGACCTTCCAGTTGGCGCTGAATGTTGCAACGCCGGTTGTGCTGGTAGGATTAGGACTGATATGCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGACGCTGGCAGATCTGACCCGTTTTTACATAATCGAGAAATGCGCGCaa                                                 <  1:409937/251‑1 (MQ=255)
                                                                                                                                                                                                            gATTTGGTGTTTTGCCTTTACCGCCTTTGCCTGTTTGACGGGGATCTTCCCGAAAATGGAAGCCTTCACTGCAGAGTGGACCTTCCAGTTGGCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGACTCTGGCAGATCTGACCCGTTTTTACATAATCGAGAAATGCGCGCAAGCCGCAg                                          <  2:10437/251‑1 (MQ=255)
                                                                                                                                                                                                             aTTTGGTGTTTTGCCTTTACCGCCTTTGCCTGTTTGACGGGGGTCTTCCCGAAAATGGAAGCCTTCACTGCAGAGTGGACCTTCCAGTTGGCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGACTCTGGCAGATCTGACCCGTTTTTACATAATCGAGAAATGCGCGCAAGCCGCAgg                                         >  2:227412/1‑251 (MQ=255)
                                                                                                                                                                                                                          ccTTTACCGCCTTTGCCTGTTTGACGGGGATCTTCCCGAAAATGGAAGCCTTCACTGCAGAGTGGACCTTCCAGTTGGCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGACTCTGGCAGATCTGACCCGTTTTTACATAATCGAGAAATGCGCGCAAGCCGCAGGACATATACTGGCg                            >  2:318288/1‑251 (MQ=255)
                                                                                                                                                                                                                                ccGCCTTTGCCTGTTTGACGGGGATCTTCCCGAAAATGGAAGCCTTCACTGCAGAGTGGACCTTCCAGTTGGCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGACTCTGGCAGATCTGACCCGTTTTTACATAATCGAGAAATGCGCGCAAGCCGCAGGACATATACTGGCGGTTTgg                     <  1:206281/251‑1 (MQ=255)
                                                                                                                                                                                                                                 cGCCTTTGCCTGTTTGACGGGGCTCTTCCCGAACAGGGAAGCCTTCACTGCAGAGTGGACCTTCCAGTTGGCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGCCTGATATTCACGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGACTCTGGCAGATCTGACCCGTTTTTACATAATCGAGAAATGCGCGCAAGCCGCAGGAAATATACAGGCGGGTTTgg                     >  2:19349/1‑251 (MQ=255)
                                                                                                                                                                                                                                        gCCTGTTTGACGGGGGTCTTCCCGAAAATGGAAGCCTTCACTGCAGAGTGGACCTTCCAGTTGGCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGACTCTGGCAGATCTGACCCGTTTTTACATAATCGAGAAATGCGCGCAAGCCGCAGGACATATACTGGCGGTTTGGATAATAGa             <  2:127255/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   ggggATCTTCCCGAAAATGGAAGCCTTCACTGCAGAGTGGACCTTCCAGTTGGCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGAATCGTAAATAATTATTGTGGTCGGTGACGCTGACTCTGGCAGATCTGACCCGTTTTTACATAATCGAGAAATGCGCGCAAGCCGCAGGACATATACTGGCGGTTTGGATAAAAGATCTTGAAGccc  >  1:592775/1‑250 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
GGTGGCGATTCTGATTATGCTGCCGACTCTCGGCATTGGTGATATGAACAATCTCTATAAATGGCTGTTGAACCTTAATTCGGTAGTGATGCCGCTGCGTTATTTATGGGTATTTGTTGCATTTATTGCAGTCGTTCGCTTGGCGCAGAAATATAAACCAGAATATGTCTTTATTCGTAATAAGCCTCTGGCAATGACCGTCGGGATTTGGTGTTTTGCCTTTACCGCCTTTGCCTGTTTGACGGGGATCTTCCCGAAAATGGAAGCCTTCACTGCAGAGTGGACCTTCCAGTTGGCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGACTCTGGCAGATCTGACCCGTTTTTACATAATCGAGAAATGCGCGCAAGCCGCAGGACATATACTGGCGGTTTGGATAATAGATCTGGAAGCCT  >  NC_000913/948281‑948774

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