Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 996,031 T→A 6.3% Q522L (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,0310TA 6.3% 190.4 / 3.0 63Q522L (CAG→CTG) PP_0860sulfite reductase flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base T (32/27);  new base A (2/2);  total (34/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

TGCAGGCCAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGCGAGCAGGGCATAGTTGA  >  NC_002947/995884‑996178
                                                                                                                                                   |                                                                                                                                                   
tGCAGGCCAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATtgctg                                                                                                                                                    >  5:222182/1‑149 (MQ=255)
  cAGGCCAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATtgctg                                                                                                                                                    <  3:16021/147‑1 (MQ=255)
     gCCAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAggg                                                                                                                                               >  1:233601/1‑149 (MQ=255)
     gCCAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCAGCAATGGCAGCAgga                                                                                                                                               <  5:82742/149‑2 (MQ=255)
     gCCAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCAGCAATGGCAGCAgga                                                                                                                                               <  6:255468/149‑2 (MQ=255)
       cAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGggcg                                                                                                                                             <  6:222182/149‑1 (MQ=255)
        aaGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGggcgg                                                                                                                                            >  3:86312/1‑149 (MQ=255)
        aaGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGggcgg                                                                                                                                            >  6:176363/1‑149 (MQ=255)
                acacCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCg                                                                                                                                    >  3:20336/1‑149 (MQ=255)
                  acCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGc                                                                                                                                  >  5:343948/1‑149 (MQ=255)
                  acCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGc                                                                                                                                  <  7:124273/149‑1 (MQ=255)
                     ggTTGGCCTTGGCTGCCGTGGTTGAGCAGTTCGCGGCGCGCAAGAGACCAGTGGGCAAAGGTTGGCGGCTCGCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGGCATTGATGCAGGGCGGCCTGGTCGGAGtt                                                                                                                               >  6:86980/1‑149 (MQ=255)
                      gTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTg                                                                                                                              <  2:233601/149‑1 (MQ=255)
                      gTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTg                                                                                                                              >  3:171626/1‑149 (MQ=255)
                          gCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTcc                                                                                                                           <  2:308603/148‑1 (MQ=255)
                           ccTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCAc                                                                                                                         >  3:116223/1‑149 (MQ=255)
                            cTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACt                                                                                                                        >  4:39143/1‑149 (MQ=255)
                                         gTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGa                                                                                                           <  8:8904/149‑1 (MQ=255)
                                             aGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATg                                                                                                        <  6:343948/148‑1 (MQ=255)
                                                 gTTCGCGGCCCACCAGTGCCCAGTTGGCAAAGGTTGGCGGCTGCCAAGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGCATGCActg                                                                                                   >  6:282471/1‑149 (MQ=255)
                                                         cgcACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCAcc                                                                                           <  3:203394/149‑1 (MQ=255)
                                                                        ttGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCa                                                                            >  8:328333/1‑149 (MQ=255)
                                                                        ttGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCa                                                                            >  5:230928/1‑149 (MQ=255)
                                                                           gCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATg                                                                          <  6:230928/148‑1 (MQ=255)
                                                                           gCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGc                                                                         >  3:248502/1‑149 (MQ=255)
                                                                             aaaGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCt                                                                       >  7:62725/1‑149 (MQ=255)
                                                                             aaaGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCt                                                                       <  3:39143/149‑1 (MQ=255)
                                                                                ggTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGaa                                                                    <  8:62725/149‑1 (MQ=255)
                                                                                      cggcTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGcgac                                                              <  5:176363/149‑1 (MQ=255)
                                                                                       ggcTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGcgac                                                              >  8:95857/1‑148 (MQ=255)
                                                                                        gcTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGt                                                                                                                             >  8:235070/1‑84 (MQ=255)
                                                                                        gcTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGt                                                                                                                             <  7:235070/84‑1 (MQ=255)
                                                                                           gCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGc                                                          >  5:68913/1‑148 (MQ=255)
                                                                                            ccAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGa                                                        <  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                                                                                                       <  1:222235/72‑1 (MQ=255)
                                                                                                                          ttGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGc                                                                                                       >  2:222235/1‑72 (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                            >  8:212994/1‑130 (MQ=255)
                                                                                                                                           gcCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATc                            <  7:212994/130‑1 (MQ=255)
                                                                                                                                           gcCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGCGAGCAGGGCa         >  2:167112/1‑149 (MQ=255)
                                                                                                                                                 gctgcAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGCGAGCAGGGCATAGtt    <  4:171626/148‑1 (MQ=255)
                                                                                                                                                  ctgcAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGCGAGCAGGGCATAGTTGa  <  4:350535/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGCAGGCCAAGCAGGTACACCGGTTGGCCTTGGCTGCCGGGGTTGAGCAGTTCGCGGCGCACCAGTGACCAGTTGGCAAAGGTTGGCGGCTGCCAGGCGGCCACTGGCTGGACACCGGTCAGTTGCGCCAGCTCCTGCTGCCATTGCTGCAGGGCGGCCTGGTCGGCGTTGTCCACTTCCACCGGGCTGAATGCACTGCTGGCACCGCGCTCTCCCAGCCATGCCTGAAGGCGACGAGCGAAACCGCAGAAGTGCGGGTACTGACGATCACCGAGGGCGAGCAGGGCATAGTTGA  >  NC_002947/995884‑996178

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