Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,054,316 G→C 5.2% P104A (CCG→GCG)  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,3160GC 5.2% 242.5 / 3.3 78P104A (CCG→GCG) PP_1829alpha/beta family hydrolase
Reads supporting (aligned to +/- strand):  ref base G (34/38);  new base C (2/2);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.18e-01

CCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCGATG  >  NC_002947/2054172‑2054461
                                                                                                                                                |                                                                                                                                                  
ccAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCAGGCTCGTCAGGCGTGGCCggctgg                                                                                                                                                >  6:322690/1‑149 (MQ=255)
 cAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCggctggc                                                                                                                                               <  4:148086/149‑1 (MQ=255)
 cAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCggctggc                                                                                                                                               >  5:188760/1‑149 (MQ=255)
 cAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCggctggc                                                                                                                                               >  7:345159/1‑149 (MQ=255)
    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                                                                                                                                     <  1:399904/149‑1 (MQ=255)
           aTCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGCCAGCCAGCCGGCCAc                                                                                                                                     <  6:219049/149‑1 (MQ=255)
           aTCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGCCAGCCAGCCGGCCAc                                                                                                                                     >  2:399904/1‑149 (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:365825/149‑1 (MQ=255)
                                         aTCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAAt                                                                                                      >  4:289262/1‑149 (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                                                                 >  6:362193/1‑149 (MQ=255)
                                                                              ggCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAgg                                                                 <  3:122363/149‑1 (MQ=255)
                                                                              ggCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGACTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAgg                                                                 <  5:155726/149‑1 (MQ=255)
                                                                                            acaTCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGGTGTCGACCCGGcc                                                   >  8:165585/1‑149 (MQ=255)
                                                                                             caTCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGccc                                                  >  1:42795/1‑149 (MQ=255)
                                                                                              aTCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGccc                                                  <  4:169320/148‑1 (MQ=255)
                                                                                              aTCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGccc                                                  >  3:169320/1‑148 (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                                         >  1:392062/1‑139 (MQ=255)
                                                                                                                cACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCggg                                         <  2:392062/139‑1 (MQ=255)
                                                                                                                cACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGt                               >  5:314602/1‑149 (MQ=255)
                                                                                                                  cGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCa                                                                            <  1:292844/102‑1 (MQ=255)
                                                                                                                  cGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCa                                                                            >  2:292844/1‑102 (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                                                                              >  7:191515/1‑90 (MQ=255)
                                                                                                                            aCGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGc                                                                              <  8:191515/90‑1 (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                 >  1:407474/1‑137 (MQ=255)
                                                                                                                                          tGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTgggg                 <  2:407474/137‑1 (MQ=255)
                                                                                                                                             ccGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCgatg  <  2:42795/149‑1 (MQ=255)
                                                                                                                                             ccGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCgatg  <  1:361549/149‑1 (MQ=255)
                                                                                                                                              cGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTcg                             <  8:297886/121‑1 (MQ=255)
                                                                                                                                              cGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTcg                             >  7:297886/1‑121 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
CCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCGATG  >  NC_002947/2054172‑2054461

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: