Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,309,066 T→C A701A (GCA→GCG)  pnp ← polynucleotide phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,309,0660TC95.1% 197.9 / ‑2.7 61A701A (GCA→GCG) pnppolynucleotide phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (27/31);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTACTTCCTGACCCATCTGCAGGTAATCGGTCACTTTCTCAAC  >  NC_000913/3308932‑3309194
                                                                                                                                      |                                                                                                                                
tCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCATATGGCAACCTTACTCGCCCTGTTCACCAGCCGGAGCTTCCGGcgctgc                                                                                                                             >  2:303423/1‑140 (MQ=255)
  ccGAAGACCGCGGTTGAATGAACGCCCTGTTCCCGGTTGCTAACAAGCCGTCCTGCCCGGTTAAAGGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTgcag                                                                                                                           <  2:179827/140‑1 (MQ=255)
  ccGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTgcag                                                                                                                           <  1:939793/140‑1 (MQ=255)
  ccGAAGACCACGGTTGAATGAACGTCCTGTTCCCGCTTGCTAACAAGGGGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTgcag                                                                                                                           <  2:436186/140‑1 (MQ=255)
   cGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATAGCAACCTTACTCACCCAGTTCAGCACCCGGAGCTTCCGGCGCTgcagc                                                                                                                          >  2:356846/1‑140 (MQ=255)
     aaGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGcagcag                                                                                                                        >  2:713478/1‑140 (MQ=255)
     aaGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGcagcag                                                                                                                        >  1:844569/1‑140 (MQ=255)
       gACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGt                                                                                                                      <  2:815731/140‑1 (MQ=255)
          cACGGTTGAATGACCGTCCTGTTCCCGGTTGCTAAAAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTga                                                                                                                   <  2:517955/140‑1 (MQ=255)
          cACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTga                                                                                                                   >  2:922153/1‑140 (MQ=255)
          cACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTga                                                                                                                   >  2:150261/1‑140 (MQ=255)
          cACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCGTACTCGCCGTGTACAGCAGCCGGAGGTTCCCGCGCTGCAGCAGGTTga                                                                                                                   >  2:324281/1‑140 (MQ=255)
          cACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCATACTCGCCATGTGCAGAAGCCGGAGTTTCCGGCGCTGCAGCAGGTTga                                                                                                                   >  2:376547/1‑140 (MQ=255)
            cGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTgaga                                                                                                                 >  1:108494/1‑140 (MQ=255)
            cGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTgaga                                                                                                                 <  2:54670/140‑1 (MQ=255)
            cGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTgaga                                                                                                                 >  1:505385/1‑140 (MQ=255)
            cGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCCCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTgaga                                                                                                                 >  1:949265/1‑140 (MQ=255)
               ttGAATGAACGTCCTTTTCCCGCTTCCTAGCAAGGCGGCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTg                                                                                                              <  1:325237/140‑1 (MQ=255)
                 gaatgaaCGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCt                                                                                                            <  2:505385/140‑1 (MQ=255)
                 gaatgaaCGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCt                                                                                                            <  1:150261/140‑1 (MQ=255)
                   atgaaCGTCCTGTTCCCGGTTGGTAACAACGCGTCCTCCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCa                                                                                                          <  1:336108/140‑1 (MQ=255)
                        cGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTTAGCAGCCGTAGCTTCCAGCGCTGCAGCAGGTTGAGACTGCTCAGTCGc                                                                                                     >  2:263457/1‑140 (MQ=255)
                                  ccGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATg                                                                                           >  1:166329/1‑140 (MQ=255)
                                       tGCTAACAAGGGGGCGTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAg                                                                                      <  2:844569/140‑1 (MQ=255)
                                        gCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGa                                                                                     <  1:108231/140‑1 (MQ=255)
                                               aGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATAc                                                                              >  2:994368/1‑140 (MQ=255)
                                                ggCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACg                                                                             >  2:987005/1‑140 (MQ=255)
                                                   gTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGcc                                                                          <  1:994368/140‑1 (MQ=255)
                                                    tCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGccc                                                                         <  1:922153/140‑1 (MQ=255)
                                                    tCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGa                                                                                     <  1:968349/128‑1 (MQ=255)
                                                    tCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGa                                                                                     >  2:968349/1‑128 (MQ=255)
                                                      cTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTg                                                                       >  2:1010604/1‑140 (MQ=255)
                                                      cTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTg                                                                       >  2:613240/1‑140 (MQ=255)
                                                          ccGGTTAAAAGCCCCTCGCCGCCGCAGAGCGCAAATGTCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGa                                                                   <  1:263457/140‑1 (MQ=255)
                                                           cGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGAt                                                                  <  1:242075/140‑1 (MQ=255)
                                                           cGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGAt                                                                  <  2:557221/140‑1 (MQ=255)
                                                             gTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCa                                                                <  2:108494/140‑1 (MQ=255)
                                                             gTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCa                                                                <  1:75835/140‑1 (MQ=255)
                                                              ttAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTgg                                                                      <  2:662339/133‑1 (MQ=255)
                                                              ttAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTgg                                                                      >  1:662339/1‑133 (MQ=255)
                                                                   aGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCGTTAGTGCTCAGACGGAGACGGAGCTGGCGATCAACTTcc                                                          >  2:444297/1‑140 (MQ=255)
                                                                       ccccGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGaa                                                      <  2:531104/140‑1 (MQ=255)
                                                                         ccgccgCAGCGGAGGGCAAATGGCATCCTCACTCGCCGTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACt                                                    <  1:303423/140‑1 (MQ=255)
                                                                         ccgccgCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACt                                                    <  2:102477/140‑1 (MQ=255)
                                                                             cgcAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCa                                                <  1:987005/140‑1 (MQ=255)
                                                                              gcAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCAc                                               <  1:713478/140‑1 (MQ=255)
                                                                               cAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCAcc                                              >  1:45783/1‑140 (MQ=255)
                                                                                         cAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTACTTCCt                                    <  1:181131/140‑1 (MQ=255)
                                                                                             tGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTACTTCCTGAcc                                <  2:45783/140‑1 (MQ=255)
                                                                                                cAACCTTACTCGCCCAGTTCAGAAGGCGGAGCTTCCGACGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTACTTCCTGACCCAt                             <  1:324281/140‑1 (MQ=255)
                                                                                                 aaCCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTACTTCCTGACCCATc                            >  1:11155/1‑140 (MQ=255)
                                                                                                             ccTGTTCAGCAGCCGGAGCTTCCGGTGCTCCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTACTTCCTGACCCATCTGCAGGTAATCg                <  1:356846/140‑1 (MQ=255)
                                                                                                             ccTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGGACTTTCAGCGGTACTTGCTGACCCATCTGCAGGTAATCg                >  2:328193/1‑140 (MQ=255)
                                                                                                             ccTGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTACTTCCTGACCCATCTGCAGGTAATCg                >  1:698295/1‑140 (MQ=255)
                                                                                                               tGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTACTTCCTGACCCATCTGCAGGTAATCGGt              <  2:949265/140‑1 (MQ=255)
                                                                                                               tGTTCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTACTTCCTGACCCATCTGCAGGTAATCGGt              <  2:886667/140‑1 (MQ=255)
                                                                                                                 ttCAGCAGCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTACCAGAACTATCACCGGCACTTCCTGGCCCAGCTGCAGGTAATCGGTCa            >  2:382116/1‑140 (MQ=255)
                                                                                                                   cagcagCCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTACTTCCTGACCCATCTGCAGGTAATCGGTCACt          <  1:613240/140‑1 (MQ=255)
                                                                                                                      cagcCGGAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGAGACGGCCCTGGCGATCAACTTCCAGAACTTTCACTGGTCCTTCCTGACCCACCTGCCGGTAATCGGTCACTTtc       >  2:428865/1‑140 (MQ=255)
                                                                                                                           ggAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTACTTCCTGACCCATCTGCAGGTAATCGGTCACTTTCTCAAc  >  2:711449/1‑140 (MQ=255)
                                                                                                                           ggAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTACTCCCTGACCCATCGGCAGGTAATCGGTCACTTTCTCAAc  >  1:322900/1‑140 (MQ=255)
                                                                                                                           ggAGCTTCCGGCGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTGCACCGGTAGTTCGTGAGCCATCTGCCGGTAATCGGTCACTTTCTCAAc  >  2:281817/1‑140 (MQ=255)
                                                                                                                                      |                                                                                                                                
TCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTACTTCCTGACCCATCTGCAGGTAATCGGTCACTTTCTCAAC  >  NC_000913/3308932‑3309194

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