Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,921,1730TC34.2% 64.4 / 12.4 42I310V (ATT→GTT) gudPgalactarate/D‑glucarate transporter GudP
Reads supporting (aligned to +/- strand):  ref base T (24/3);  new base C (0/14);  total (24/18)
Fisher's exact test for biased strand distribution p-value = 1.93e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CATGTTGAACAGGCCACCGGAAAGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGC‑AGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAA‑TAATCCCACCC‑AGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGT  >  NC_000913/2920928‑2921407
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                           
cATGTTGAACAGGCCACCGGAAAGACCGCTGATCTCTTTTGGCGCGGTATCCGCCCTTACTGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGC‑AGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAA‑TAATcccac                                                                                                                                                                                                                                        >  1:410361/1‑251 (MQ=255)
             ccACCGGAAAGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGC‑AGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAA‑TAATCC‑ACCC‑AGCACAccgccg                                                                                                                                                                                                                         <  2:376084/251‑1 (MQ=255)
                    aaaGACCGATGCTCTCTTTTGGCGCGAAATCTGCAATTACTGCCCAAACCAGCGCACCGCTGCCCTTACCGAAGACGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGC‑AGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAA‑TAATCCCCCCC‑AGCACACCGCCGATaaa                                                                                                                                                                                                                    >  2:253334/1‑250 (MQ=255)
                       gACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGC‑AGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAA‑TAATCCCACCC‑AGCACACCGCCGATAaaaccg                                                                                                                                                                                                                >  2:343508/1‑251 (MQ=255)
                          cGCTGATCTCTTTGGGCGCGGTATCTGCCATTCCTGCCCACCCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACTTTGAGGTAGTGGC‑AAAACCCCATCACCATCGACAGCAACATGCCGATTACTATCGGTGTTGTACGCTCAATGTTCGGCGACCCCGTGGGGCGCAGCAGCCAGTGGAAAA‑TAAGCCCCCCC‑AGCACACCGCCGATAaaaccgcaa                                                                                                                                                                                                             <  1:238239/251‑1 (MQ=255)
                              gATCTCTTTTGGCGCGTTATCTGCCATTACTCCCCACCCCACCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCACCCACCCAACGTTGACGTATTTGC‑AGAACCCCACCCCCACCGACAGCAACATCCCCATAACGATCGGTGTTTTACGCGAAATGTCCAGCGACCCCGTGCGGCGCATCACCCAGTCGGAAA‑TAACCCCACCC‑AGCACACCGCCGATAaaaccgcaaacc                                                                                                                                                                                                          <  1:93377/250‑1 (MQ=255)
                                               tATCTGCCATTACTGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGC‑AGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAA‑TAATCCCACCC‑AGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCcac                                                                                                                                                                                         >  1:298533/1‑250 (MQ=255)
                                                          aCTGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGC‑AGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAA‑TAATCCCACCC‑AGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCttt                                                                                                                                                                             >  1:76950/1‑251 (MQ=255)
                                                                  ccccgccccccctaccccctccccGAAAAAGGCAAGCCCCACAAAGCCAATGATCATCACCACAACGTTGACGGAGTTGC‑AGAACCCCACCACCATCAACAGCACCACGCCCTTTCCTATGGGTGTTTTCCGCGAAATGTCCAGCGCCCCCGGGCGGCCCACCCGCCAGTCGAAAA‑TAACCCCACCC‑ACCCCCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCg                                                                                                                                                                     <  2:214497/229‑1 (MQ=255)
                                                                         gcCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGC‑AGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAA‑TAATcccaccc                                                                                                                                                                                                                                      >  1:247594/1‑180 (MQ=255)
                                                                         gcCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTGGACGTAGTTGC‑AGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTGACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAA‑TAATcccaccc                                                                                                                                                                                                                                      <  2:247594/180‑1 (MQ=255)
                                                                           ccagATGCCCTAACCGAAGAAGGCCAGCCCAAAAAAGCCAATGATCACCCACTCAACGTTGACG‑AATTGCCAGAACCCCATCCCCATCGACACCAACATGCCCATTACGATGGGTGTTTTCCGCGCAATGTTCAGCATTCCCGGGCGGCGCATCA‑CCAGTTGGAAAATACCCCCCCCC‑AGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCcgc                                                                                                                                                            <  1:309274/248‑1 (MQ=255)
                                                                                tGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGC‑AGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAA‑TAATCCCACCC‑AGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGcc                                                                                                                                                       >  1:399613/1‑251 (MQ=255)
                                                                                       aCCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGCCGTAGTTGC‑AGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAA‑TAATCCCACCC‑AGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCAcc                                                                                                                                                 >  1:380476/1‑250 (MQ=255)
                                                                                           aagaagGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGC‑AGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAA‑TAATCCCACCC‑AGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGAt                                                                                                                                             >  2:30785/1‑250 (MQ=255)
                                                                                               agGCCAGCCCCATAAAGCCGATGATCATCCCCTCAACGTGGACGTAGTTGC‑AGAACACCATCACCATCGACAGCAACATCCCCATTACGATCGGTGTTTTACGCGAAATGTTCAGCGACCCCGGGCGGCGCACCAGCCAGCCGAAAA‑TAACCCCCCCC‑ACCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAAcc                                                                                                                                        <  1:106491/251‑1 (MQ=255)
                                                                                                          aTAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGC‑AGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAA‑TAATCCCACCC‑AGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGt                                                                                                                             >  2:189461/1‑251 (MQ=255)
                                                                                                                  aatgatCACCCACTCACCGTTGACGTAGTTGC‑AAACCACCATCACCACCGACAGCAACAGGCCCATTACGATCGGTTTTTTACGCGCAATGTTCAGCGATCCCGTGGGGCGCATCAGCCAGTCGGAAA‑AAACCCCCCCC‑ACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAaaga                                                                                                                     <  1:256434/250‑1 (MQ=255)
                                                                                                                      accaccccctaaccgtggACGTAGTTGC‑AAAACCCAATCCCCATCGACAGCAACAGGCCCTTGACGTTCGGTGTTTGACCCGCAAGGTTCAGCGAGCCGGGGCGGCGACTCACCAACTCGAAAA‑AAACCCCCCCC‑ACCAACCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGaagta                                                                                                                 <  1:109185/234‑1 (MQ=255)
                                                                                                                                      tGACGTAGTTGC‑AGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAA‑TAATCCCACCC‑AGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGAt                                                                                                  >  1:123977/1‑250 (MQ=255)
                                                                                                                                        aCGTAGTTGC‑AAAACCCCATCCCCACCGACAGCAACATGCCCATTACGACCGGTGTTTTACGCCCAATGTCCAGCGACCCCGTGCGGCGCATCAGCCAGCCGAAAA‑AAACCCCACCC‑ACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACa                                                                                               <  1:30785/251‑1 (MQ=255)
                                                                                                                                                    aaaCACCATCACCATCGAAAGAAACATCCCAATTCCGATCGGTGTTTTACGCGAAATGTTCAGCGCCCCCGTGCGGCGCACCACCCAGTCGAAAA‑TAACCCCACCC‑AGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCGa                                                                                    <  1:379045/250‑1 (MQ=255)
                                                                                                                                                         ccatcaccatcGACAGCAACATGCCCATTACGATTGGTGTTTTACGCGCAATGTTCAGAGATCCCGTGCGGCGCATCAGCCAGTCGGAAA‑TAATCCCACCC‑AGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATACCGCGTGCCTGCACCAGATAAACCGGGAACAAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCGAgata                                                                                >  1:4501/1‑250 (MQ=255)
                                                                                                                                                                catcGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAA‑TAATCCCACCC‑AGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCGATATAAACGCCg                                                                         >  2:160724/1‑250 (MQ=255)
                                                                                                                                                                  ccgacataaacttGCCCAGTACGAGCGGTTTTTTCCGCGCATTGTCCAGCGATCCCTTGCGGCGCAACACCAAGTCGAAAA‑AAACCCCCCCC‑ACAACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCGATATAAACGCCGatc                                                                      <  1:48494/239‑1 (MQ=255)
                                                                                                                                                                     aCAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAA‑TAATCCCACCC‑AGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCGATATAAACGCCGatcatc                                                                   >  1:77492/1‑251 (MQ=255)
                                                                                                                                                                                      aCGATCGGTATTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAA‑TAATACCACCC‑AGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTACTAAAGAAGTCAGTCAGGGCGTTGATACAGTACTGACCGATATAAACGCCACTCATCATCCGTGACCCTagc                                                    >  2:210740/1‑249 (MQ=255)
                                                                                                                                                                                         aTCGGTGTTTTACGCGAAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAA‑TAATCCCACCC‑AGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCACGCGTGCCTGCAACAGATAAAACGGGAACCAGGAAATAAAGAAGAAAGTCAGGGCGTTGATACAGTACTGACCGATATAAAAGCAGATCATCATCCGTGACCCTAGCAGCTg                                               >  2:231520/1‑251 (MQ=255)
                                                                                                                                                                                          ggggTGTTTTCCGCGCAATGTCCAGCGACCCCGTGCGGCGCACCAGCCAGTCGAAAA‑AAACCCCCCCC‑AGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGt                                              <  1:343508/249‑1 (MQ=255)
                                                                                                                                                                                           cGGTGTTTTGCGCGCAATTTCCAGCGACCCCGTGCGGCCCAACACCCAGTCG‑AAA‑AAATCCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGt                                              <  1:148356/250‑1 (MQ=255)
                                                                                                                                                                                           cGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAA‑TAATCCCACCC‑AGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGtt                                             >  1:330600/1‑251 (MQ=255)
                                                                                                                                                                                                   tACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAA‑TAATCCCACCC‑AGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGc                                     >  1:47397/1‑251 (MQ=255)
                                                                                                                                                                                                       cgcAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAA‑TAATCCCACCC‑AGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGccccc                                 >  1:248036/1‑251 (MQ=255)
                                                                                                                                                                                                         cAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAA‑TAATCCCACCC‑AGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCa                                                                                                                                >  2:29711/1‑154 (MQ=255)
                                                                                                                                                                                                             gTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAA‑TAATCCCACCC‑AGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTc                            >  1:97395/1‑250 (MQ=255)
                                                                                                                                                                                                                  gCGATCCCGTGCGCCGCATCAGCCAGTCGAAAA‑AAACCCCCCCC‑AGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTg                      <  1:377410/251‑1 (MQ=255)
                                                                                                                                                                                                                     cccccGTGCGGCGCACCAGCCAGTCGAAAA‑AAACCCCCCCC‑ACCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAAc                   <  1:236686/249‑1 (MQ=255)
                                                                                                                                                                                                                              ggcgCATCAGCCAGTCGGAAA‑TAATCCCACCC‑AGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACttt           >  2:314279/1‑250 (MQ=255)
                                                                                                                                                                                                                                cgcACCAGCCATCCGGAAA‑TAACCCCACCC‑AGAGCGCCGCCGATAAATCCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAAGCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCATGTAATAAAGAAGTAAGTCAACGCGTTGATACAGTACTGACCGATATAAACGCCTATCATCATGCGTGACCCTATCAGCTGTTTGATCTGCCTCCACTTCACTCTGAACGGCACTTTAAc        <  1:136136/251‑1 (MQ=255)
                                                                                                                                                                                                                                cacaacAGCCCGTAGGAAA‑AAACCCCCCCC‑AGCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACAACCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAAc        <  2:307038/246‑1 (MQ=255)
                                                                                                                                                                                                                                 gcATCAGCCAGTCGGAAA‑TAATCCCACCC‑AGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACt       >  2:396633/1‑251 (MQ=255)
                                                                                                                                                                                                                                      aGCCAGTCGGAAA‑TAATCCCACCC‑AGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTggg  >  2:101659/1‑250 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                           
CATGTTGAACAGGCCACCGGAAAGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGC‑AGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAA‑TAATCCCACCC‑AGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCGATATAAACGCCGATCATCATCCGTGACCCTAGCAGCTGTTTGATCTGCCCCCACTTCACGCTGAACGGCACTTTAACTTTGGT  >  NC_000913/2920928‑2921407

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