Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 996,028 T→G 5.5% Q523P (CAA→CCA)  PP_0860 ← sulfite reductase flavoprotein subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947996,0280TG 5.5% 292.0 / 7.5 91Q523P (CAA→CCA) PP_0860sulfite reductase flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base T (41/45);  new base G (3/2);  total (44/47)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.69e-01

AGGCCAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGCGAGCAGGGCAT  >  NC_002947/995887‑996172
                                                                                                                                             |                                                                                                                                                
aGGCCAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGCCACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAg                                                                                                                                           >  8:388596/1‑149 (MQ=255)
  gCCAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGAGCACCAGTGACCAGTTGGCAAAGTTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAggg                                                                                                                                         <  7:347593/149‑1 (MQ=255)
    cAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGggcg                                                                                                                                       >  5:102336/1‑149 (MQ=255)
     aaGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCAGCAATGGCAGCAggagct                                                                                                                                      <  6:321208/149‑5 (MQ=255)
       gCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGC‑AAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCt                                                                                                                                   >  4:160812/1‑149 (MQ=255)
         aGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTg                                                                                                                                  >  6:162172/1‑149 (MQ=255)
         aGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTg                                                                                                                                  >  5:278570/1‑149 (MQ=255)
         aGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTg                                                                                                                                  >  6:118956/1‑149 (MQ=255)
          ggTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTgg                                                                                                                                 <  7:293975/149‑1 (MQ=255)
           gTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGt                                                                                                                                <  7:340276/149‑1 (MQ=255)
               acCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGC‑AAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCg                                                                                                                           <  3:160812/149‑1 (MQ=255)
                ccGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCg                                                                                                                           <  2:198382/149‑1 (MQ=255)
                ccGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCg                                                                                                                           <  6:102336/149‑1 (MQ=255)
                   gTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTAGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTg                                                                                                                        <  2:386523/149‑1 (MQ=255)
                      ggCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTc                                                                                                                      >  7:107607/1‑148 (MQ=255)
                      ggCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTc                                                                                                                      <  8:166269/148‑1 (MQ=255)
                       gCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCa                                                                                                                    <  6:278570/149‑1 (MQ=255)
                          ttGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACtt                                                                                                                 >  5:62656/1‑149 (MQ=255)
                               tGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCAcc                                                                                                            >  3:169654/1‑149 (MQ=255)
                                  cGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCAGCAATGGCAGCAgga                                                                                                                                         >  4:281925/1‑116 (MQ=255)
                                  cGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCAGCAATGGCAGCAgga                                                                                                                                         <  3:281925/117‑2 (MQ=255)
                                    gggTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTCCCACCGGGCt                                                                                                       >  1:52767/1‑149 (MQ=255)
                                     ggTTGAGCAGTTCGCGGCGCACCAGTGACTAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCAtt                                                                                                                                                  <  2:183122/105‑1 (MQ=255)
                                     ggTTGAGCAGTTCGCGGCGCACCAGTGACTAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCAtt                                                                                                                                                  >  1:183122/1‑105 (MQ=255)
                                       ttGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGaa                                                                                                    >  1:200530/1‑149 (MQ=255)
                                         gAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAAt                                                                                                   <  5:118956/148‑1 (MQ=255)
                                         gAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAAt                                                                                                   <  4:203260/148‑1 (MQ=255)
                                            ccggtCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCAc                                                                                               <  8:92393/145‑1 (MQ=255)
                                            cAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCAc                                                                                               <  6:252366/149‑1 (MQ=255)
                                             aGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCAct                                                                                              <  6:62656/149‑1 (MQ=255)
                                               ttCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCActgc                                                                                            >  6:324025/1‑149 (MQ=255)
                                               ttCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCActgc                                                                                            >  3:86864/1‑149 (MQ=255)
                                                 cgcgGCGCACCAGTGACCAGTTGGCATAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGtt                                                                                                                         <  3:142317/118‑1 (MQ=255)
                                                 cgcgGCGCACCAGTGACCAGATGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGACAGATGCGCCAGCTCCTGCTGCCATTGCTGCAGCGCCGCCTGGTCCGCGtt                                                                                                                         >  4:142317/1‑118 (MQ=255)
                                                  gcggcgCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTgg                                                                                         <  8:210333/149‑1 (MQ=255)
                                                  gcggcgCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTgg                                                                                         >  5:224826/1‑149 (MQ=255)
                                                    ggcgCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGTGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCa                                                                                       >  4:2716/1‑149 (MQ=255)
                                                       gcACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACcg                                                                                    >  3:174495/1‑149 (MQ=255)
                                                           cAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGct                                                                                >  3:300793/1‑149 (MQ=255)
                                                             gTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGctct                                                                              >  3:325141/1‑149 (MQ=255)
                                                                aCCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGTGCTGAAt                                                                                                   <  8:158251/125‑1 (MQ=255)
                                                                aCCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAAt                                                                                                   >  7:158251/1‑125 (MQ=255)
                                                                 ccAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTcc                                                                            >  7:307699/1‑147 (MQ=255)
                                                                 ccAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCa                                                                          <  5:324025/149‑1 (MQ=255)
                                                                   aGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGc                                                                        <  4:169654/149‑1 (MQ=255)
                                                                        gCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGc                                                                   <  8:307699/149‑1 (MQ=255)
                                                                          aaaGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGCATGCACTGCTGGCACCGCGCTCTCCCCGCCATGCCt                                                                 >  4:175015/1‑149 (MQ=255)
                                                                          aaaGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAAt                                                                                                   >  3:198031/1‑115 (MQ=255)
                                                                          aaaGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAAt                                                                                                   <  4:198031/115‑1 (MQ=255)
                                                                             ggTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGaa                                                              <  4:300793/149‑1 (MQ=255)
                                                                             ggTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGaa                                                              <  5:162172/149‑1 (MQ=255)
                                                                              gTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAg                                                             <  4:174495/149‑1 (MQ=255)
                                                                                   cggcTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGcgac                                                        <  1:67889/149‑1 (MQ=255)
                                                                                    ggcTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGcgac                                                        >  7:30482/1‑148 (MQ=255)
                                                                                        gCCAGCCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCAACTTCAACCGGGCTGAATGCACTGCTGGCTCCGCGCTCTCCCAGCCATGCCTGAAGGcgacgatc                                                    >  7:185988/1‑146 (MQ=255)
                                                                                         ccAGGTGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGGAGGCGACGAGc                                                    <  3:2716/147‑1 (MQ=255)
                                                                                           aGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGc                                                                   <  8:14576/130‑1 (MQ=255)
                                                                                           aGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGc                                                                   >  7:14576/1‑130 (MQ=255)
                                                                                               ggcCACTGGCTGGACACCGGGCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGc                                            <  3:314585/149‑1 (MQ=255)
                                                                                                gcCTCTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCa                                                                                                >  4:60649/1‑96 (MQ=255)
                                                                                                gcCTCTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCa                                                                                                <  3:60649/96‑1 (MQ=255)
                                                                                                       gctggACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGc                                    >  5:36096/1‑149 (MQ=255)
                                                                                                       gctggACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGc                                    <  6:224826/149‑1 (MQ=255)
                                                                                                        ctggACTCCGTTCAGTTGCGACAGCTCCTGCTGACATTGCTGCAGGGCGGCCTTGTCGGCGTTGTCCACTTCCACCTGGCTGAATGCACTGCTTGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCg                                   <  8:185988/149‑1 (MQ=255)
                                                                                                        ctggACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCg                                   >  1:92786/1‑149 (MQ=255)
                                                                                                        ctggACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCg                                   >  6:392628/1‑149 (MQ=255)
                                                                                                        ctggACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCg                                   <  1:176035/149‑1 (MQ=255)
                                                                                                              acCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACt                             >  1:23085/1‑149 (MQ=255)
                                                                                                              acCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACt                             >  1:320485/1‑149 (MQ=255)
                                                                                                                cGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGa                           <  8:30482/149‑1 (MQ=255)
                                                                                                                cGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGACCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGaa                                                              >  3:372275/1‑114 (MQ=255)
                                                                                                                cGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGACCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGaa                                                              <  4:372275/114‑1 (MQ=255)
                                                                                                                  gTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGAc                          >  4:141456/1‑148 (MQ=255)
                                                                                                                   tCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGa                        <  8:107607/149‑1 (MQ=255)
                                                                                                                       ttGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGAAGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCAc                    >  5:215734/1‑149 (MQ=255)
                                                                                                                          cgcCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCAACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTTAAGGCGACGAGCGAAACAGAAGAAGTGCGGGTACTGACGATCAcc                   >  8:160074/1‑147 (MQ=255)
                                                                                                                             cAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAggg              >  2:367774/1‑149 (MQ=255)
                                                                                                                             cAGCTCCTGCTGCCATTGCTGAAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAggg              <  7:160074/149‑1 (MQ=255)
                                                                                                                               gCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGTCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGa                           >  6:290933/1‑134 (MQ=255)
                                                                                                                               gCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGTCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGa                           <  5:290933/134‑1 (MQ=255)
                                                                                                                               gCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGCg            <  2:200530/149‑1 (MQ=255)
                                                                                                                               cgcccTGCAGCAATGGCAGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGc             >  6:217366/4‑148 (MQ=255)
                                                                                                                                  ccTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGCGAGc         <  7:388596/149‑1 (MQ=255)
                                                                                                                                  ccTGCAGCAATGGCAGCAGGGCGGCCTGGTCGGCGTTGTCCCCTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCTGAAGTGCGGGTCCTGACGATCACCGAGGGCTAg          >  3:308067/1‑148 (MQ=255)
                                                                                                                                    tgctgcCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGCGAGCAg       >  1:244526/1‑149 (MQ=255)
                                                                                                                                    tgctgcCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGCGAGCAg       >  8:371319/1‑149 (MQ=255)
                                                                                                                                      ctgcCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGTGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGCGAGCAg       <  4:86864/147‑1 (MQ=255)
                                                                                                                                      ctgcCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGCGAGCAgg      <  2:92786/148‑1 (MQ=255)
                                                                                                                                      ctgcCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGCGAGCAgg      <  6:36096/148‑1 (MQ=255)
                                                                                                                                         ccATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGCGAGCAGGGCa   <  3:141456/148‑1 (MQ=255)
                                                                                                                                         ccATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGCGAGCAGGGCAt  <  3:175015/149‑1 (MQ=255)
                                                                                                                                           aTTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAg                                          <  1:277764/107‑1 (MQ=255)
                                                                                                                                           aTTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAg                                          >  2:277764/1‑107 (MQ=255)
                                                                                                                                             |                                                                                                                                                
AGGCCAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGCGAGCAGGGCAT  >  NC_002947/995887‑996172

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