Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913445,5580AC32.0% 74.9 / 30.9 50V370G (GTA→GGA) yajRputative major facilitator superfamily transport protein YajR
Reads supporting (aligned to +/- strand):  ref base A (24/10);  new base C (0/16);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 1.23e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAA  >  NC_000913/445425‑445690
                                                                                                                                     |                                                                                                                                    
aCGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTTATGGTACTGGCGACTGTCAGCCACACTGCGGCCCGCATTTCGCCAGCGAGAAATAcccccc                                                                                                                                >  2:1170334/1‑139 (MQ=255)
  gCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGa                                                                                                                              <  1:1598197/140‑1 (MQ=255)
     cTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCg                                                                                                                           <  1:946920/140‑1 (MQ=255)
       tttAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTc                                                                                                                         <  2:1673180/140‑1 (MQ=255)
            cGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACa                                                                                                                    <  1:895915/140‑1 (MQ=255)
            cGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGGTATTTTATGGTACATGCGACTGTCAGCCACACCGCCACCAGCCTTCCGCCCGCGAAAAACACCCCCCCACCCGCAAACa                                                                                                                    >  2:577874/1‑140 (MQ=255)
            cGCCTCGTTTGCGGCAATGTCCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCGCACTGCGGCCAGCATTGCGCCAGCGAGAAATCCCCCCTGACCGTCAAACa                                                                                                                    <  1:1248691/140‑1 (MQ=255)
             gCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGTCAGCGAGAAATACCCCCTGACCGTCAAACAt                                                                                                                   <  1:1461170/140‑1 (MQ=255)
              ccTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCCTTGCGCCAGCGAGAAATACCCCCCGACCGTCAAACATg                                                                                                                  >  1:960026/1‑140 (MQ=255)
                tCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAAAACCCCCTGACCGTCAAACATGcc                                                                                                                >  2:83965/1‑140 (MQ=255)
                tCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGGTCTTTCATGGGACCGGCGGCCGTCAGCCACACTGCCGGCCGCCTTTCGCCAGCGAGGAAAAACCCCCGACCCGCCAACAAGcc                                                                                                                >  2:337005/1‑140 (MQ=255)
                     tGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTaa                                                                                                           <  2:826423/140‑1 (MQ=255)
                       cGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATc                                                                                                         <  1:1321673/140‑1 (MQ=255)
                          cAATGTTCGCTGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTGCTGGCGACTGACAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAg                                                                                                      <  2:1242680/140‑1 (MQ=255)
                          cAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAg                                                                                                      <  2:163991/140‑1 (MQ=255)
                            aTGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAgcc                                                                                                    <  1:83965/140‑1 (MQ=255)
                                 cGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTTTTTTATGGGACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCCGCGAGAAATACCCCCCGACCGTCAAACAAGCCGTTTATCCAGCCGcccc                                                                                               >  2:979701/1‑139 (MQ=255)
                                 cGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGGACTGGCGACTGTCAGCCACACCGCGGCCCGCATTGCGCCAGCGGGAAATACCCCCCGACCGTCAAACATGCCGCTAATCCAGCCGcccc                                                                                               >  2:171166/1‑139 (MQ=255)
                                                 cgcAAATTGTTGACATGGGGCGTTTTTTTCATGGTACTGCCGACTTTCCGCCACCCCGGGGCCCGCATGCCCCCAGCGAGAAATACCCCCTGACCGCCAACCATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATcgc                                                                               <  2:328232/140‑1 (MQ=255)
                                                   cAAACTGCTGCCATCCGGCGTTTTTTTCATGGTACGGGCGACTGTCGCCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATcgcca                                                                             <  2:655740/140‑1 (MQ=255)
                                                          ctgACATACGGCGTTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAg                                                                      <  1:1482280/140‑1 (MQ=255)
                                                                  cggcggGTCTTTCCTGGTACAGGCGAGCGTCAGCCACACCTCCGCCATCAGTGCGTCCACCTGAAATATCCCCCGCCCCCCAAAAATGCCGTTAATCCAGCCGCCCAGCGTACCGCCCATCCCCACGCCGATaaaaacgg                                                              >  2:215743/1‑135 (MQ=255)
                                                                     ggtttCTTTCAGGGTACTGGCGACTGTCAGCCACACGGCGGCCAGCATTGCGCCAGCGAAAAATCCCCCCTGACCTTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAActggctgg                                                           <  1:1096708/137‑1 (MQ=255)
                                                                      gtttCTTTCATGGTACGGGGGACTGCCGGCCACACTGCGGCCAGCTTTGCGCCAGCAAAAAATCCCCCCTGCCCGTCAAACATGCCGTTAACCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCtggt                                                          <  2:129256/138‑1 (MQ=255)
                                                                       tttttttCATGGTACGGGCGACTGTCAGCCACACGGCGGCAAGCATTGCGCCAGCGAGAAATCCCCCTGGCCCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCtggtg                                                         <  2:1673085/136‑1 (MQ=255)
                                                                       gTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCtggtg                                                         <  1:1408258/140‑1 (MQ=255)
                                                                            ttaagggTACTGGCGACTGTCAGCCACACTGGGGCCAGCATTGCGCCGGCGAGAAATCCCCCCTGACCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTa                                                    <  1:1762148/135‑1 (MQ=255)
                                                                                  gtccgggCAACTTTCACCCACACGGCGGCCAGCATTGCCCCAGCGAAAAAACCCCCCTACCCGCAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACAcc                                              <  2:1518244/135‑1 (MQ=255)
                                                                                   tACGGGCGACTGCCCGCCCCACGGCGGCCACCTTTGCGCGAGGGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACAccc                                             <  2:156579/140‑1 (MQ=255)
                                                                                    ccgggCGACTGTCAGCCACACGGCGGCCAGCATTGCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAACCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCa                                            <  1:337005/137‑1 (MQ=255)
                                                                                     cTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCAt                                           >  1:840844/1‑140 (MQ=255)
                                                                                        gCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGc                                        >  2:1051890/1‑140 (MQ=255)
                                                                                            ctgctaGCCCCACAGCGGCGAGCATCGCAAGAGGGAAAAATGCGCGCTGCCCACAAAACATCGCGTTAATCCAGCCGCCCAGGGAGCCGCCAATCGCCCCGCCAAGAAACTGGCTGGTGGAGTAAACGCCCATCGCCGTa                                    <  2:1529992/135‑1 (MQ=255)
                                                                                            cTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTa                                    >  2:580078/1‑140 (MQ=255)
                                                                                            cTGCCAGCCCCACTGCGGCCAGCATTGCGCCAGGGAAAAATCCCCCCGGACCGTAAAACATGCCGTTAACCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTa                                    <  2:1268813/140‑1 (MQ=255)
                                                                                             tGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTAc                                   >  2:917253/1‑140 (MQ=255)
                                                                                             tGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGCGGAGTAAACACCCATCGCCGTAc                                   >  2:110321/1‑140 (MQ=255)
                                                                                               tCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCt                                 >  1:13772/1‑140 (MQ=255)
                                                                                               tCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCtt                                >  1:917856/1‑140 (MQ=255)
                                                                                                  gCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTg                              >  2:868082/1‑140 (MQ=255)
                                                                                                       acTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAAcc                         >  2:1092652/1‑140 (MQ=255)
                                                                                                         tGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACctg                       >  2:1687852/1‑140 (MQ=255)
                                                                                                         ggCGGCCAGCTTTGCGCCAGCGAAAAATCCCCCCTGCCCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACctg                       <  2:960026/139‑1 (MQ=255)
                                                                                                         ggCGGCCAGCATTGCGCCAGCGAGAAATCCCCCCTGCCCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACctg                       <  2:509687/139‑1 (MQ=255)
                                                                                                           cGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACctgct                     >  1:1543523/1‑140 (MQ=255)
                                                                                                           cGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACctgct                     >  2:295826/1‑140 (MQ=255)
                                                                                                             gCCAGCATTGCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTgg                   <  1:1092652/140‑1 (MQ=255)
                                                                                                             gCCAGCATTGCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTgg                   <  1:1051890/140‑1 (MQ=255)
                                                                                                             gCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTgg                   >  2:226478/1‑140 (MQ=255)
                                                                                                             gCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTgg                   >  1:1177671/1‑140 (MQ=255)
                                                                                                             gCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTgg                   >  1:205391/1‑140 (MQ=255)
                                                                                                               cAGAATTGCGCCAGCGAAAAATCCCCCTGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGa                <  2:917856/140‑1 (MQ=255)
                                                                                                                aGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGa                >  1:823975/1‑140 (MQ=255)
                                                                                                                    ttGCGCCAGCGAGAAATCCCCCCTGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGActct            <  1:908086/140‑1 (MQ=255)
                                                                                                                    ttGCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGActct            <  1:29917/140‑1 (MQ=255)
                                                                                                                     tGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCtt           >  1:785399/1‑140 (MQ=255)
                                                                                                                      gcgcCAGGGAGAAATCCCCCTTGACCGCCAAACATGCAGTTAATCAGGCGGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCttt          <  2:1462623/140‑1 (MQ=255)
                                                                                                                      gccccAGCGAAAAATCCCCCCTGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCttt          <  2:13772/137‑1 (MQ=255)
                                                                                                                       cgcCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTa         >  1:306285/1‑140 (MQ=255)
                                                                                                                       cgcCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTa         >  1:355200/1‑140 (MQ=255)
                                                                                                                           aGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGa     <  1:295826/140‑1 (MQ=255)
                                                                                                                           aGCGAGAAATCCCCCCTGACCGCAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGa     <  1:1529425/140‑1 (MQ=255)
                                                                                                                            gCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGAt    >  1:1646753/1‑140 (MQ=255)
                                                                                                                            gCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGAt    >  2:879256/1‑140 (MQ=255)
                                                                                                                            gCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGAt    >  1:11737/1‑140 (MQ=255)
                                                                                                                            gCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCACCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGAt    >  2:750301/1‑140 (MQ=255)
                                                                                                                              gagaAATCCCCCCTGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATaa  <  1:171166/140‑1 (MQ=255)
                                                                                                                              aaaaaaTCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATaa  <  1:979701/137‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                    
ACGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAA  >  NC_000913/445425‑445690

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