Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,666,4660TG21.6% 205.6 / 58.9 97Y384* (TAT→TAG) clcBH(+)/Cl(‑) exchange transporter
Reads supporting (aligned to +/- strand):  ref base T (41/35);  new base G (11/10);  total (52/45)
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.89e-01

CTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCAT  >  NC_000913/1666368‑1666565
                                                                                                  |                                                                                                   
cTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGTTTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATc                                                                                                    >  1:757708/1‑100 (MQ=255)
cTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATc                                                                                                    >  2:2918024/1‑100 (MQ=255)
cTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATc                                                                                                    >  1:1232630/1‑100 (MQ=255)
 ttttGCTCGGATTGACCGGGATGGCGACACTGTTGGCGGAGACCACGCACGCGACGATTATGTCGACGTTGATGAGAAGTGAAATGACCGGGGAGTATCa                                                                                                   <  2:952364/100‑1 (MQ=255)
      cTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTAc                                                                                              <  2:663220/100‑1 (MQ=255)
      cTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTAc                                                                                              >  1:47819/1‑100 (MQ=255)
           aTTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTcccc                                                                                         <  2:158322/100‑1 (MQ=255)
           aTTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTcccc                                                                                         <  2:816371/100‑1 (MQ=255)
           aTTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTcccc                                                                                         <  1:1948512/100‑1 (MQ=255)
             tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCgg                                                                                       <  1:2361262/100‑1 (MQ=255)
             tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCgg                                                                                       >  1:211790/1‑100 (MQ=255)
             tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                       <  1:295192/100‑1 (MQ=255)
             tGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                       <  1:2542404/100‑1 (MQ=39)
              gACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGt                                                                                      >  1:278026/1‑100 (MQ=255)
              gACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                       <  1:2395836/99‑1 (MQ=255)
              gACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCGGt                                                                                      <  1:1727387/100‑1 (MQ=39)
              gACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCGGt                                                                                      <  1:85000/100‑1 (MQ=39)
                ccGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGttt                                                                                    >  2:2195645/1‑100 (MQ=255)
                   ggATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTatt                                                                                 <  1:781997/100‑1 (MQ=255)
                    gATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattg                                                                                <  1:170427/100‑1 (MQ=255)
                    gATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCGGtcatttc                                                                                <  2:1734389/100‑7 (MQ=37)
                    cttGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCGGtcatttc                                                                                <  1:730548/98‑7 (MQ=37)
                      tGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgat                                                                              >  1:2008600/1‑100 (MQ=255)
                      tGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCGGtcatttcag                                                                              >  2:730544/1‑92 (MQ=37)
                       ggCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgatt                                                                             <  2:2690237/100‑1 (MQ=255)
                        gCGACACTGTTGGCGGCGACCACGCGCGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgattg                                                                            <  2:488062/100‑1 (MQ=255)
                        gCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgattg                                                                            >  1:814609/1‑100 (MQ=255)
                        gCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgattg                                                                            <  2:2505892/100‑1 (MQ=255)
                        gCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgattg                                                                            <  2:2551073/100‑1 (MQ=255)
                            cacTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTg                                                                        <  2:1209565/100‑1 (MQ=255)
                            cacTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTg                                                                        >  2:376914/1‑100 (MQ=255)
                            cacTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTg                                                                        <  2:2643155/100‑1 (MQ=255)
                              cTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCg                                                                      <  1:2255843/100‑1 (MQ=255)
                              cTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCg                                                                      <  2:2624074/100‑1 (MQ=255)
                              cTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCg                                                                      <  1:2151345/100‑1 (MQ=255)
                               tGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGt                                                                     >  1:396232/1‑100 (MQ=255)
                                 ttGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTaa                                                                   >  1:2725878/1‑100 (MQ=255)
                                     cggcgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGc                                                               >  2:2222200/1‑100 (MQ=255)
                                        cgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTc                                                            >  1:199440/1‑100 (MQ=255)
                                           ccACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGt                                                         <  1:1935500/100‑1 (MQ=255)
                                            cacgcacgCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTa                                                        >  1:1681043/1‑100 (MQ=255)
                                            cacgcacgCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTa                                                        >  2:1194688/1‑100 (MQ=255)
                                              cgcacgcGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAAt                                                      >  2:1040903/1‑100 (MQ=255)
                                              cgcacgcGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAAt                                                      <  1:1194513/100‑1 (MQ=255)
                                                  cgcgCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcg                                                  <  1:2797371/100‑1 (MQ=255)
                                                   gcgcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTct                                                  <  2:529773/99‑2 (MQ=255)
                                                    cgcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcgcg                                                >  1:1755017/1‑100 (MQ=255)
                                                    cgcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcgcg                                                <  1:207185/100‑1 (MQ=255)
                                                    cgcCGATTATGTCGACGTTGAGGATACGTGAAATGCCCGGGGAGTATCAGCTACTCCCCGGTTTATTGCTTGCCTGCGTAATTGCGTCGGTAATTTcgcg                                                >  2:874972/1‑100 (MQ=255)
                                                     gcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcgcg                                                <  2:230622/99‑1 (MQ=255)
                                                     gcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCgg                                               <  1:1333213/100‑1 (MQ=255)
                                                             tGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTaca                                       >  1:1905135/1‑100 (MQ=255)
                                                                  aCGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTg                                  >  2:1498256/1‑100 (MQ=255)
                                                                  aCGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTg                                  >  2:920647/1‑100 (MQ=255)
                                                                    gTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGAc                                >  2:2511926/1‑100 (MQ=255)
                                                                    gTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGAc                                <  2:278025/100‑1 (MQ=255)
                                                                      tgatgaTATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGActc                              <  1:1040914/100‑1 (MQ=255)
                                                                       gatgatATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGActct                             >  2:2330000/1‑100 (MQ=255)
                                                                        atgatATGTGAAATGACCGGGGAGTATCAGCTACTCTCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtcta                            >  2:879466/1‑100 (MQ=255)
                                                                        atgatATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtcta                            >  2:879467/1‑100 (MQ=255)
                                                                        atgatATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtcta                            >  2:1884023/1‑100 (MQ=255)
                                                                          gatATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctatc                          >  2:1177862/1‑100 (MQ=255)
                                                                          gatATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctatc                          >  2:2775126/1‑100 (MQ=255)
                                                                             atGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTAc                       >  2:563897/1‑100 (MQ=255)
                                                                              tgtgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTAcc                      <  1:3014152/100‑1 (MQ=255)
                                                                               gtgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCTTGACTCTATCTACCg                     >  2:2836323/1‑100 (MQ=255)
                                                                                 gAAATGACCGGGGAGTATCGGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGcc                   >  1:456627/1‑100 (MQ=255)
                                                                                 gAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGcc                   >  2:2094558/1‑100 (MQ=255)
                                                                                 gAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGcc                   <  1:577620/100‑1 (MQ=255)
                                                                                  aataaaCCGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGcc                   <  1:714901/94‑1 (MQ=38)
                                                                                  aataaaCCGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCCTCGGTAATGTCGCGGACGTTACACCGTGACTCTATCTACCGCCa                  >  2:714898/6‑100 (MQ=37)
                                                                                    aTGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGc                <  1:2626976/100‑1 (MQ=255)
                                                                                    aTGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGc                >  1:2537283/1‑100 (MQ=255)
                                                                                    aTGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTAAACAATGCCTCTATCTCCCGCCAGc                >  1:2575896/1‑100 (MQ=255)
                                                                                     tGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGca               >  1:2454010/1‑100 (MQ=255)
                                                                                     aaaCCGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGca               >  1:2206991/3‑100 (MQ=38)
                                                                                     aaaCCGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGca               >  2:2066035/3‑100 (MQ=38)
                                                                                        ccGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac            >  1:599473/1‑100 (MQ=255)
                                                                                        ccGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac            >  1:487376/1‑100 (MQ=255)
                                                                                        ccGGGGAGTATCAGATTCTCCCCGGTGTAGTGATTGCCTGCGTAATTGCGTCGGGAATGTCGCGGACTTTACACCGTGACTCTATCTACCGCCAGca               <  2:757703/97‑1 (MQ=255)
                                                                                        ccGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGTacac            >  1:243029/1‑100 (MQ=38)
                                                                                        ccGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac            >  1:2638275/1‑100 (MQ=39)
                                                                                         cGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACt           >  1:2879579/1‑100 (MQ=255)
                                                                                         cGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACt           >  1:1340882/1‑100 (MQ=255)
                                                                                          tggcaGTATCATTTTCTCCCCGGTTTATTGATTGCCTCCGTAATTGCTTCGGTAATTTCTCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTg          <  2:2780069/96‑1 (MQ=255)
                                                                                          ggggAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTg          >  1:536801/1‑100 (MQ=255)
                                                                                          ggggAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTg          <  1:1031274/100‑1 (MQ=255)
                                                                                          ggggAGTAGCTGATACTCCGCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACTGTGACTCTATCTACCGCCAGCACACt           >  2:616518/1‑99 (MQ=255)
                                                                                          ggggAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGCCGTTACACCGTGACTCTATCTACCGCCAGCACACt           >  1:2679722/1‑99 (MQ=255)
                                                                                          ggggAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTg          >  1:1044777/1‑100 (MQ=255)
                                                                                             gAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGATGTTACACCGTGACTCTATCTACCGTCAGCACACTgcgc       <  2:1905114/100‑1 (MQ=255)
                                                                                              aGTATCAGCTACTCCCCTGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgca      <  2:1414441/100‑1 (MQ=255)
                                                                                              aGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgca      <  1:2354894/100‑1 (MQ=255)
                                                                                              aGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgca      >  1:817536/1‑100 (MQ=255)
                                                                                               gTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcag     <  2:2332276/100‑1 (MQ=255)
                                                                                                tATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcagc    <  1:2584398/100‑1 (MQ=255)
                                                                                                 aTCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcagca   >  2:3053087/1‑100 (MQ=255)
                                                                                                  tCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCAt  <  1:1151997/100‑1 (MQ=255)
                                                                                                  tCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCAt  <  1:1152957/100‑1 (MQ=255)
                                                                                                  tCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCAt  <  1:2775153/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
CTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCAT  >  NC_000913/1666368‑1666565

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