Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 271,316 A→G 5.5% H406R (CAC→CGC)  PP_0218 → sensory box protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947271,3160AG 5.5% 217.5 / 2.0 73H406R (CAC→CGC) PP_0218sensory box protein
Reads supporting (aligned to +/- strand):  ref base A (36/33);  new base G (2/2);  total (38/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 = 7.00e-01

TACAAGGCCAAAAGCAACGGGCGGGCCTGCTATGCGCTGTACACCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCG‑GGTGAGTTCATCCCGATTGCCGAACGCACCGGCCTGATTGCCGAGATCGATACCTGGGTGCTGCGCCAG  >  NC_002947/271171‑271464
                                                                                                                                                 |                                                                                                                                                      
tACAAGGCCAAAAGCAACGGGCGGGCCTGCTATGCGCTGTACACCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGa                                                                                                                                                     <  5:335705/149‑1 (MQ=255)
tACAAGGCCAAAAGCAACGGGCGGGCCTGCTATGCGCTGTACACCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGa                                                                                                                                                     >  3:73688/1‑149 (MQ=255)
 aCAAGGCCAAAAGCAACGGGCGGGCCTGCTATGCGCTGTACACCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGAc                                                                                                                                                    >  4:223859/1‑149 (MQ=255)
  cAAGGCCAAAAGCAACGGGCGGGCCTGCTATGCGCTGTACACCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGAcc                                                                                                                                                   >  6:2799/1‑149 (MQ=255)
   aaGGCCAAAAGCAACGGGCGGGCCTGCTATGCGCTGTACACCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCt                                                                                                                                                  >  8:263597/1‑149 (MQ=255)
   aaGGCCAAAAGCAACGGGCGGGCCTGCTATGCGCTGTACACCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAAATGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCt                                                                                                                                                  >  5:61558/1‑149 (MQ=255)
      gCCAAAAGCAACGGGCGGGCCTGCTATGCGCTGTACACCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGc                                                                                                                                               <  3:223859/149‑1 (MQ=255)
       ccAAAAGCAACGGGCGGGCCTGCTATGCGCTGTACACCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCg                                                                                                                                              >  2:377043/1‑149 (MQ=255)
        cAAAAGCAACGGGCGGGCCTGCTATGCGCTGTACACCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGt                                                                                                                                             <  7:209843/149‑1 (MQ=255)
        cAAAAGCAACGGGCGGGCCTGCTATGCGCTGTACACCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGCACTTCCAGCCGGTACACGACCTGGCGt                                                                                                                                             >  5:295725/1‑149 (MQ=255)
            aGCAACGGGCGGGCCTGCTATGCGCTGTACACCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAgg                                                                                                                                         >  7:198185/1‑149 (MQ=255)
             gCAACGGGCGGGCCTGCTATGCGCTGTACACCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGcctgacgc                                                                                                                                        <  1:227408/149‑8 (MQ=255)
              cAACGGGCGGGCCTGCTATGCGCTGTACACCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGcc                                                                                                                                       >  3:348832/1‑149 (MQ=255)
                aCGGGCGGGCCTGCTATGCGCTGTACACCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGa                                                                                                                                     <  7:328968/149‑1 (MQ=255)
                 cgggcgggCCTGCTATGCGCTGTACACCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGAc                                                                                                                                    <  4:73688/149‑1 (MQ=255)
                  gggcgggcCTGCTATGCGCTGTACACCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGCATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGa                                                                                                                                     >  2:375027/1‑147 (MQ=255)
                     caggcCTGCTATGCGCTGTACACCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGt                                                                                                                                <  2:233557/147‑1 (MQ=255)
                            gCTATGCGCTGTACACCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTc                                                                                                                        >  1:398640/1‑149 (MQ=255)
                             cTATGCGCTGTACACCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTc                                                                                                                        <  4:348832/148‑1 (MQ=255)
                             cTATGCGCTGTACACCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGCACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTc                                                                                                                        <  6:295725/148‑1 (MQ=255)
                                  cgcTGTACACCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACg                                                                                                                  >  4:397583/1‑149 (MQ=255)
                                       tACCCCGAAGAGCTGACCGCCCACGCCCTGCATCGGGGCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGCGTCGAGACGCTg                                                                                                               <  3:54713/147‑1 (MQ=255)
                                          acCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACGCGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTgcg                                                                                                          >  5:74887/1‑149 (MQ=255)
                                          acCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACGCGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTgcg                                                                                                          >  1:379686/1‑149 (MQ=255)
                                            cGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCt                                                                                                        >  4:200499/1‑149 (MQ=255)
                                             gAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTg                                                                                                       <  5:2799/149‑1 (MQ=255)
                                               agagCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGgc                                                                                                     <  7:263597/149‑1 (MQ=255)
                                               agagCGGACCGCCCACGCCCATCATCGGGTCGATTCCGCCGCTGAATTGCGTTCGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGgc                                                                                                     <  1:332695/149‑1 (MQ=255)
                                                   cTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGgcagca                                                                                                 >  2:123057/1‑149 (MQ=255)
                                                    tGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCAc                                                                                                <  6:61558/149‑1 (MQ=255)
                                                      aCCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGcc                                                                                                                                       >  1:209879/1‑109 (MQ=255)
                                                      aCCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGcc                                                                                                                                       <  2:209879/109‑1 (MQ=255)
                                                      aCCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGATGCTGGTGCGCTGgca                                                                                                    <  1:358189/143‑1 (MQ=255)
                                                      aCCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGATGCTGGTGCGCTGgca                                                                                                    >  2:358189/1‑143 (MQ=255)
                                                           ccacgcccaGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGgc                                                                                         <  3:397583/149‑1 (MQ=255)
                                                            cacgcccaGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGgcg                                                                                        <  2:398640/149‑1 (MQ=255)
                                                                 ccaGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTc                                                                                   >  8:312928/1‑149 (MQ=255)
                                                                  caGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCt                                                                                  >  4:404230/1‑149 (MQ=255)
                                                                     cATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACGCGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGt                                                                               <  6:74887/149‑1 (MQ=255)
                                                                     cATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACGCGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGt                                                                               <  2:379686/149‑1 (MQ=255)
                                                                     cATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTgg                                                                                >  4:346004/1‑148 (MQ=255)
                                                                      aTCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCCCCCCAGGCGCGGGCTGGTa                                                                              >  4:82252/1‑149 (MQ=255)
                                                                       tCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTAc                                                                             <  3:404230/149‑1 (MQ=255)
                                                                           gTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCa                                                                                            <  4:396455/130‑1 (MQ=255)
                                                                           gTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCa                                                                                            >  3:396455/1‑130 (MQ=255)
                                                                              gagaCCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTc                                                                                   >  7:361702/1‑136 (MQ=255)
                                                                              gagaCCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTc                                                                                   <  8:361702/136‑1 (MQ=255)
                                                                              gagaCCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCg‑gg                                                                     >  8:169164/1‑149 (MQ=255)
                                                                                 aCCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCTACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCG‑GGTGa                                                                  <  6:188738/149‑1 (MQ=255)
                                                                                  ccgcctCGGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGATCTGGTACCGCCG‑GGTGa                                                                  <  2:136026/148‑1 (MQ=255)
                                                                                      cgcTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCG‑GGTGAGTTCa                                                             >  5:138710/1‑149 (MQ=255)
                                                                                         tGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCG‑GGTGAGTTCATcc                                                          >  2:293851/1‑149 (MQ=255)
                                                                                            aTTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCG‑GGTGAGTTCATCCCGa                                                       <  1:291872/149‑1 (MQ=255)
                                                                                            aTTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCG‑GGTGAGTTCATCCCGa                                                       <  3:200499/149‑1 (MQ=255)
                                                                                               gCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCG‑GGTGAGTTCATCCCGAtt                                                     <  1:375027/148‑1 (MQ=255)
                                                                                                  tCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCG‑GGTGAGTTCATCCCGATTGCCg                                                 <  6:138710/149‑1 (MQ=255)
                                                                                                     ggCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGTCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCG‑GGTGAGTTCATCCCGATTg                                                    >  5:22787/1‑143 (MQ=255)
                                                                                                     ggCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGTCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCG‑GGTGAGTTCATCCCGATTg                                                    <  6:22787/143‑1 (MQ=255)
                                                                                                         cTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCG‑GGTGAGTTCATCCCGATTGCCGAACGCAc                                          >  5:218384/1‑149 (MQ=255)
                                                                                                            gAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGTAGGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCG‑GGTGAGTTCATCCCGATTGCCGAACGCACCg                                        >  1:121316/1‑149 (MQ=255)
                                                                                                                 ggATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCG‑GGTGAGTTCATCCCGATTGCCGAACGCACCGGCCTGa                                  <  7:312928/149‑1 (MQ=255)
                                                                                                                  gATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCG‑GGTGAGTTCATCCCGATTGCCGAACGCACCGGCCTGAt                                 >  1:244392/1‑149 (MQ=255)
                                                                                                                     gAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCG‑GGTGAGTTCATCCCGATTGCCGAACGCACCGGCCTGATTGc                              >  5:282869/1‑149 (MQ=255)
                                                                                                                            gggTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCGCGGTGAGTTCATCCCGATTGCCGAACGCACCGGCCTGATTGCCGTGAt                        <  1:293851/149‑1 (MQ=255)
                                                                                                                            gggTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCG‑GGTGAGTTCATCCCGATTGCCGAACGCACCGGCCTGATTGCCGAGATc                       <  3:346004/149‑1 (MQ=255)
                                                                                                                              gtgtACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCG‑GGTGAGTTCATCCCGATTGCCGAACGCACCGGCCTGATTGCCGAGATCGa                     >  5:63351/1‑149 (MQ=255)
                                                                                                                                    ttCCAGCCGGTACACGACCTGGCGTCAGGTCGACAGGT‑CGGGGTCGAGACGCTGGTGTGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCG‑GGTGAGTTCATCCCGATTGCCGAACGCACCGGCCTGATTGCCGAGATCGATACCTg               >  8:198478/1‑149 (MQ=255)
                                                                                                                                    ttCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCg‑g                                                                      <  6:173827/94‑1 (MQ=255)
                                                                                                                                    ttCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCg‑g                                                                      >  5:173827/1‑94 (MQ=255)
                                                                                                                                    ttCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCG‑GGTGAGTTCATCCCGATTGCCGAACGCACCGGCCTGATTGCCGAGATCGATACCTg               >  8:106719/1‑149 (MQ=255)
                                                                                                                                    ttCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCG‑GGTGAGTTCATCCCGATTGCCGAACGCACCGGCCTGATTGCCGAGATCGATACCTg               >  8:273401/1‑149 (MQ=255)
                                                                                                                                    ttCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCG‑GGTGAGTTCATCCCGATTGCCGAACGCACCGGCCTGATTGCCGAGATCGATACCTg               >  1:207380/1‑149 (MQ=255)
                                                                                                                                      ccAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCG‑GGTGAGTTCATCCCGATTGCCGAACGCACCGGCCTGATTGCCGAGATCGATACCTgg              <  6:63351/148‑1 (MQ=255)
                                                                                                                                                cacGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCG‑GGTGAGTTCATCCCGATTGCCGAACGCACCGGCCTGATTGCCGAGATCGATACCTGGGTGCTgcgc     >  6:196577/1‑147 (MQ=255)
                                                                                                                                                cacGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCG‑GGTGAGTTCATCCCGATTGCCGAACGCACCGGCCTGATTGCCGAGATCGATACCTGGGTGCTGCGCCa   <  7:169164/149‑1 (MQ=255)
                                                                                                                                                 acGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCG‑GGTGAGTTCATCCCGATTGCCGAACGCACCGGCCTGATTGCCGAGATCGATACCTGGGTGCTGCGCCAg  >  6:317858/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                      
TACAAGGCCAAAAGCAACGGGCGGGCCTGCTATGCGCTGTACACCGAAGAGCTGACCGCCCACGCCCAGCATCGGGTCGAGACCGCCGCTGAATTGCGTCGGGCCCTGGAGCAGGATGAACTGCGGGTGTACTTCCAGCCGGTACACGACCTGGCGTCAGGCCGACAGGT‑CGGGGTCGAGACGCTGGTGCGCTGGCAGCACCCCAGGCGCGGTCTGGTACCGCCG‑GGTGAGTTCATCCCGATTGCCGAACGCACCGGCCTGATTGCCGAGATCGATACCTGGGTGCTGCGCCAG  >  NC_002947/271171‑271464

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: ‑