Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,225,493 C→T 7.5% G47S (GGC→AGC)  yohC ← putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,225,4930CT 7.5% 155.0 / 3.6 62G47S (GGC→AGC) yohCputative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (24/33);  new base T (2/3);  total (26/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.76e-01

CATCCACCAGATGACCCGCCCCATCACCGCGACCCCAGCCAGCATAACGCCGTAAAACAAGACAGCCAGCGCCAGTCCGGTAAACCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCAATCTGTGTAGTGCCAATGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACAGTCCCCAGACATGGCTCATA  >  NC_000913/2225353‑2225632
                                                                                                                                             |                                                                                                                                           
cATCCAGCAAATCCCCCGCCCCACCACCGCGCCACCAGACACAATAACGCCGTCCAATAAGACGGCCAGCGCCAGTCCGGTAAACAACGATAACTTCAAGATAGTTCAATCGCCAAAATTCCACCCAATCTGTGTAGTGCCaa                                                                                                                                           <  2:582735/143‑1 (MQ=255)
cATCCACCAGATGACCCGCCCCATCACCGCGACCCCAGCCAGCATAACGCCGTAAAACAAGACAGCCAGCGCCAGTCCGGTAAACCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCAATCTGTGTAGTGCCaa                                                                                                                                           <  1:231533/143‑1 (MQ=255)
    caccaGATGACCCTCCCCGCGCCCGCGACCCCAGCCAGCATAACGCCGTAACACAAGGCAGCCAGCGCCAGTCCGGTAAACCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCAATCTGTGTAGTGCCAATGaa                                                                                                                                       <  2:529632/143‑1 (MQ=255)
    caccaGATGACCCGCCCCATCACCGCGACCCCAGCCAGCATAACGCCGTAAAACAAGACAGCCAGCGCCAGTCCGGTAAACCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCAATCTGTGTAGTGCCAATGaa                                                                                                                                       >  1:854147/1‑143 (MQ=255)
     accaGATGACCCGCCCCATCACCGCGACCCCAGCCAGCATAACGCCGTAAAACAAGACAGCCAGCGCCAGTCCGGTAAACCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCAATCTGTGTAGTGCCAATGAAg                                                                                                                                      >  1:178734/1‑143 (MQ=255)
       caGATGCCCCCCCCCATCACCCCGACCCCCGCCAGCATAACGCCGTAAAACAAGACAGCCAGCGCCAGTCCGGTAAACCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCAATCTGTGTAGTGCCAATGAAGgc                                                                                                                                    <  1:178787/143‑1 (MQ=255)
         gaggACCCGCACCATCGCGGCGACCGCAGCCAGCATAACGCCGTAAAACAAGGCAGCCAGCGCCAGTCCGGTAAACCAGGATAACTTCAGGATAGTGCGATCGCCAAAATTCCAGCCAATCTGTGTAGTGCCAATGAAGgcgc                                                                                                                                  <  2:627915/140‑1 (MQ=255)
           tGACCCACTCCATCCCCGCGACCCCAGCCAGCATAACGCCGTAAAACAAGACAGCCAGCGCCAGTCCGGTAAACCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCAATCTGTGTAGTGCCAATGAAGGCGCaa                                                                                                                                <  1:581212/143‑1 (MQ=255)
            gccccGCCCCACCACCGCGGCCCCAGCCAGCATAACGCCGTAGGCGGAGACAGCCAGCGCCAGTCCGGTAAACCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCAATCTGTGTAGTGCCAATGAAGGCGCaaa                                                                                                                               <  1:369107/141‑1 (MQ=255)
             aCCCGCCCCGCCACCGCGCCCCCAGCCAGCATAACGCCGTAAAACAAGACAGCCCGCGCCAGTCCGGTAAACCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCAATCTGTGTAGTGCCAATGAAGGCGCAAAt                                                                                                                              <  1:284173/143‑1 (MQ=255)
                cGCCCCATCACCGCGACCCCAGCCAGCATAACGCCGTAAAACAAGACAGCCAGCGCCAGTCCGGTAAACCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCAATCTGTGTAGTGCCAATGCAGGCGCAAATCAc                                                                                                                           >  1:297339/1‑143 (MQ=255)
                cGCCCCATCACCGCGACCCCAGCCAGCATAACGCCGTAAAACAAGACAGCCAGCGCCAGTCCGGTAAACCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCAATCTGTGTAGTGCCAATGAAGGCGCAAATCAc                                                                                                                           >  1:766904/1‑143 (MQ=255)
                      aTCACCGCGACCCCAGCCAGCATAACGCCGTAAAACAGGACAGCCAGCGCCAGTCCGGTAAACCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCAATCTGTGTAGTGCCAATGAAGGCGCAAATCACCGGAAt                                                                                                                     <  1:335130/143‑1 (MQ=255)
                       tCACCGCGACCCCAGCCAGCATAACGCCGTAAAACAAGACAGCCAGCGCCAGTCCGGTAAACCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCAATCTGTGTAGTGCCAATGAAGGCGCAAATCACCGGAATc                                                                                                                    >  1:664304/1‑143 (MQ=255)
                        cACCGCGACCCCAGCCAGCATAACGCCGTAAAACAAGACAGCCAGCGCCAGTCCGGTAAACCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCAATCTGTGTAGTGCCAATGAAGGCGCAAATCACCGGAATcg                                                                                                                   <  1:309304/143‑1 (MQ=255)
                                ccccAGCCAGCATAACGCCGTAAAACAAGACAGCCAGCGCCAGTCCGGTAAACCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCAATCTGTGTAGTGCTAATGAAGGCGCAAATCACCGGAATCGCCGCCATc                                                                                                           >  1:538210/1‑143 (MQ=255)
                                ccccAGCCAGCATAACGCCGTAAAACAAGACAGCCAGCGCCAGTCCGGTAAACCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCAATCTGTGTAGTGCTAATGAAGGCGCAAATCACCGGAATCGCCGCCATc                                                                                                           >  1:110295/1‑143 (MQ=255)
                                  ccagccagCATAACGCCGTAAAACAAGACAGCCAGCGCCAGTCCGGTAACCCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCCATCTGTGTAGTGCCAATGAAGGCGCAAATCACCGGAATCGCCGCCATCAg                                                                                                         >  2:245293/1‑143 (MQ=255)
                                    agccagcATAACGCCGTAAAACAAGACAGCCAGCGCCAGTCCGGTAAACCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCAATCTGTGTAGTGCCAATGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCa                                                                                                       >  1:52055/1‑143 (MQ=255)
                                      ccagcATAACGCCGTAAAACAAGACAGCCAGCGCCAGTCCGGTAAACCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCAATCTGTGTAGTGCCAATGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCaaa                                                                                                     >  2:210634/1‑143 (MQ=255)
                                      ccagcATAACGCCGTAAAACAAGACAGCCAGCGCCAGTCCGGTAAACCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCAATCTGTGTAGTGCCAATGAAGGCGCAAATCACCGGAATCGCCGCCACCAGCaac                                                                                                     >  1:171683/1‑142 (MQ=255)
                                                 ccGTAAAACAAGACAGCCAGCGCCAGTCCGGTAAACCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCAGTCTGTGTAGTGCCAATGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCAAAACGTggtgggt                                                                                          >  1:309019/1‑143 (MQ=255)
                                                   gTAAAACAAGACAGCCAGCGCCAGTCCGGTAAACCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCAATCTGTGTAGTGCCAATGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCAAAACGTGGTGGgtgt                                                                                        >  2:474641/1‑143 (MQ=255)
                                                     aaaaCGAGAGAGCCAGCGCCAGTCCGGTAAACCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCAATCTGTGTAGTGCCAATGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCAAAACGTGGTGGGTGTaa                                                                                      <  2:32254/143‑1 (MQ=255)
                                                      aaaCAAGACAGCCAGCGCCAGTCCGGTAAACCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCAATCTGTGTAGTGCCAATGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCAAAACGTGGTGGGTGTAAt                                                                                     >  1:238329/1‑143 (MQ=255)
                                                       aaCAAGACAGCCAGCGCCAGTCCGGTAAACCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCAATCTGTGTAGTGCTAATGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCAAAACGTGGTGGGTGTAATg                                                                                    <  2:538210/143‑1 (MQ=255)
                                                       aaCAAGACAGCCAGCGCCAGTCCGGGAAACCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCAATCTGTGTAGTGCTAATGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCAAAACGTGGTGGGTGTAATg                                                                                    <  2:110295/143‑1 (MQ=255)
                                                          aaGACAGCCAGCGCCAGACCGGTAAACCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCAATCTGTGTAGTGCCAATGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCAAAACGTGGTGGGTGTAATGGTg                                                                                 <  1:99348/143‑1 (MQ=255)
                                                             aCAGCCAGCGCCAGTCCGGTAAACCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCAATCTGTGTAGTGCCAATGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCAAAACGTGGTGGGTGTAATGGTgaga                                                                              >  1:342405/1‑143 (MQ=255)
                                                             aCAGCCAGCGCCAGTCCGGTAAACCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCAAACTGTGTAGTGCCAATGAAGGAGCAAATCACCGGCATGGCCACCATCAGCAAAACGTGGTAGGGGGAGTGTTgaga                                                                              >  2:529592/1‑143 (MQ=255)
                                                               agccagcGCCAGTCCGGTAAACCAGGATAACTTGGGGCAACCACCAGCGCCAAAATTCCAGCCAATCTGTGTAGTGCCAATGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCAAAACGTGGTGGGTGTAATGGTGAGaaa                                                                            <  2:854147/143‑1 (MQ=255)
                                                                 ccagcGCCAGTCCGGTAAACCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCAATCTGTGTAGTGCCAATGAAGGCGCAAATCACCGGCATCGCCGCCATCAGCAAAACGTGGTGGGTGTAATGGTGCGAAAtc                                                                          >  1:615957/1‑143 (MQ=255)
                                                                 ccagcGCCAGTCCGGTAAACCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCAATCTGTGTAGTGCCAATGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCAAAACGTGGTGGGTGTAATGGTGAGgaaac                                                                          >  2:517780/1‑141 (MQ=255)
                                                                   agcGCCAGTCCGGTAAACCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCAATCTGGGTAGTGCCAATGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCAAAACGTGGTGGGTGGAGTGGTGAGAGAtcag                                                                        >  2:82164/1‑141 (MQ=255)
                                                                    gcgcCTGTCCGGTAAACCAGGATAACGTCAGGATAGTGCCATCGCCAAAATACCTGCTCTATTGTGTAGTGCCAATGAAGGCGCAAATCACCGGTATCGACGCCATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAAtcgtt                                                                       <  2:101765/143‑1 (MQ=255)
                                                                     cgcCAGTCCGGTAAACCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCAATCTGTGTAGTGCCAATGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAAtcgttt                                                                      >  2:791596/1‑143 (MQ=255)
                                                                             cGGTAAACCAGGATATCTTCAGGATAGTGCCACCGCCCAAATTCCAGCCAATCTGTGTAGTGCCACTGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTcg                                                              <  1:92405/143‑1 (MQ=255)
                                                                             cGGTAAACCAGGAGAACTTCAGGATAGGGCCATCGCCAAAATGCCAGCCAATCTGTGTAGTGCCAATGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTcg                                                              <  1:261082/143‑1 (MQ=255)
                                                                               ggaaaCCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCAATCTGTGTAGTGCCAATGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTcgcg                                                            <  1:97812/141‑1 (MQ=255)
                                                                                        gatacattCAGGATGGCGCCATCGCAAAAAGTCCAGCCAATCTGTGTAGTGCCAATGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCAc                                                   <  2:178734/137‑1 (MQ=255)
                                                                                           aaCTTCAGGATAGGGCCAGCGCCAAAATTCCAGCCAATCTGTGTAGTGCCAATGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTg                                                <  2:586856/143‑1 (MQ=255)
                                                                                                cAGGATAGTGCCATCGCCAAAATTCCAGCCAATCTGTGTAGTGCCAATGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCAttt                                           <  1:426143/143‑1 (MQ=255)
                                                                                                    aTAGTGCCATCGCCAAAATTCCAGCCGATCTGTGTAGTGCCAATGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACg                                       <  2:746113/143‑1 (MQ=255)
                                                                                                         gCCATCGCCAAAATTCCAGCCAATCTGTGTAGTGCCAATGAAGGCGCAAATCACCGGAATCGCCGCCAGCAGCAAAACGTGGTGGGAGTAATGGTGAGAACTCGTTGCGTTTGCGCGATTGATAACCTGCATTCCACGATCgg                                  >  2:350035/1‑143 (MQ=255)
                                                                                                           cATCGCCAAAAGTCCAGCCAATCTGTGTAGTGCCAATGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGa                                <  1:285846/143‑1 (MQ=255)
                                                                                                              cGCCAAAATTCCAGCCAATCTGTGTAGTGCCAATGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATgg                             <  2:342405/143‑1 (MQ=255)
                                                                                                                ccAAAATTCCAGCCAATCTGTGTAGTGCCAATGAAGGCGCAAATCACCGGAATCGCCTCCATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGAATGGGa                           >  2:244052/1‑143 (MQ=255)
                                                                                                                  caaaTTCCAGCCAACCCGCCTATGTCTAATGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCAAATCGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGaaa                         <  2:138698/142‑1 (MQ=255)
                                                                                                                  aaaaTTCCAGCCAATCTGTGTAGTGCCAATGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGAACGGGATGGGaaa                         >  2:505517/1‑143 (MQ=255)
                                                                                                                   aaaTTCCAGCCAATCTGTGTAGTGCCAATGAAGGCGCAAATCACCGGAATCGCCGCCATCTGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGaaaa                        <  1:292569/143‑1 (MQ=255)
                                                                                                                    tgttcgtgCCAAGCTGCGTAGTTCCACTGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCAAAGCGTGGTGTGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAAc                       <  2:601916/136‑1 (MQ=255)
                                                                                                                    aaTTCCTGCCAAGCTGTGTAGTGCCAATGAAGGCGTCAATCTGCGGAATCGCCGCCATCTGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAAc                       <  2:432148/143‑1 (MQ=255)
                                                                                                                     aTTCCAGCCAATCTGTGTAGTGCCAATGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACa                      <  1:791596/143‑1 (MQ=255)
                                                                                                                        ccAGCCAATCTGTGTAGTGCCAATGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACAGTc                   >  1:659392/1‑143 (MQ=255)
                                                                                                                           aggcgaTCTGTGTAGTGCCAATGTAGGCGCAAATCACCGGAATCGCTGTCCTCAGCAAATCGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCTCGATTAATCCCCTGCATTTCACGATCGGGATGGGAAAACAGTccc                 <  1:614937/138‑1 (MQ=255)
                                                                                                                                  gtgtgTAGTGCCAATGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACAGTCCCCAGACAt          <  2:688363/142‑1 (MQ=255)
                                                                                                                                  cTGTGTAGTGCCAATGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACAGTCCCCAGACAt          <  2:221947/143‑1 (MQ=255)
                                                                                                                                    gtgtAGTGCCAATGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAACCAGTCCCCAGACATgg        >  2:357029/1‑143 (MQ=255)
                                                                                                                                    gtgccatGCCCATGACGCCGCAAACCACCGGGATCCTCCCCATCAGCACAACGTGCTGGGTGTAATGGTGCGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACAGTCCCCAGACATgg        <  2:685172/137‑1 (MQ=255)
                                                                                                                                    gcgtggTGCCAATGAAGGGGGCAAGCACCGGAATCGCCGCCATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTGCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACAGTCCCCAGACATgg        <  2:750231/138‑1 (MQ=255)
                                                                                                                                        aGTGCCAATGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACAGTCCCCAGACATGGCTCa    <  1:13627/143‑1 (MQ=255)
                                                                                                                                          tGCCAATGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACAGGCCCCAGACATGGGTcgtg  >  1:419860/1‑140 (MQ=255)
                                                                                                                                             |                                                                                                                                           
CATCCACCAGATGACCCGCCCCATCACCGCGACCCCAGCCAGCATAACGCCGTAAAACAAGACAGCCAGCGCCAGTCCGGTAAACCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCAATCTGTGTAGTGCCAATGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACAGTCCCCAGACATGGCTCATA  >  NC_000913/2225353‑2225632

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: