Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,132,787 A→C I204S (ATC→AGC)  wcaA ← putative glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,132,7870AC98.4% 192.2 / NA 63I204S (ATC→AGC) wcaAputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (30/32);  minor base G (1/0);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTT  >  NC_000913/2132648‑2132926
                                                                                                                                           |                                                                                                                                           
gATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGTGGTATAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGctct                                                                                                                                          <  1:336252/143‑1 (MQ=255)
 aTCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTg                                                                                                                                         >  1:609316/1‑143 (MQ=255)
  tCTGGTAGGGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCTGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGc                                                                                                                                        <  2:417610/143‑1 (MQ=255)
     ggTAGAGGTCAGACAGCGGATCTTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTGCCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTc                                                                                                                                     <  2:566122/143‑1 (MQ=255)
     ggTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAACGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTc                                                                                                                                     >  1:910452/1‑143 (MQ=255)
      gTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCCGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCg                                                                                                                                    >  1:919829/1‑143 (MQ=255)
      gTAGAGGGTAAACAGCTGATATTTTTTGCAGCAGGGTACGAATGTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTGCACCGTGATGGATAGGCAGGCTCTGCGTCg                                                                                                                                    >  2:913694/1‑143 (MQ=255)
         gaccgTAAACACCTGATAGTTTTTGCTGGCGCGGTCCAATTTGTGATTGGGTTTGGGGTGCAATCGAAAATACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCct                                                                                                                                 <  1:1108024/139‑1 (MQ=255)
             gTAAACAGCTGATATTTTTTGCTGGCGCGGTCGCATTTGTCTTTGTGCTTGCGGTAAAGGTGGAAGTACCCGGAGAGTTTTTTCGGTGACGAGGGGATCTGCATTTCCCCGTGATTGGTATGCCGGGCCTGCGtggtcgccgt                                                                                                                             >  2:194038/1‑134 (MQ=255)
             gTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtcttc                                                                                                                             <  1:853759/143‑1 (MQ=255)
              tAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCCCCGCGTCCCCtcttcc                                                                                                                            >  2:356076/1‑142 (MQ=255)
                  cAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGGTTGATATGCAGGCTCTGCGTCGCCTCTTCCAcgt                                                                                                                        >  1:1018159/1‑141 (MQ=255)
                      tGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGAGGAGGTAATCTGCATCTGACCGTGATTTATATTAAGGCGCGGCGTCGCCtcttgcactgccca                                                                                                                    >  2:477099/1‑133 (MQ=255)
                      tGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCa                                                                                                                    >  2:1072241/1‑143 (MQ=255)
                       gATATGTTCGGCTGGCGCGGTCGACTTTGTCTTTGTGCTTGCGGTATAAGTGGAAGTTCCCCGTGAATTTTTTCGTTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACttcacac                                                                                                                   <  2:925959/143‑6 (MQ=255)
                        cgatTTGTTTCTGGAGCGGTAGAACTTGACGGTTTGCTGGGGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACg                                                                                                                  <  1:201797/141‑1 (MQ=255)
                        atatTTTGTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTATTTTCCACg                                                                                                                  <  2:1039317/143‑1 (MQ=255)
                           tttttgGCGGGGGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGtt                                                                                                               <  2:704371/143‑1 (MQ=255)
                              gttGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGc                                                                                                            <  1:401989/142‑1 (MQ=255)
                                tGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCg                                                                                                          <  1:293971/143‑1 (MQ=255)
                                      cgcgGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTTTTTTGGACGCACCGCCGTACTcc                                                                                                    >  2:959650/1‑143 (MQ=255)
                                         ggTCGATTTTGTCTTTGTGCTTGCGGTAATATTGGTATTTCTCGTAGAATTTTTTCGGTGACGAGGTACTCTGCAGGGCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccacc                                                                                                 <  1:833134/143‑1 (MQ=255)
                                         ggTCGAATTTGTCTTTGTGCTTACGCTAAATGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCCGGCTCTGCGTCACCCCATCTACTTTCCACGGTTCGACCTACCcaacc                                                                                                 >  2:137084/1‑143 (MQ=255)
                                          gTCGTATTGGTCTTTGTGCTTGCGGTAAAAGTGGAAGTTCCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccacca                                                                                                <  2:106252/143‑1 (MQ=255)
                                               agtggccTTGGTGCCTGTGGCAATAGTAGAACGCCCCGCAGTATTTTTTCGGGGACGAGGTGATCTCCATTTCTCCGTGATTTCTATGCAGGCTCTGCGTCGCCCCTTCTACTTTCCTCGGTTCGCCGTACTCCACCACCATc                                                                                           <  2:99529/137‑1 (MQ=255)
                                                    tCTTTGTGCTTGCGGTAAAGGGGGGAGTACGCGGGGGTCTATTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAg                                                                                      <  2:228068/143‑1 (MQ=255)
                                                       ttGTGCGTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTGACCGAGATTAACATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGaa                                                                                   <  2:971126/143‑1 (MQ=255)
                                                           ggtggCGGTGAAAGGGGCAGTAGCTGGAGAATTTTTTCGGTGACGAGGTGATCTGCAGGTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTGGCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATg                                                                               <  2:973645/139‑1 (MQ=255)
                                                            cTTGCGGCACAGGGGGGAGTGCCCGGAGAATTTTTTCGGTGCCGAGGTGATCTGCATTTCACCGTGATCGATATTCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGt                                                                              <  2:1107480/143‑1 (MQ=255)
                                                              tGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCg                                                                            <  1:290043/143‑1 (MQ=255)
                                                               gCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGt                                                                           <  1:669160/143‑1 (MQ=255)
                                                                 ggTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTaa                                                                         <  1:884827/143‑1 (MQ=255)
                                                                        gTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgtg                                                                  >  1:817514/1‑143 (MQ=255)
                                                                         tGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGc                                                                 >  1:270065/1‑143 (MQ=255)
                                                                                 gtcgcagcaTTTTTTCGGTGACGAGGTGCTCTGCATTTCACCGTGATTGATCTGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCa                                                         <  2:916445/135‑1 (MQ=255)
                                                                                 cccGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCa                                                         <  1:13903/143‑1 (MQ=255)
                                                                                    ggacactgctctCGCTGACGCGGTGCTCTGCATTTCACCGTGATTGATACTCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCCCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGtt                                                      <  1:382896/132‑1 (MQ=255)
                                                                                    ggAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGtt                                                      <  1:593787/143‑1 (MQ=255)
                                                                                         aTTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATGCGCCGGCTCTTCGTCGCATATTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTa                                                 <  1:465215/143‑1 (MQ=255)
                                                                                         aTTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTCCCACGGTTCGCAGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTa                                                 <  1:466253/143‑1 (MQ=255)
                                                                                          ttttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCCGGCTCTGCGTCGCCTCCTCTACTTTCCACGGTTCGCCGTACCTCACCACCATCCGCAGGACGATGTCGCACCCCTGTGCCGCTTACAGtagagcat                                                >  2:377833/1‑136 (MQ=255)
                                                                                          ttttctCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTAt                                                >  2:703189/1‑143 (MQ=255)
                                                                                               tCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATGTGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGCTACCGaac                                           <  2:851512/143‑1 (MQ=255)
                                                                                               tCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaac                                           <  1:651079/143‑1 (MQ=255)
                                                                                                             aTCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACAGTCTTTGaa                             >  1:736222/1‑143 (MQ=255)
                                                                                                               cTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGACGATGTCGTAATCCTGTGCCGCTTTCCGTTCGGTATCGCTCCACCCCTCTTTgctac                           >  1:255626/1‑139 (MQ=255)
                                                                                                               cTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAAc                           >  1:1098637/1‑143 (MQ=255)
                                                                                                                tGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGACGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTgcagcg                          >  2:191405/1‑138 (MQ=255)
                                                                                                                tGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTCTgcatcg                          >  1:518527/1‑138 (MQ=255)
                                                                                                                tGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCATGTGCAGCTTTCCGTTCGGTATCGAAAAAAACATATTTGAAACg                          >  2:619439/1‑143 (MQ=255)
                                                                                                                  tattccacctttcttgatATGTAGTCTCTGCGTCGCCTCTTCTACATTCCACGGCGCTCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGcc                        <  1:197493/130‑1 (MQ=255)
                                                                                                                        aCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCCGTTCGGTATCGAACAAACACTCTTTGAAACGCCCTGccc                  >  2:373457/1‑143 (MQ=255)
                                                                                                                        aCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTACCGAACAAACACTCTTTGAAGCGCCATGacc                  >  1:134803/1‑143 (MQ=255)
                                                                                                                         ccGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTGTCGAACAAACACTCTTTGCCACGCCATGCCCa                 >  1:46398/1‑143 (MQ=255)
                                                                                                                          cGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGACGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAg                >  2:1055770/1‑143 (MQ=255)
                                                                                                                          cGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCACTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAg                >  1:287778/1‑143 (MQ=255)
                                                                                                                           gTGATTGATATGCAGGCTATGCGTCGCCTCGTCTACTGTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAgg               <  1:300326/143‑1 (MQ=255)
                                                                                                                               ttgatCTGCAGCCTGTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGGACTCCACCACCACCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGCAACGCCATGCCCAGGTaaa           <  2:546371/143‑1 (MQ=255)
                                                                                                                               ttgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTaaa           <  1:913449/143‑1 (MQ=255)
                                                                                                                                  atatGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTCTCGACCAAACACTCCTTGAAACGCCATGCCCAGGTAAAGAc        >  1:877191/1‑143 (MQ=255)
                                                                                                                                  atatGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTGGGACCGCCATGCCCAGGTAAAGAc        >  2:578088/1‑143 (MQ=255)
                                                                                                                                   tatGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACCAACCCTCTTTGTCCCGCCCTGCCCCCGTAATGAcc       >  1:117751/1‑143 (MQ=255)
                                                                                                                                    atGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACCCCACCACCATCCGCCGGAAGATGTCGTCATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACCCTCTTTGAAACGCCATGCCCCGGTCCAGAccc      >  1:27369/1‑142 (MQ=255)
                                                                                                                                        aGGCTCTGCGTCGCCTCTTCTACTTTCCACTGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTAGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTCAAGACCTggct  >  1:844076/1‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                           
GATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTT  >  NC_000913/2132648‑2132926

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: