Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 271,335 A→C 8.5% R412R (CGA→CGC)  PP_0218 → sensory box protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947271,3350AC 8.5% 179.4 / 5.4 59R412R (CGA→CGC) PP_0218sensory box protein
Reads supporting (aligned to +/- strand):  ref base A (28/26);  new base C (3/2);  total (31/28)
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.88e-01

ACGGGCGGGCCTGCTATGCGCTGTACACCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCGGGTGAGTTCATCCCGATTGCCGAACGCACCGGCCTGATTGCCGAGATCGATACCTGGGTGCTGCGCCAGGCCTGCCGGCAAATGGT  >  NC_002947/271187‑271481
                                                                                                                                                    |                                                                                                                                                  
aCGGGCGGGCCTGCTATGCGCTGTACACCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGTGTCAGGCCGa                                                                                                                                                    >  8:26703/1‑149 (MQ=255)
aCGGGCGGGCCTGCTATGCGCTGTACACCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGa                                                                                                                                                    <  1:368/149‑1 (MQ=255)
    gcgggcCTGCTATGCGCTGTACACCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAgg                                                                                                                                                >  6:165340/1‑149 (MQ=255)
    gcgggcCTGCTATGCGCTGTACACCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAgg                                                                                                                                                <  5:24637/149‑1 (MQ=255)
         ccTGCTATGCGCTGTACACCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCTGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCggg                                                                                                                                           <  1:194636/149‑1 (MQ=255)
             cTATGCGCTGTACACCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTc                                                                                                                                        >  3:260926/1‑148 (MQ=255)
             cTATGCGCTGTACACCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTc                                                                                                                                        <  2:268852/148‑1 (MQ=255)
                  cgcTGTACACCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACg                                                                                                                                  >  3:366745/1‑149 (MQ=255)
                       tACACCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGTGTCAGGCCGACAGGTCGGGGTCGAGACGCTg                                                                                                                               <  7:26703/147‑1 (MQ=255)
                       tACACCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTg                                                                                                                               <  1:360613/147‑1 (MQ=255)
                       tACACCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTg                                                                                                                               <  5:240133/147‑1 (MQ=255)
                         cacCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTgc                                                                                                                           >  2:231241/1‑149 (MQ=255)
                          acCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTgcg                                                                                                                          >  2:238901/1‑149 (MQ=255)
                               agagCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGgc                                                                                                                     <  4:269442/149‑1 (MQ=255)
                                        cgcccacgcccaGCATCGGGTCGAGAACGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCAccc                                                                                                              >  2:229060/1‑147 (MQ=255)
                                          cccacgcccaGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAgg                                                                                                          >  8:247838/1‑149 (MQ=255)
                                           ccacgcccaGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCNGGTGCGCTGGCAGCACCCCNGgc                                                                                                         >  6:160924/1‑149 (MQ=255)
                                               gcccaGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCGGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCgg                                                                                                     <  2:188629/149‑1 (MQ=255)
                                                   aGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTg                                                                                                 <  5:165340/149‑1 (MQ=255)
                                                              gagaCCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCggg                                                                                      >  6:107003/1‑149 (MQ=255)
                                                              gagaCCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCggg                                                                                      >  3:39063/1‑149 (MQ=255)
                                                               agaCCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCGGGt                                                                                     <  1:231241/149‑1 (MQ=255)
                                                               agaCCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCGGGt                                                                                     <  7:247838/149‑1 (MQ=255)
                                                                 aCCGCCGCTGAATTGCGTCGGGCCCTGGAGGAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCGGGt                                                                                     >  8:141012/1‑147 (MQ=255)
                                                                        cTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTgg                                                                                                >  7:308703/1‑129 (MQ=255)
                                                                        cTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTgg                                                                                                <  8:308703/129‑1 (MQ=255)
                                                                        cTGAATTGCGGCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGGCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCGGGTGAGTTCATc                                                                            <  6:254339/149‑1 (MQ=255)
                                                                           aaTTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCGGGTGAGTTCATccc                                                                          >  6:344576/1‑148 (MQ=255)
                                                                             ttGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCGGGTGAGTTCATCCCGAt                                                                       <  1:229060/149‑1 (MQ=255)
                                                                                cGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCGGGTGAGTTCATCCCGATTg                                                                     <  4:366745/148‑1 (MQ=255)
                                                                                 gTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTAcc                                                                                            >  1:250148/1‑124 (MQ=255)
                                                                                 gTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTAcc                                                                                            <  2:250148/124‑1 (MQ=255)
                                                                                  tCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCGGGTGAGTTCATCCCGATTGcc                                                                   >  1:44752/1‑148 (MQ=255)
                                                                                          tGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGgcg                                                                                                        <  2:356405/103‑1 (MQ=255)
                                                                                          tGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGgcg                                                                                                        >  1:356405/1‑103 (MQ=255)
                                                                                             aGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGGCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCGGGTGAGTTCATCCCGATTGCCGAACGCACCGGc                                                       <  1:238901/149‑1 (MQ=255)
                                                                                              gCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCGGGTGAGTTCATCCCGATTGCCGAACGCACCGGcc                                                      <  4:39063/149‑1 (MQ=255)
                                                                                              gCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCGGGTGAGTTCATCCCGATTGCCGAACGCACCGGcc                                                      >  3:344422/1‑149 (MQ=255)
                                                                                                    tGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCGGGTGAGTTCATCCCGAtt                                                                      >  6:53417/1‑127 (MQ=255)
                                                                                                    tGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCGGGTGAGTTCATCCCGAtt                                                                      <  5:53417/127‑1 (MQ=255)
                                                                                                     gAACTGCTGGTGTACTTCCAGCAGGTACACGACATCGACTCAGGCCGACAGCTCCCGCTCGACAAGCTGGTTCGATGGCATCACCCCAGGCGCCTTCTGGTACCGCCGGGTGAGATCAACACGATTGACGAACGCAACGGCCTGATTg                                                >  4:10957/1‑148 (MQ=255)
                                                                                                     gAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCGGGTGAGTTCATCCCGATTGCCGAACGCACCGGCCTGATTg                                                <  4:344422/148‑1 (MQ=255)
                                                                                                      aaCTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCGGGTGAGTTCATCCCGATTGCCGAACGCACCGGCCTGATTGcc                                              <  1:56990/149‑1 (MQ=255)
                                                                                                      aaCTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCGGGTGAGTTCATCCCGATTGCCGAACGCACCGGCCTGATTGcc                                              <  7:333810/149‑1 (MQ=255)
                                                                                                       aCTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCGGGTGAGTTCATCCCGATTGCCGAACGCACCGGCCTGATTGcc                                              <  3:10957/148‑1 (MQ=255)
                                                                                                       aCTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCGGGTGAGTTCATCCCGATTGCCGAACGCACCGGCCTGATTGcc                                              >  8:202250/1‑148 (MQ=255)
                                                                                                                  aCTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCGGGTGAGTTCATCCCGATTGCCGAACGCACCGGCCTGATTGCCGAGAt                                         <  5:82173/142‑1 (MQ=255)
                                                                                                                  aCTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCGGGTGAGTTCATCCCGATTGCCGAACGCACCGGCCTGATTGCCGAGAt                                         >  6:82173/1‑142 (MQ=255)
                                                                                                                    ttCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCGGGTGAGTTAATCCCGATTGCCGAACGCACCGGCCTGATTGCCGAGATCGATACCTg                                >  3:225727/1‑149 (MQ=255)
                                                                                                                       cAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCGGGTGAGTTCATCCCGATTGCCGAACGCACCGGCCTGATTGCCGAGATCGATACCTGGGt                             >  3:82051/1‑149 (MQ=255)
                                                                                                                                cccGACCTGTCGGCCTGACGCCAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCGGGTGAg                                                                                  <  5:225747/85‑1 (MQ=37)
                                                                                                                                cccGACCTGTCGGCCTGACGCCAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCGGGTGAg                                                                                  >  6:225747/3‑87 (MQ=37)
                                                                                                                                cccGACCTGTCGGCCTGACGCCAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCGGGTGAGTTCATCCCGATTGCCGAACGCACCGGCCTGATTGCCGAGATCGATACCTGGGTGCTgcgc                      >  8:179870/3‑147 (MQ=255)
                                                                                                                                cccGACCTGTCGGCCTGACGCCAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCGGGTGAGTTCATCCCGATTGCCGAACGCACCGGCCTGATTGCCGAGATCGATACCTGGGTGCTgcgc                      >  8:267676/3‑147 (MQ=255)
                                                                                                                                  cGACCTGTCGGCCTGACGCCAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCGGGTGAGTTCATCCCGATTGCCGAACGCACCGGCCTGATTGCCGAGATCGATACCTGGGTGCTGCGCCAgg                  <  7:267676/149‑1 (MQ=255)
                                                                                                                                      cTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCGGGTGAGTTCATCCCGATTGCCGAACGCACCGGCCTGATTGCCGAGATCGATACCTGGGTGCTGCGCCAGGCCTg              >  4:153001/1‑149 (MQ=255)
                                                                                                                                         gCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCGGGTGAGTTCATCCCGATTGCCGAACGCACCGGCCT‑ATTGCCGAGATCGATACCTGGGTGCTGCGCCAGGCCTGCCgg          >  8:52604/1‑149 (MQ=255)
                                                                                                                                           gTCAGGCCGACAGGTCGGGGTTGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCGGGTGAGTTCATCCCGATTGCCGAACGCACCGGCCTGATTGCCGAGATCGATACCTGGGTGCTGCGTCAGGCCGGCAGGc         >  4:122311/1‑149 (MQ=255)
                                                                                                                                                gCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCGGGTGAGTTCATCCCGATTGCCGAACGCACCGGCCTGATTGCCGAGATCGATACCTGGGTGCTGCGCCAGGCCTGCCGGCAAATg    >  2:354480/1‑149 (MQ=255)
                                                                                                                                                 ccGACAGGTCGGGGTTGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCGGGTGAGTTCATCCCGATTGCCGAACGCACCGGCCTGATTGCCGAGATCGATACCTGGGTGCTGCGTCAGGCCGGCAGGcctggcg   <  3:122311/149‑7 (MQ=255)
                                                                                                                                                  cGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCGGGTGAGTTCATCCCGATTGCCGAACGCACCGGCCTGATTGCCGAGATCGATACCTGGGTGCTGCGCCAGGCCTGCCGGCAAATGGt  >  8:216471/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ACGGGCGGGCCTGCTATGCGCTGTACACCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTCGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCGGGTGAGTTCATCCCGATTGCCGAACGCACCGGCCTGATTGCCGAGATCGATACCTGGGTGCTGCGCCAGGCCTGCCGGCAAATGGT  >  NC_002947/271187‑271481

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