Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 996,022 T→A 5.9% Q525L (CAG→CTG)  PP_0860 ← sulfite reductase flavoprotein subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947996,0220TA 5.9% 263.6 / 5.3 85Q525L (CAG→CTG) PP_0860sulfite reductase flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base T (40/40);  new base A (3/2);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.21e-01

TTCTGCAGCCTGCAGGCCAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGCGAGCAGG  >  NC_002947/995874‑996168
                                                                                                                                                    |                                                                                                                                                  
ttCTGCAGCCTGCAGGCCAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCctgct                                                                                                                                                    <  3:236817/149‑1 (MQ=255)
   tGCAGCCTGCAGGCCAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCtgctgc                                                                                                                                                  >  7:210333/1‑148 (MQ=255)
     cAGCCTGCAGGCCAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCtgctgc                                                                                                                                                  <  3:220960/146‑1 (MQ=255)
     cAGCCTGCAGGCCAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCAt                                                                                                                                               >  5:252366/1‑149 (MQ=255)
             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)
                    gCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCt                                                                                                                               >  4:160812/1‑149 (MQ=255)
                      aGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTg                                                                                                                              >  5:278570/1‑149 (MQ=255)
                      aGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTg                                                                                                                              >  6:162172/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)
                            acCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCg                                                                                                                       <  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                                                                                                                                     <  3:281925/117‑2 (MQ=255)
                                               cGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCAGCAATGGCAGCAgga                                                                                                                                     >  4:281925/1‑116 (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                                                                                        >  3:86864/1‑149 (MQ=255)
                                                            ttCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCActgc                                                                                        >  6:324025/1‑149 (MQ=255)
                                                              cgcgGCGCACCAGTGACCAGTTGGCATAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGtt                                                                                                                     <  3:142317/118‑1 (MQ=255)
                                                              cgcgGCGCACCAGTGACCAGATGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGACAGATGCGCCAGCTCCTGCTGCCATTGCTGCAGCGCCGCCTGGTCCGCGtt                                                                                                                     >  4:142317/1‑118 (MQ=255)
                                                               gcggcgCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTgg                                                                                     >  5:224826/1‑149 (MQ=255)
                                                               gcggcgCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTgg                                                                                     <  8:210333/149‑1 (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                                                                                               <  4:198031/115‑1 (MQ=255)
                                                                                       aaaGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAAt                                                                                               >  3:198031/1‑115 (MQ=255)
                                                                                          ggTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGaa                                                          <  5:162172/149‑1 (MQ=255)
                                                                                          ggTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGaa                                                          <  4:300793/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                                                                                            <  3:60649/96‑1 (MQ=255)
                                                                                                             gcCTCTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCa                                                                                            >  4:60649/1‑96 (MQ=255)
                                                                                                                    gctggACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGc                                >  5:36096/1‑149 (MQ=255)
                                                                                                                    gctggACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGc                                <  6:224826/149‑1 (MQ=255)
                                                                                                                     ctggACTCCGTTCAGTTGCGACAGCTCCTGCTGACATTGCTGCAGGGCGGCCTTGTCGGCGTTGTCCACTTCCACCTGGCTGAATGCACTGCTTGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCg                               <  8:185988/149‑1 (MQ=255)
                                                                                                                     ctggACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCg                               >  6:392628/1‑149 (MQ=255)
                                                                                                                     ctggACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCg                               <  1:176035/149‑1 (MQ=255)
                                                                                                                     ctggACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCg                               >  1:92786/1‑149 (MQ=255)
                                                                                                                           acCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACt                         >  1:320485/1‑149 (MQ=255)
                                                                                                                           acCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACt                         >  1:23085/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)
                                                                                                                                                    |                                                                                                                                                  
TTCTGCAGCCTGCAGGCCAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGCGAGCAGG  >  NC_002947/995874‑996168

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: