Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,815,599 T→G K281Q (AAA→CAA)  gshA ← glutamate‑cysteine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,815,5990TG98.3% 214.4 / ‑4.2 59K281Q (AAA→CAA) gshAglutamate‑cysteine ligase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (17/41);  total (18/41)
Fisher's exact test for biased strand distribution p-value = 3.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

CACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGA  >  NC_000913/2815453‑2815741
                                                                                                                                                  |                                                                                                                                              
cACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTgctc                                                                                                                                             >  2:566842/1‑150 (MQ=255)
    tCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATAc                                                                                                                                        >  1:917392/1‑151 (MQ=255)
     cAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATAcc                                                                                                                                       <  2:643896/151‑1 (MQ=255)
      aaTATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCa                                                                                                                                      <  2:4619/151‑1 (MQ=255)
       atatatTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCa                                                                                                                                      >  1:26299/1‑150 (MQ=255)
         atatTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAAt                                                                                                                                    <  1:701778/150‑1 (MQ=255)
             tCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTc                                                                                                                                <  2:678963/150‑1 (MQ=255)
              cAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTcgc                                                                                                                              >  2:374253/1‑151 (MQ=255)
                  gccgccACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTAc                                                                                                                          >  1:121621/1‑151 (MQ=255)
                  gccgccACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTAc                                                                                                                          <  2:751835/151‑1 (MQ=255)
                      ccACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCtt                                                                                                                      >  2:506322/1‑151 (MQ=255)
                      ccACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCtt                                                                                                                      >  2:250842/1‑151 (MQ=255)
                        aCGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTcc                                                                                                                    <  1:537229/151‑1 (MQ=255)
                        aCGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTcc                                                                                                                    <  1:846061/151‑1 (MQ=255)
                            aaCAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGAt                                                                                                                 <  1:652928/150‑1 (MQ=255)
                              cAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGc                                                                                                              <  1:874252/151‑1 (MQ=255)
                                   cATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGtttt                                                                                                         >  1:766916/1‑151 (MQ=255)
                                   cATCAGAAGGCGACTCGCCACTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGtttt                                                                                                         <  2:248500/151‑1 (MQ=255)
                                            gCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTg                                                                                                <  2:869299/151‑1 (MQ=255)
                                             cGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGt                                                                                               >  2:562135/1‑151 (MQ=255)
                                                     cgcTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGcc                                                                                       <  2:514751/151‑1 (MQ=255)
                                                          cGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTa                                                                                  <  2:163468/151‑1 (MQ=255)
                                                            ggTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACg                                                                                <  1:350842/151‑1 (MQ=255)
                                                            ggTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACg                                                                                <  2:531992/151‑1 (MQ=255)
                                                            ggTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACg                                                                                <  2:758941/151‑1 (MQ=255)
                                                                 cgcgTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTc                                                                           <  2:886877/151‑1 (MQ=255)
                                                                 cgcgTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTc                                                                           <  1:710746/151‑1 (MQ=255)
                                                                    gTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGGCCGCTACGTACTCGt                                                                         >  1:384197/1‑150 (MQ=255)
                                                                     gtttGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTa                                                                        >  2:640565/2‑150 (MQ=255)
                                                                       ttGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTaaa                                                                      <  1:250847/150‑1 (MQ=255)
                                                                        tGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTaa                                                                       >  2:495111/1‑148 (MQ=255)
                                                                         ggACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAg                                                                     <  2:213395/149‑1 (MQ=255)
                                                                         ggACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGa                                                                    <  2:32637/150‑1 (MQ=255)
                                                                         ggACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGAt                                                                   <  1:343925/151‑1 (MQ=255)
                                                                         ggACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCGCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGAt                                                                   <  2:457436/151‑1 (MQ=255)
                                                                          gACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTaa                                                                       >  2:882173/1‑146 (MQ=255)
                                                                             gAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGtt                                                               <  2:435067/151‑1 (MQ=255)
                                                                                      cgTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAAt                                                      <  1:584980/151‑1 (MQ=255)
                                                                                       gTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATa                                                     <  1:135953/151‑1 (MQ=255)
                                                                                               cGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGAtt                                             <  2:474423/151‑1 (MQ=255)
                                                                                               cGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGAt                                              <  1:562150/150‑1 (MQ=255)
                                                                                                  tttCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGc                                           <  1:795077/150‑1 (MQ=255)
                                                                                                    tCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCttt                                        >  1:882326/1‑151 (MQ=255)
                                                                                                     cAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCttt                                        <  1:480151/150‑1 (MQ=255)
                                                                                                     cAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCttt                                        >  2:824005/1‑150 (MQ=255)
                                                                                                     cAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTg                                       <  1:757741/151‑1 (MQ=255)
                                                                                                     cAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCGTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCttt                                        <  2:86587/150‑1 (MQ=255)
                                                                                                           gCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGAtt                                  <  1:242211/150‑1 (MQ=255)
                                                                                                                 cGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGt                           <  1:312519/151‑1 (MQ=255)
                                                                                                                    tGCTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGtat                         <  2:240484/150‑1 (MQ=255)
                                                                                                                      cTGTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGc                      >  1:846636/1‑151 (MQ=255)
                                                                                                                        gTTGATTTGCAGCCTCTTACCGTCTTGCTCAATCCCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGca                     <  2:843891/150‑2 (MQ=255)
                                                                                                                        gTTGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGc                      <  2:215126/149‑1 (MQ=255)
                                                                                                                          tGATTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCg                    >  1:745322/1‑149 (MQ=255)
                                                                                                                            aTTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCg                    >  1:483648/1‑147 (MQ=255)
                                                                                                                            aTTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGAt                <  1:105538/151‑1 (MQ=255)
                                                                                                                            aTTTGCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGAt                <  1:847309/151‑1 (MQ=255)
                                                                                                                                gCAGCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGc             <  2:889743/150‑1 (MQ=255)
                                                                                                                                   gCCTCTTACCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCa          >  1:285011/1‑150 (MQ=255)
                                                                                                                                          aCCGTCTTGCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAaga  <  1:224740/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
CACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGA  >  NC_000913/2815453‑2815741

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑