Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,399,742 A→C D164A (GAT→GCT)  miaA → delta(2)‑isopentenylpyrophosphate tRNA‑adenosine transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,399,7420AC98.3% 219.7 / ‑4.6 59D164A (GAT→GCT) miaAdelta(2)‑isopentenylpyrophosphate tRNA‑adenosine transferase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (24/34);  total (24/35)
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.91e-01

GCTGGAAGGGTTGTCGCCGCTACCGTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGATCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGACGCAAACGTCAGGAGACGCTCTACCGTATCAGGTGCATCAGTTCGCCATCGCCCCGGCGAGCCGTGAACTGCTCCATCAACG  >  NC_000913/4399596‑4399889
                                                                                                                                                  |                                                                                                                                                   
gCTGGAAGGGTTGTCGCCGCTACCGTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGATCca                                                                                                                                                 <  1:10705/151‑1 (MQ=255)
gCTGGAAGGGTTGTCGCCGCTACCGTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGACATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTcc                                                                                                                                                  >  2:592129/1‑150 (MQ=255)
   ggAAGGGTTGTCGCCGCTACCGTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCcaca                                                                                                                                               <  2:509782/150‑1 (MQ=255)
    gAAGGGTTGTCGCCGCTACCGTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACaaa                                                                                                                                             <  1:71149/151‑1 (MQ=255)
     aaGGGTTGTCGCCGCTACCGTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACaa                                                                                                                                              >  1:282468/1‑149 (MQ=255)
      aGGGTTGTCGCCGCTACCGTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAgg                                                                                                                                           >  1:505280/1‑151 (MQ=255)
          ttGTCGCCGCTACCGTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCttt                                                                                                                                       >  1:900318/1‑151 (MQ=255)
             tCGCCGCTACCGTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTccc                                                                                                                                    >  2:83289/1‑151 (MQ=255)
                 cgcTACCGTCGGCAGAGCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGATGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGc                                                                                                                                >  2:199873/1‑151 (MQ=255)
                   cTACCGTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCAc                                                                                                                              >  1:511827/1‑151 (MQ=255)
                     aCCGTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTg                                                                                                                            <  2:351810/151‑1 (MQ=255)
                      ccGTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTgg                                                                                                                           >  1:924438/1‑151 (MQ=255)
                        gTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGaa                                                                                                                         <  2:77348/151‑1 (MQ=255)
                        gTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGaa                                                                                                                         <  2:526374/151‑1 (MQ=255)
                        gTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGaa                                                                                                                         <  2:284418/151‑1 (MQ=255)
                        gTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGaa                                                                                                                         <  2:226304/151‑1 (MQ=255)
                        gTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGa                                                                                                                          >  2:345063/1‑150 (MQ=255)
                          cGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAg                                                                                                                        <  2:878352/150‑1 (MQ=255)
                              aGACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGttt                                                                                                                     <  2:197155/149‑1 (MQ=255)
                               gACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGtttttt                                                                                                                  <  1:62398/151‑1 (MQ=255)
                               gACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGtttttt                                                                                                                  <  1:467381/151‑1 (MQ=255)
                                 cccGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTc                                                                                                                >  2:68516/1‑151 (MQ=255)
                                  ccGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCa                                                                                                               <  1:282569/151‑1 (MQ=255)
                                   cGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCAt                                                                                                              <  2:848578/151‑1 (MQ=255)
                                       aGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCg                                                                                                          <  1:587108/151‑1 (MQ=255)
                                           cGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTa                                                                                                      <  2:181904/151‑1 (MQ=255)
                                                cAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACt                                                                                                 <  2:356610/151‑1 (MQ=255)
                                                   aaTTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACttt                                                                                               >  2:469057/1‑150 (MQ=255)
                                                              aGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGAc                                                                                   >  1:267333/1‑151 (MQ=255)
                                                               ggcggcAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGACg                                                                                  >  1:628531/1‑151 (MQ=255)
                                                                     agagCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGACGCAAACg                                                                            <  2:260807/151‑1 (MQ=255)
                                                                         cAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGACGCAAACGTc                                                                          >  1:305134/1‑149 (MQ=255)
                                                                         cAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGACGCAAACGTCAg                                                                        <  1:922151/151‑1 (MQ=255)
                                                                               tGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGACGCAAACGTCAGGAGAc                                                                   >  2:370952/1‑150 (MQ=255)
                                                                                  gAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGACGCAAACGTCAGGAGACGctc                                                               <  1:339470/151‑1 (MQ=255)
                                                                                    gTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGACGCAAACGTCAGGAGACGctct                                                              >  2:794857/1‑150 (MQ=255)
                                                                                       aTTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGACGCAAACGTCAGGAGACGCTCTACCg                                                          <  1:847319/151‑1 (MQ=255)
                                                                                         tGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGACGCAAACGTCAGGGGACGCTCTACCGt                                                         <  2:328898/150‑1 (MQ=255)
                                                                                         tGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGACGCAAACGTCAGGAGACGCTCTACCGt                                                         >  1:837029/1‑150 (MQ=255)
                                                                                            aTCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGACGCAAACGTCAGGAGACGCTCTACCGTATc                                                      >  1:18056/1‑150 (MQ=255)
                                                                                            aTCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGACGCAAACGTCAGGAGACGCTCTACCGTATCa                                                     <  1:592144/151‑1 (MQ=255)
                                                                                                 cAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGACGCAAACGTCAGGAGACGCTCTACCGTATCAGGTGc                                                <  2:481566/151‑1 (MQ=255)
                                                                                                   aCTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGACGCAAACGTCAGGAGACGCTCTACCGTATCAGGTGCAt                                              <  2:756980/151‑1 (MQ=255)
                                                                                                       cAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGACGCAAACGTCAGGAGACGCTCTACCGTATCAGGTGCATCAg                                           <  2:636065/150‑1 (MQ=255)
                                                                                                        aggaggTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGACGCAAACGTCAGGAGACGCTCTACCGTATCAGGTGCATCAGtt                                         >  1:886069/1‑151 (MQ=255)
                                                                                                            ggTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGACGCAAACGTCAGGAGACGCTCTACCGTATCAGGTGCATCAGTTCGcc                                     <  2:665175/151‑1 (MQ=255)
                                                                                                             gTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGACGCAAACGTCAGGAGACGCTCTACCGTATCAGGTGCATCAGTTCGCCa                                    <  1:262388/151‑1 (MQ=255)
                                                                                                               aGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGACGCAAACGTCAGGAGACGCTCTACCGTATCAGGTGCATCAGTTCGCCATc                                  >  1:801396/1‑151 (MQ=255)
                                                                                                                gATCCGGTTGCGGGAGCAAGGATACATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGACGCAAACGTCAGGAGACGCTCTTCCGTATCAGGTGCAACAGTACGCCATCg                                 <  2:715935/151‑1 (MQ=255)
                                                                                                                   ccGGTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGACGCAAACGTCAGGAGACGCTCTACCGTATCAGGTGCATCAGTTCGCCATCGccc                              >  1:119537/1‑151 (MQ=255)
                                                                                                                     ggTTGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGACGCAAACGTCAGGAGACGCTCTACCGTATCAGGTGCATCAGTTCGCCATCGCCCCg                            >  1:742574/1‑151 (MQ=255)
                                                                                                                        tGCGGCAGCAAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGACGCAAACGTCAGGAGACGCTCTACCGTATCAGGTGCATCAGTTCGCCATCGCCCCGGc                          <  1:742482/150‑1 (MQ=255)
                                                                                                                               gcaAGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGACGCAAACGTCAGGAGACGCTCTACCGTATCAGGTGCATCAGTTCGCCATCGCCCCGGCGAGCCGTg                  <  1:83289/151‑1 (MQ=255)
                                                                                                                                 aaGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGACGCAAACGTCAGGAGACGCTCTACCGTATCAGGTGCATCAGTTCGCCATCGCCCCGGCGAGCCGTGa                 >  2:37411/1‑150 (MQ=255)
                                                                                                                                 aaGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGACGCAAACGTCAGGAGACGCTCTACCGTATCAGGTGCATCAGTTCGCCATCGCCCCGGCGAGCCGTGa                 >  1:737953/1‑150 (MQ=255)
                                                                                                                                  aGGATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGACGCAAACGTCAGGAGACGCTCTACCGTATCAGGTGCATCAGTTCGCCATCGCCCCGGCGAGCCGTGAAc               <  1:446244/151‑1 (MQ=255)
                                                                                                                                   ggATTCATCCAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGACGCAAACGTCAGGAGACGCTCTACCGTATCAGGTGCATCAGTTCGCCATCGCCCCGGCGAGCCGTGAACt              <  1:68516/151‑1 (MQ=255)
                                                                                                                                            cAAATGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGACGCAAACGTCAGGAGACGCTCTACCGTATCAGGTGCATCAGTTCGCCATCGCCCCGGCGAGCCGTGAACTGCTCCATCa     <  1:551225/151‑1 (MQ=255)
                                                                                                                                               aTGCTCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGACGCAAACGTCAGGAGACGCTCTACCGTATCAGGTGCATCAGTTCGCCATCGCCCCGGCGAGCCGTGAACTGCTCCATCAACg  <  1:140847/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCTGGAAGGGTTGTCGCCGCTACCGTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCATCGTCAACTTCAGGAGGTAGATCCGGTTGCGGCAGCAAGGATTCATCCAAATGATCCACAAAGGCTTTCCCGGGCACTGGAAGTTTTTTTCATTTCGGGTAAAACTTTAACGGAACTGACGCAAACGTCAGGAGACGCTCTACCGTATCAGGTGCATCAGTTCGCCATCGCCCCGGCGAGCCGTGAACTGCTCCATCAACG  >  NC_000913/4399596‑4399889

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: