Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,503,465 G→T 7.4% E90* (GAG→TAG)  rplM → 50S ribosomal protein L13

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,503,4650GT 7.4% 213.0 / 5.9 68E90* (GAG→TAG) rplM50S ribosomal protein L13
Reads supporting (aligned to +/- strand):  ref base G (34/29);  new base T (2/3);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01

CAGAATACACCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCCACACACTGCTCAGCAGCCTCAA  >  NC_002947/1503322‑1503611
                                                                                                                                               |                                                                                                                                                   
cAGAATACACCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCTCGAAg                                                                                                                                                >  8:154230/1‑149 (MQ=255)
 aGAATACACCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGt                                                                                                                                               <  8:301099/149‑1 (MQ=255)
  gAATACACCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGtt                                                                                                                                              <  4:221956/149‑1 (MQ=255)
       cacCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCg                                                                                                                                         >  8:82385/1‑149 (MQ=255)
       cacCCCTCACGGTGACCCCGGCTACTACATCGTCGTTATCAACGCCGAGCATGTTCGTGTTAATGTTGCCAAGTCTTCCCACAAAATATACTACACCCACTCCGGTTTCCCATGGCGTATCAAGTAAATCAACTTCTAGAAGTTGAttt                                                                                                                                         <  6:247761/149‑3 (MQ=255)
        acCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGc                                                                                                                                        >  3:271479/1‑149 (MQ=255)
        acCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACACCGGTTTCCCAGGCGGGATGAAGGAAATCAACGTCGAGAAGGTGATCGc                                                                                                                                        >  8:279168/1‑149 (MQ=255)
         ccccTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCTCGAAGTTGAtttcc                                                                                                                                       <  7:154230/149‑5 (MQ=255)
          cccTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGcc                                                                                                                                       >  1:165739/1‑148 (MQ=255)
           acTCACGTTGACACCGGCGACTACATCGTCGTTATCAACTCCGAGCAGGTTCGTGTGAGTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCaa                                                                                                                                     <  5:335875/148‑1 (MQ=255)
            cTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAg                                                                                                                                    >  1:138889/1‑149 (MQ=255)
             tCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAgg                                                                                                                                   >  3:136087/1‑149 (MQ=255)
               aCGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGcc                                                                                                                                 >  2:293590/1‑149 (MQ=255)
                 gTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCTCGAAGTTGAtttcctt                                                                                                                                     <  7:68579/143‑7 (MQ=255)
                 gTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCTCGAAGTTGAtttcctt                                                                                                                                     >  8:68579/1‑137 (MQ=255)
                     acacCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGa                                                                                                                            <  2:272644/148‑1 (MQ=255)
                       acCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCg                                                                                                                         >  4:77975/1‑149 (MQ=255)
                        ccGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCgt                                                                                                                        <  3:77975/149‑1 (MQ=255)
                         cGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCgtg                                                                                                                       <  4:271479/149‑1 (MQ=255)
                         cGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCgtg                                                                                                                       >  7:160694/1‑149 (MQ=255)
                              aCTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATc                                                                                                                  <  1:3110/149‑1 (MQ=255)
                                    tcgtcgTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAgg                                                                                                                                   <  8:123658/126‑1 (MQ=255)
                                    tcgtcgTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACt                                                                                                            <  7:276277/149‑1 (MQ=255)
                                       tcgtTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑Cg                                                                                                        >  5:5565/1‑149 (MQ=255)
                                       tcgtTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑Cg                                                                                                        <  7:82385/149‑1 (MQ=255)
                                       tcgtTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGAACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑Cg                                                                                                        >  3:123576/1‑149 (MQ=255)
                                                    cGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCAtgc                                                                                           >  1:188471/1‑149 (MQ=255)
                                                         aGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCt                                                                                      <  2:188471/149‑1 (MQ=255)
                                                         aGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCt                                                                                      <  2:112159/149‑1 (MQ=255)
                                                            ttCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCt                                                                                      >  2:17544/1‑146 (MQ=255)
                                                                        gTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGt                                                                       >  6:292514/1‑149 (MQ=255)
                                                                        gTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGt                                                                       >  2:255034/1‑149 (MQ=255)
                                                                           ccAAGTCTTCCGACAAAATGTACTACTCGCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTcgc                                                                    >  1:192088/1‑149 (MQ=255)
                                                                           ccAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTcgc                                                                    >  4:17324/1‑149 (MQ=255)
                                                                            cAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTcgc                                                                    >  2:37657/1‑148 (MQ=255)
                                                                             aaGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCGCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGa                                                                  <  1:293590/149‑1 (MQ=255)
                                                                             aaGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGa                                                                  <  1:37657/149‑1 (MQ=255)
                                                                             aaGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGa                                                                  <  8:160694/149‑1 (MQ=255)
                                                                              aGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGAc                                                                 <  1:17544/149‑1 (MQ=255)
                                                                                  ttCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGACCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATg                                                              >  7:216777/1‑148 (MQ=255)
                                                                                      gACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCg                                                         >  5:102773/1‑149 (MQ=255)
                                                                                        cAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGc                                                        >  8:11450/1‑148 (MQ=255)
                                                                                         aaaaTGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGc                                                        >  3:310811/1‑147 (MQ=255)
                                                                                         aaaaTGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGc                                                        >  1:17559/1‑147 (MQ=255)
                                                                                           aaTGTAAAACTACCACTCCGGTTTCCCAGGCGGTATCAAGGAAATAAACTTCGATAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCCTGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAg                                                     <  1:255034/149‑1 (MQ=255)
                                                                                             tGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGc                                                    >  2:147307/1‑147 (MQ=255)
                                                                                             tGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGc                                                    <  1:147307/147‑1 (MQ=255)
                                                                                                 ctactCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGaaa                                               <  7:216581/148‑1 (MQ=255)
                                                                                                  tactCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGTCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAAc                                                                                >  7:218458/1‑114 (MQ=255)
                                                                                                  tactCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAAc                                                                                <  8:218458/114‑1 (MQ=255)
                                                                                                     tCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTa                                          >  6:238983/1‑149 (MQ=255)
                                                                                                     tCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTa                                          >  3:136590/1‑149 (MQ=255)
                                                                                                       ccACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAgtgt                                           <  6:69985/146‑1 (MQ=255)
                                                                                                       ccACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAgtgt                                           >  5:69985/1‑146 (MQ=255)
                                                                                                       ccACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTAcg                                        >  7:128835/1‑149 (MQ=255)
                                                                                                          cTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCgg                                     <  4:77135/149‑1 (MQ=255)
                                                                                                          cTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTCAAAGTGTACGCgg                                     >  1:323741/1‑149 (MQ=255)
                                                                                                             cGGTTTCCCAGGCGGTATGAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGTGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGtg                                  <  7:279168/149‑1 (MQ=255)
                                                                                                             cGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGt                                   <  1:135645/148‑1 (MQ=255)
                                                                                                                    ccAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCa                           <  8:128835/149‑1 (MQ=255)
                                                                                                                       ggCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCAccc                        <  4:136087/149‑1 (MQ=255)
                                                                                                                        gCGGTATCAAGGTAATCAACTTCGACAAGTTGATCGCCAACGCCCCTGAGCGTTTTATCGAAACTT‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCca                       <  2:138889/149‑1 (MQ=255)
                                                                                                                           gTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCcaca                     >  1:47802/1‑148 (MQ=255)
                                                                                                                            tATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCcaca                     <  6:251056/147‑1 (MQ=255)
                                                                                                                            tATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCcaca                     >  5:251056/1‑147 (MQ=255)
                                                                                                                                aaGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCCACACACTGc               >  4:205034/1‑149 (MQ=255)
                                                                                                                                           cTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGt                                   >  4:297645/1‑118 (MQ=255)
                                                                                                                                              cGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCCACACACTGCTCAGCAGCCTCaa  <  3:17324/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CAGAATACACCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTG‑CGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCCACACACTGCTCAGCAGCCTCAA  >  NC_002947/1503322‑1503611

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