Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,184,542 C→A P1100Q (CCG→CAG)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,184,5420CA97.2% 207.9 / ‑6.6 72P1100Q (CCG→CAG) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base A (26/43);  total (26/46)
Fisher's exact test for biased strand distribution p-value = 5.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCA  >  NC_000913/4184403‑4184676
                                                                                                                                            |                                                                                                                                      
cTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCt                                                                                                                                     <  1:804173/143‑1 (MQ=255)
  ggCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGAGCAACCCGATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTgg                                                                                                                                   >  1:295134/1‑143 (MQ=255)
   gCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTggg                                                                                                                                  >  1:773066/1‑143 (MQ=255)
    cGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGc                                                                                                                                 <  2:1164900/143‑1 (MQ=255)
    cGGCTAAACGCCGAATCCAGCCTGGTGACACGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCTATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGc                                                                                                                                 <  1:366939/143‑1 (MQ=255)
       ttAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTa                                                                                                                              >  1:289931/1‑143 (MQ=255)
        tAAACTCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATAGTACTGAACCAGCTGGGCGTAc                                                                                                                             >  2:784578/1‑143 (MQ=255)
        tAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTAc                                                                                                                             <  1:827668/143‑1 (MQ=255)
         aaaCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATGTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGACAACGG‑CACGCCGGTCGACCTCGTACTGAACCAGCTGCGCCTAcc                                                                                                                            >  2:27065/1‑143 (MQ=255)
           aCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGt                                                                                                                          >  1:713046/1‑143 (MQ=255)
                  aTCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTAt                                                                                                                   <  1:720752/143‑1 (MQ=255)
                   tCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGAGATGCCTTACGATGAAAACGG‑GACGCCGGTAGACATCGTCCTGAACCAGCTGGGCGTGCCGTCTCGTATg                                                                                                                  >  2:936846/1‑143 (MQ=255)
                     cAGCCCGGTGACAAGAGGGCAGTTCGGCACGGTAACCAGGGTGCAAGTTCTAAGATCAACCCGATCGAAGATATGCCTGACGATGAAAACGG‑TACGCCGGTAGACAGCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGaa                                                                                                                <  1:23727/143‑1 (MQ=255)
                      agctcggTGACAAGTTGGCACGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAc                                                                                                               <  2:631433/138‑1 (MQ=255)
                        ccTGGTGACAAGATGGCAGGTCGTCACGGAAACAAGGGTGTACTTATCAACGCGGCCCCGATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACAt                                                                                                             <  1:855602/143‑1 (MQ=255)
                         cTGGTGACAAGATGGCAGGTCGTCACGGGAACAAGGGTGTAACTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATc                                                                                                            <  1:54797/143‑1 (MQ=255)
                            gTGACAAGAGGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATgctcatttgg                                                                                                         <  2:784451/143‑10 (MQ=255)
                             tGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTc                                                                                                        >  2:159080/1‑143 (MQ=255)
                              gACAAGATGGCAGGTCGTCACGGGAACAAGGGTGTAATCTCTACGTTCCACCCCATCGAAGATGCGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCa                                                                                                       <  1:44359/143‑1 (MQ=255)
                               aCAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCATCTGGGCGTACCGTCTCGTACGAACATCGGTCAg                                                                                                      >  2:348511/1‑143 (MQ=255)
                                  aGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATc                                                                                                   <  1:368610/143‑1 (MQ=255)
                                     tGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTc                                                                                                <  1:1099082/143‑1 (MQ=255)
                                     tGGCAGGTCGTCACGGTAACAAGGGGGTGATTTCGAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTc                                                                                                <  2:468149/143‑1 (MQ=255)
                                        cAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCATCCCGATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaa                                                                                             <  1:804609/143‑1 (MQ=255)
                                         agggcgttgcGGTAACAGGTATGTAATTTCTAAGGTCAAGCCGATCGAAGATATGCCTAACGAT‑ACAGCGGCCACCGAGATAGACACCGAACTGAACAAGCTTGGCGTACCGTCTGGTATGAAGATCGGTCAGATCCTCGaaa                                                                                            <  2:759678/134‑1 (MQ=255)
                                         aGGTCGTCCGGGTAACCAGGGTGTAATTTCTAAGATCAACCCGATCGAAGGTATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaaa                                                                                            <  1:1128057/143‑1 (MQ=255)
                                         aGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACAGCGTTCAGATCCTCGaac                                                                                            >  1:484027/1‑142 (MQ=255)
                                         aGGGCGTCACCGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaaa                                                                                            <  2:968233/143‑1 (MQ=255)
                                             cgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCa                                                                                        <  1:1091086/143‑1 (MQ=255)
                                             cgtcACGGCACCAAGGGTGTAGTGTCTACGATCAACCGGATCGAATATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCa                                                                                        <  1:517362/143‑1 (MQ=255)
                                              gtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGACACAGCTGGGCGTACCGTCTCGTATGCACATCGGTCAGCTCCCCGGAAcaccc                                                                                       >  2:665856/1‑139 (MQ=255)
                                               tcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGCTAAAGCATCAGTTTGACACGG‑TACGCAGTCACACATCGTACGGAACAACCGCGGAGAACCGACTAGTAGGAACATCAAACAGATCCTAGAAACCCaac                                                                                      >  2:141162/1‑141 (MQ=255)
                                               tcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTAGGAACATCGGTCAGATCCTCGAAACCCAcc                                                                                      >  1:196508/1‑143 (MQ=255)
                                                     tAACAGGGGTGTCATTTCTAAGATCAACCCGGTCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTa                                                                                <  2:673124/143‑1 (MQ=255)
                                                      aaCAAGGGTGTAAGTTCTAAGAGCAACCCGATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTAt                                                                               <  2:361702/143‑1 (MQ=255)
                                                      aaCAAGGGCGGAATTCCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGCACCAGCTGGGCGTACCGTGTCGTATGAAGATCGGTCAGATCCTCGAAACCCACCTGGGTAt                                                                               <  1:839117/143‑1 (MQ=255)
                                                        cAAGCGTGTCATTTCTAAGCTCAACCCGATCGAAGATATGCCTTACGACGAAAACGG‑TACGCCGGTAGACATCGTACTGCTCCCCCTGGGGGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATgg                                                                             <  1:674930/143‑1 (MQ=255)
                                                         aaGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGCACATCGGTCAGATCCGCGAAACCCACCCGGGCACGGc                                                                            >  2:423306/1‑143 (MQ=255)
                                                           gggTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTg                                                                          >  2:1052655/1‑143 (MQ=255)
                                                            ggTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGc                                                                         <  1:702260/143‑1 (MQ=255)
                                                               gtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaa                                                                      >  2:1102644/1‑143 (MQ=255)
                                                                 aaTGTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTGCTGAACCGTCTGGGCGTCCCCACACGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAg                                                                    <  2:834697/143‑1 (MQ=255)
                                                                  aTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCACGAAACCCACCTGGGTATGGCTGCGAAAgg                                                                   >  2:857060/1‑143 (MQ=255)
                                                                   tttCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTGTCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGt                                                                  <  1:321476/143‑1 (MQ=255)
                                                                                cccGATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCGACCTGGGTATGGCTGCGAAAGGGATCGGCGagcaga                                                     >  1:78261/1‑138 (MQ=255)
                                                                                 ccGCTCGAAGATATGCCTTACGAGGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGAt                                                    <  1:507500/143‑1 (MQ=255)
                                                                                  gtaTCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATc                                                   <  2:1023480/141‑1 (MQ=255)
                                                                                  ctaTCGATGATATGCCCTACGATGATAAGGC‑TACGGCCGTAGACATCGGACTGAACCAGCTGGGCGTACCGTCTAGTATGCACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATc                                                   <  1:751526/141‑1 (MQ=255)
                                                                                  cGATCGAAGACATGCCGTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATc                                                   <  2:318437/143‑1 (MQ=255)
                                                                                   gATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGGATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCa                                                  <  2:773066/143‑1 (MQ=255)
                                                                                     tCGAAGATAGCGCTTACGCTGAACCCGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGGACCGTCTCGTATGAACATCGGTCAGACCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAAc                                                <  1:565559/143‑1 (MQ=255)
                                                                                        aaGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACCTCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACAAGGGGACGGGATGGAAAGGTATCCGCGACAAGAGCAACGcc                                             >  2:514929/1‑143 (MQ=255)
                                                                                        aaGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGcc                                             >  1:923791/1‑143 (MQ=255)
                                                                                            tatGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGc                                         <  1:401738/143‑1 (MQ=255)
                                                                                            tatGCCTTACGAGGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGc                                         <  1:256037/143‑1 (MQ=255)
                                                                                                 cTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGGCGGGGATCTGGGAGAAGATCAACGCCATGCTGaaa                                    >  1:910598/1‑143 (MQ=255)
                                                                                                    aCGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTCTGGCTGCGAAAGGTAACGGCGACAAGATCAACGCCATGCTGAACcag                                 >  2:542729/1‑143 (MQ=255)
                                                                                                    aCGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcag                                 <  2:295134/143‑1 (MQ=255)
                                                                                                     ggagcaagaCGG‑TACGGCGGTAGACATCGGACTGCACCAGCTGGGCGTACCGTCTCGTATGAGCATCGGTCAGCGCCTCGAAACCGACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagc                                <  2:1115434/135‑1 (MQ=255)
                                                                                                     cGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGCCCCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGCGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagc                                <  1:990510/143‑1 (MQ=255)
                                                                                                      gATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagca                               <  1:1102644/143‑1 (MQ=255)
                                                                                                      gATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGGATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagca                               >  1:875088/1‑143 (MQ=255)
                                                                                                           aaaCGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACAGTCTCGTATGAACATCGGTCAGATCCTCGAGACCGACCTGGGTATGGCTGCGAAAGGGGTCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCcag                          >  1:455853/1‑143 (MQ=255)
                                                                                                            aaCGG‑TACGCCACTAGACCTCGTACTGAACCATCTGGCCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaaga                         <  2:760432/143‑1 (MQ=255)
                                                                                                             ccGG‑TACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAGCGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaagaa                        <  2:645080/142‑1 (MQ=255)
                                                                                                             aCGG‑TATGCTGTTAGACTTCGTACTGATCCAGCTGGGCGTACCGTCTCGTATGAACATCTGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaagaa                        <  2:1048340/143‑1 (MQ=255)
                                                                                                                       cGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGaag               <  1:286599/143‑2 (MQ=255)
                                                                                                                         gTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACt             <  1:421600/143‑1 (MQ=255)
                                                                                                                             aCATCGTACTGATCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTgcgc         <  1:907682/143‑1 (MQ=255)
                                                                                                                             aCATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGGCAGATCCTCGAAACCCACCTGGGGGTGGCTGCGAAAGGTATCGGCGACAAGCTCAACGCCATGCGGAACAAGCAGCAAGAAGGAGCGAAACTgcgc         >  1:381134/1‑143 (MQ=255)
                                                                                                                               aTCGTACTGAACCAGCTGGGCGTACGGTCTCGTAGTAGCGGCGGGCAGATCCTCTAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGACACTGCGCGa       >  1:551327/1‑143 (MQ=255)
                                                                                                                                    aCTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCa  <  1:84693/143‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                      
CTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGG‑TACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCA  >  NC_000913/4184403‑4184676

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