Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,054,324 G→C 5.4% A101G (GCC→GGC)  PP_1829 ← alpha/beta family hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,054,3240GC 5.4% 244.0 / 3.2 74A101G (GCC→GGC) PP_1829alpha/beta family hydrolase
Reads supporting (aligned to +/- strand):  ref base G (32/38);  new base C (2/2);  total (34/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

GCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCGATGATCAGCAC  >  NC_002947/2054176‑2054469
                                                                                                                                                    |                                                                                                                                                  
gCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGctggctgg                                                                                                                                                    >  8:326270/1‑149 (MQ=255)
    gCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTggcggc                                                                                                                                                >  3:365825/1‑149 (MQ=255)
    gCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTggcg                                                                                                                                                  <  3:163414/147‑1 (MQ=255)
    gCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGACGAAGCGCGGGTCGGCTTTGAGATTGCGGCACCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTggcgcc                                                                                                                                                >  2:361549/1‑147 (MQ=255)
       aTCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGCCAGCCAGCCGGCCAc                                                                                                                                             >  2:399904/1‑149 (MQ=255)
       aTCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGCCAGCCAGCCGGCCAc                                                                                                                                             <  1:399904/149‑1 (MQ=255)
       aTCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGCCAGCCAGCCGGCCAc                                                                                                                                             <  6:219049/149‑1 (MQ=255)
           gTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCcg                                                                                                                                        >  7:378386/1‑149 (MQ=255)
                ccTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCGCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGttt                                                                                                                                   <  8:293428/149‑1 (MQ=255)
                   tCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTc                                                                                                                                <  8:345159/149‑1 (MQ=255)
                    cGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTc                                                                                                                                <  8:238469/148‑1 (MQ=255)
                    cGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTc                                                                                                                                <  5:170234/148‑1 (MQ=255)
                         tgGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAac                                                                                                                           <  4:123274/148‑1 (MQ=255)
                                tCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTg                                                                                                                   >  7:162569/1‑149 (MQ=255)
                                  aGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTgg                                                                                                                  <  1:136704/148‑1 (MQ=255)
                                     aTCAGCGGGCCGAAGCGCGGGTCGTCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAAt                                                                                                              <  2:26590/149‑1 (MQ=255)
                                     aTCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAAt                                                                                                              <  6:10288/149‑1 (MQ=255)
                                     aTCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAAt                                                                                                              <  1:291722/149‑1 (MQ=255)
                                     aTCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAAt                                                                                                              >  4:289262/1‑149 (MQ=255)
                                     aTCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAAt                                                                                                              <  4:365825/149‑1 (MQ=255)
                                     aTCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGAGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCaag                                                                                                              >  6:315950/1‑148 (MQ=255)
                                           gggCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTc                                                                                                        <  2:43635/149‑1 (MQ=255)
                                              ccGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGtt                                                                                                     <  8:378386/149‑1 (MQ=255)
                                               cGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTg                                                                                                    >  6:291759/1‑149 (MQ=255)
                                                gAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTgg                                                                                                   <  6:123430/149‑1 (MQ=255)
                                                gAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTgg                                                                                                   >  6:51789/1‑149 (MQ=255)
                                                 aaGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGc                                                                                                  <  4:413/149‑1 (MQ=255)
                                                  aGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACACCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGc                                                                                                  <  1:78431/149‑1 (MQ=255)
                                                       gggTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCcagcag                                                                                            <  5:51789/149‑1 (MQ=255)
                                                                    gCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGc                                                                               >  1:152588/1‑149 (MQ=255)
                                                                        gCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGcc                                                                                     <  6:277767/139‑1 (MQ=255)
                                                                        gCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGcc                                                                                     >  5:277767/1‑139 (MQ=255)
                                                                        gCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCtt                                                                             >  3:329214/1‑147 (MQ=255)
                                                                        gCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGACTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCtt                                                                             >  6:155726/1‑147 (MQ=255)
                                                                          ggCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAgg                                                                         <  3:122363/149‑1 (MQ=255)
                                                                          ggCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAgg                                                                         >  6:362193/1‑149 (MQ=255)
                                                                          ggCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGACTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAgg                                                                         <  5:155726/149‑1 (MQ=255)
                                                                                        acaTCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGGTGTCGACCCGGcc                                                           >  8:165585/1‑149 (MQ=255)
                                                                                         caTCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGccc                                                          >  1:42795/1‑149 (MQ=255)
                                                                                          aTCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGccc                                                          >  3:169320/1‑148 (MQ=255)
                                                                                          aTCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGccc                                                          <  4:169320/148‑1 (MQ=255)
                                                                                          aTCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGccc                                                          >  1:214969/1‑148 (MQ=255)
                                                                                          aTCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCg                                                         >  5:224136/1‑149 (MQ=255)
                                                                                           tCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGa                                                        >  1:271873/1‑149 (MQ=255)
                                                                                                 acaTAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCgg                                                  <  2:152588/149‑1 (MQ=255)
                                                                                                  caTAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCggg                                                 <  7:326270/149‑1 (MQ=255)
                                                                                                   aTAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGgt                                                >  1:96841/1‑149 (MQ=255)
                                                                                                   aTAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGgt                                                >  8:78644/1‑149 (MQ=255)
                                                                                                        gcTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGgttgaca                                           <  5:291759/149‑4 (MQ=255)
                                                                                                            cACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCggg                                                 <  2:392062/139‑1 (MQ=255)
                                                                                                            cACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCggg                                                 >  1:392062/1‑139 (MQ=255)
                                                                                                            cACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGt                                       >  5:314602/1‑149 (MQ=255)
                                                                                                              cGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCa                                                                                    >  2:292844/1‑102 (MQ=255)
                                                                                                              cGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCa                                                                                    <  1:292844/102‑1 (MQ=255)
                                                                                                              cGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTcg                                     <  2:214969/149‑1 (MQ=255)
                                                                                                               gCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGTCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTcgc                                    >  7:85516/1‑149 (MQ=255)
                                                                                                               gCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTcgc                                    >  7:188153/1‑149 (MQ=255)
                                                                                                                     gTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGcggtcgg                              >  3:326535/1‑149 (MQ=255)
                                                                                                                      tCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGcggtcggt                             <  6:224136/149‑1 (MQ=255)
                                                                                                                      tCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGGTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGcggtcggt                             <  5:322690/149‑1 (MQ=255)
                                                                                                                        aCGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGc                                                                                      <  8:191515/90‑1 (MQ=255)
                                                                                                                        aCGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGc                                                                                      >  7:191515/1‑90 (MQ=255)
                                                                                                                        aCGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGgtt                                               <  2:408757/129‑1 (MQ=255)
                                                                                                                        aCGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTgg                           >  6:125361/1‑149 (MQ=255)
                                                                                                                         cGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGTCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTggg                          <  8:85516/149‑1 (MQ=255)
                                                                                                                          gCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTgggg                         >  5:116536/1‑149 (MQ=255)
                                                                                                                          gCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTgg                           <  7:60599/147‑1 (MQ=255)
                                                                                                                          gCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTgg                           >  8:60599/1‑147 (MQ=255)
                                                                                                                              gTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCg                      >  5:119660/1‑148 (MQ=255)
                                                                                                                              acgcggggTGGCCGCCAGCCAGCCGGCCACGCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCgg                     >  3:158929/8‑149 (MQ=255)
                                                                                                                               tCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGa                    <  8:66466/149‑1 (MQ=255)
                                                                                                                                 aGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAg                   <  3:289262/148‑1 (MQ=255)
                                                                                                                                      tGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTgggg                         <  2:407474/137‑1 (MQ=255)
                                                                                                                                      tGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTgggg                         >  1:407474/1‑137 (MQ=255)
                                                                                                                                         ccGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCgatg          <  2:42795/149‑1 (MQ=255)
                                                                                                                                         ccGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCgatg          <  1:361549/149‑1 (MQ=255)
                                                                                                                                          cGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTcg                                     >  7:297886/1‑121 (MQ=255)
                                                                                                                                          cGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTcg                                     <  8:297886/121‑1 (MQ=255)
                                                                                                                                                gctggcGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCGATGATCAGc    <  2:271873/148‑1 (MQ=255)
                                                                                                                                                 ctggcGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCGATGATCAGcac  <  4:329214/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCGATGATCAGCAC  >  NC_002947/2054176‑2054469

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: