Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 996,025 C→A 7.5% W524L (TGG→TTG)  PP_0860 ← sulfite reductase flavoprotein subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947996,0250CA 7.5% 214.3 / 6.3 67W524L (TGG→TTG) PP_0860sulfite reductase flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base C (34/28);  new base A (3/2);  total (37/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.80e-01

CAGCCTGCAGGCCAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGCGAGCAGGGCA  >  NC_002947/995879‑996171
                                                                                                                                                   |                                                                                                                                                  
cAGCCTGCAGGCCAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCAGCAAt                                                                                                                                                  >  5:255468/1‑149 (MQ=255)
     tGCAGGCCAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATtgctg                                                                                                                                             >  5:222182/1‑149 (MQ=255)
       cAGGCCAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATtgctg                                                                                                                                             <  3:16021/147‑1 (MQ=255)
          gCCAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAggg                                                                                                                                        >  1:233601/1‑149 (MQ=255)
          gCCAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCAGCAATGGCAGCAgga                                                                                                                                        <  5:82742/149‑2 (MQ=255)
          gCCAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCAGCAATGGCAGCAgga                                                                                                                                        <  6:255468/149‑2 (MQ=255)
            cAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGggcg                                                                                                                                      <  6:222182/149‑1 (MQ=255)
             aaGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGggcgg                                                                                                                                     >  6:176363/1‑149 (MQ=255)
             aaGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGggcgg                                                                                                                                     >  3:86312/1‑149 (MQ=255)
                     acacCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCg                                                                                                                             >  3:20336/1‑149 (MQ=255)
                       acCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGc                                                                                                                           >  5:343948/1‑149 (MQ=255)
                       acCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGc                                                                                                                           <  7:124273/149‑1 (MQ=255)
                          ggTTGGCCTTGGCTGCCGTGGTTGAGCAGTTCGCGGCGCGCAAGAGACCAGTGGGCAAAGGTTGGCGGCTCGC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGGCATTGATGCAGGGCGGCCTGGTCGGAGtt                                                                                                                        >  6:86980/1‑149 (MQ=255)
                           gTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTg                                                                                                                       >  3:171626/1‑149 (MQ=255)
                           gTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTg                                                                                                                       <  2:233601/149‑1 (MQ=255)
                               gCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTcc                                                                                                                    <  2:308603/148‑1 (MQ=255)
                                ccTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCAc                                                                                                                  >  3:116223/1‑149 (MQ=255)
                                 cTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACt                                                                                                                 >  4:39143/1‑149 (MQ=255)
                                              gTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGa                                                                                                    <  8:8904/149‑1 (MQ=255)
                                                  aGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATg                                                                                                 <  6:343948/148‑1 (MQ=255)
                                                      gTTCGCGGCCCACCAGTGCCCAGTTGGCAAAGGTTGGCGGCTGCC‑AAGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGCATGCActg                                                                                            >  6:282471/1‑149 (MQ=255)
                                                       ttCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCGAGGCGCCCACTGGCTGGACACCGGTAATTTGCCCCAGATCCTCCTCCCAt                                                                                                                                                  >  1:76069/1‑95 (MQ=255)
                                                              cgcACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCAcc                                                                                    <  3:203394/149‑1 (MQ=255)
                                                                             ttGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCa                                                                     >  8:328333/1‑149 (MQ=255)
                                                                             ttGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCa                                                                     >  5:230928/1‑149 (MQ=255)
                                                                                gCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATg                                                                   <  6:230928/148‑1 (MQ=255)
                                                                                gCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGc                                                                  >  3:248502/1‑149 (MQ=255)
                                                                                  aaaGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCt                                                                >  7:62725/1‑149 (MQ=255)
                                                                                  aaaGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCt                                                                <  3:39143/149‑1 (MQ=255)
                                                                                     ggTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGaa                                                             <  8:62725/149‑1 (MQ=255)
                                                                                           cggcTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGcgac                                                       <  5:176363/149‑1 (MQ=255)
                                                                                            ggcTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGcgac                                                       >  8:95857/1‑148 (MQ=255)
                                                                                             gcTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGt                                                                                                                      >  8:235070/1‑84 (MQ=255)
                                                                                             gcTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGt                                                                                                                      <  7:235070/84‑1 (MQ=255)
                                                                                                gCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGc                                                   >  5:68913/1‑148 (MQ=255)
                                                                                                 cc‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGa                                                 <  7:328333/149‑1 (MQ=255)
                                                                                                    aGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGaaa                                               <  4:248502/149‑1 (MQ=255)
                                                                                                       cggcCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGc                                                   >  5:261506/1‑142 (MQ=255)
                                                                                                       cggcCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGc                                                   <  6:261506/142‑1 (MQ=255)
                                                                                                           cACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAg                                         >  1:206737/1‑148 (MQ=255)
                                                                                                                gctggACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTg                                    <  7:263443/148‑1 (MQ=255)
                                                                                                                    gACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCg                                            <  4:172168/136‑1 (MQ=255)
                                                                                                                    gACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCg                                            >  3:172168/1‑136 (MQ=255)
                                                                                                                    gACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCggg                                >  3:99665/1‑148 (MQ=255)
                                                                                                                    gACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCgg                                 <  2:206737/147‑1 (MQ=255)
                                                                                                                       acCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACt                            >  2:223565/1‑149 (MQ=255)
                                                                                                                         cGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCgg                                 >  3:51161/1‑142 (MQ=255)
                                                                                                                         cGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCgg                                 <  4:51161/142‑1 (MQ=255)
                                                                                                                         cGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGa                          <  5:282471/149‑1 (MQ=255)
                                                                                                                         cGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGa                          <  6:68913/149‑1 (MQ=255)
                                                                                                                         cGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTAGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGa                          >  8:12555/1‑149 (MQ=255)
                                                                                                                                ttGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGc                                                                                                >  2:222235/1‑72 (MQ=255)
                                                                                                                                ttGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGc                                                                                                <  1:222235/72‑1 (MQ=255)
                                                                                                                                ttGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCAc                   >  7:228805/1‑149 (MQ=255)
                                                                                                                                   cgcCAGCTCCTTCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCATGAAGGCGACTAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGa                <  5:86980/149‑1 (MQ=255)
                                                                                                                                      cAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAggg             >  4:223893/1‑149 (MQ=255)
                                                                                                                                      cAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAggg             >  3:350535/1‑149 (MQ=255)
                                                                                                                                         gcccTGCAGCAATGGCAGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGc            >  5:62193/3‑147 (MQ=255)
                                                                                                                                           ccTGCAGCAATGGCAGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGCGAg         >  6:241068/1‑148 (MQ=255)
                                                                                                                                            ctgctgCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGc            <  7:305401/144‑1 (MQ=255)
                                                                                                                                            ctgctgCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGc            >  8:305401/1‑144 (MQ=255)
                                                                                                                                            ctgctgCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGCGAg         >  6:82758/1‑147 (MQ=255)
                                                                                                                                               ctgcCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGCGAGCAgg     <  8:228805/148‑1 (MQ=255)
                                                                                                                                               ctgcCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGCGAGCAg      <  4:99665/147‑1 (MQ=255)
                                                                                                                                                 gcCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATc                     <  7:212994/130‑1 (MQ=255)
                                                                                                                                                 gcCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATc                     >  8:212994/1‑130 (MQ=255)
                                                                                                                                                 gcCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGCGAGCAGGGCa  >  2:167112/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CAGCCTGCAGGCCAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCC‑AGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGCGAGCAGGGCA  >  NC_002947/995879‑996171

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