Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 516,694 C→A 10.2% E100* (GAA→TAA) ‡ PP_0429 ← short‑chain dehydrogenase/reductase family oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947516,6940CA10.2% 218.2 / 7.9 73E100* (GAA→TAA) ‡PP_0429short‑chain dehydrogenase/reductase family oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (30/34);  new base A (5/3);  total (36/37)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATGGGCCAGGCGCCACCCCGGCCAGGCCGGCGACGTTGAGAATCTGCCCGCCACCCTGTACAGCCATCAGGTTACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGATACGCCGCTCGCCTTCTTCCAGCTCTAGCACGAAGCCTGAAAGGCGCAGCGGCTGGCTGAGGTCGCAGGCGCGGAACAGCACCTCGA  >  NC_002947/516548‑516837
                                                                                                                                                  |                                                                                                                                                 
ccatGGGCCAGGCGCCACCCCGGCCAGGCCGGCGACGTTGAGAATCTGCCCGCCACCCTGTACAGCCATCAGGTTACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCAt                                                                                                                                                 >  1:51545/1‑149 (MQ=255)
ccatGGGCCAGGCGCCACCCCGGCCAGGCCGGCGACGTTGAGAATCTGCCCGCCACCCTGTACAGCCATCAGGTTACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCAt                                                                                                                                                 >  7:51125/1‑149 (MQ=255)
 catGGGCCAGGCGCCACCCCGGCCAGGCCGGCGACGTTGAGAATCTGCCCACCACCCTGTACAGCCATCAGGTTACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCAt                                                                                                                                                 <  5:255006/148‑1 (MQ=255)
   tGGGCCAGGCGCCACCCCGGCCAGGCCGGCGACGTTGAGAATCTGCCCGCCACCCTGTACAGCCATCAGGTTACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTAATggg                                                                                                                                              >  7:261228/1‑149 (MQ=255)
     ggCCAGGCGCCACCCCGGCCAGGCCGGCGACGTTGAGAATCTGCCCGCCACCCTGTACAGCCATCAGGTTACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGc                                                                                                                                             <  4:35016/148‑1 (MQ=255)
        cAGGCGCCACCCCGGCCAGGCCGGCGACGTTGAGAATCTGCCCGCCACCCTGTACAGCCATCAGGTTACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATGAATGGGccgac                                                                                                                                         <  8:130416/149‑4 (MQ=255)
        cAGGCGCCACCCCGGCCAGGCCGGCGACGTTGAGAATCTGCCCGCCACCCTGTACAGCCATCAGGTTACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATGAATGGGccgac                                                                                                                                         <  2:37021/149‑4 (MQ=255)
              ccACCCCGGCCAGGCCGGCGACGTTGAGACTCTTCCCGCCACCCTGTACAGCCATCAGGTTCCCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGc                                                                                                                                 <  4:230824/149‑1 (MQ=255)
               cACCCCGGCCAGGCCGGCGACGTTGAGAATCTGCCCGCCACCCTGTACAGCCATCAGGTTACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGcc                                                                                                                                >  6:378062/1‑149 (MQ=255)
                aCCCCGGCCAGGCCGGCGACGTTGAGAATCTGCCCGCCACCCTGTACAGCCATCAGGTTACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTAATGGGc                                                                                                                                             >  7:143336/1‑137 (MQ=255)
                aCCCCGGCCAGGCCGGCGACGTTGAGAATCTGCCCGCCACCCTGTACAGCCATCAGGTTACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTAATGGGc                                                                                                                                             <  8:143336/137‑1 (MQ=255)
                 ccccGGCCAGGCCGGCGACGTTGAGAATCTGCCCGCCACCCTGTACAGCCATCAGGTTACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCa                                                                                                                              <  4:236562/149‑1 (MQ=255)
                 ccccGGCCAGGCCGGCGACGTTGAGAATCTGCCCGCCACCCTGTACAGCCATCAGGTTACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCa                                                                                                                              <  6:352948/149‑1 (MQ=255)
                 ccccGGCCAGGCCGGCGACGTTGAGAATCTGCCCGCCACCCTGTACAGCCATCAGGTTACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCa                                                                                                                              <  8:202580/149‑1 (MQ=255)
                  cccGGCCAGGCCGGCGACGTTGAGAATCTGCCCGCCACCCTGTACAGCCATCAGGTTACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCa                                                                                                                              >  2:231842/1‑148 (MQ=255)
                   ccGGCCAGGCCGGCGACGTTGAGAATCTGCCCGCCACCCTGTACAGCCATCAGGTTACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATa                                                                                                                            <  2:51545/149‑1 (MQ=255)
                      gCCAGGCCGGCGACGTTGAGAATCTGCCCGCCACCCTGTACAGCCATCAGGTTACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGt                                                                                                                         <  6:148855/149‑1 (MQ=255)
                                    gTTGAGAATCTGCCCGCCACCCTGTACAGCCATCAGGTTACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCcgcg                                                                                                            >  2:71052/1‑148 (MQ=255)
                                     ttGAGAATCTGCCCGCCACCCTGTACAGCCATCAGGTTACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCcgcg                                                                                                            >  5:188301/1‑147 (MQ=255)
                                       tagaATCTGGCCGCCACCCTGTAGTGCCTTGAGGTTACCGATTGCATGGCAAAGGCGGCTCAGCGCTAGCACGTTGATTTCCACGAGATCCTGTTCTTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCTCAg                                                                                                         <  8:261228/147‑1 (MQ=255)
                                        agaATCTGCCCGCCACCCTGTACAGCCATCAGGTTACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGtt                                                                                                       >  5:171831/1‑149 (MQ=255)
                                          aaTCTGCCCGCCACCCTGTACAGCCATCAGGTTACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGa                                                                                                     <  6:171831/149‑1 (MQ=255)
                                          aaTCTGCCCGCCACCCTGTACAGCCATCAGGTTACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGa                                                                                                     >  6:109032/1‑149 (MQ=255)
                                           aTCTGCCCGCCACCCTGTACAGCCATCAGGTTACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGAc                                                                                                    <  2:123127/149‑1 (MQ=255)
                                           aTCTGCCCGCCACCCTGTACAGCCATCAGGTTACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGAc                                                                                                    <  2:118018/149‑1 (MQ=255)
                                           aTCTGCCCGCCACCCTGTACAGCCATCAGGTTACCGATGGCATGGCACAGGCCGCACAGGCTCAGCACTTAGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGa                                                                                                     <  1:325248/148‑1 (MQ=255)
                                            tCTGCCCGCCACCCTGTACAGCCATCAGGTTACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGAc                                                                                                    >  3:110780/1‑148 (MQ=255)
                                                  cGCCACCCTGTACAGCCATCAGGTTACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCaa                                                                                              >  4:79213/1‑148 (MQ=255)
                                                        ccTGTACAGCCATCAGGTTACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGcc                                                                                                                                            <  8:16536/98‑1 (MQ=255)
                                                        ccTGTACAGCCATCAGGTTACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGcc                                                                                                                                            >  7:16536/1‑98 (MQ=255)
                                                            tACAGCCATCAGGTTACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGATAcgc                                                                                    <  8:284108/148‑1 (MQ=255)
                                                             aCAGCCATCAGGTTACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGATAcg                                                                                     <  2:142238/146‑1 (MQ=255)
                                                             aCAGCCATCAGGTTACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGATAcg                                                                                     >  1:142238/1‑146 (MQ=255)
                                                                gCCATCAGGTTACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGATACGCCGCTc                                                                               <  4:149122/149‑1 (MQ=255)
                                                                 ccATCAGGTTACCGATGGCATGGCACAGGCGGCTAAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGATACGCCGCTCg                                                                              >  4:75533/1‑149 (MQ=255)
                                                                   aTCAGGTTACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGATACGCCGCTCGcc                                                                            <  4:110780/149‑1 (MQ=255)
                                                                   aTCAGGTTACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCAATTTTTAATGTGGAATAGGCCCATAGGTGCGCAACCCCGCGCAGTTGACCAGCAAAAAGATACGCCGCTCGc                                                                             >  1:127155/1‑149 (MQ=255)
                                                                        gTTACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGATACGCCGCTCGCcttctt                                                                       >  5:229469/1‑149 (MQ=255)
                                                                         ttACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGATACGCCGCTCGCCttcttc                                                                      <  6:188301/149‑1 (MQ=255)
                                                                             cGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGAc                                                                                                    <  5:290710/115‑1 (MQ=255)
                                                                             cGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGAc                                                                                                    >  6:290710/1‑115 (MQ=255)
                                                                             cGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGAt                                                                                        <  5:72098/127‑1 (MQ=255)
                                                                             cGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGAt                                                                                        >  6:72098/1‑127 (MQ=255)
                                                                                    atggcaCAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTAATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGATACGCCGCTCGCCTTCTTCCAGCTCTAGCa                                                           >  7:250287/1‑149 (MQ=255)
                                                                                        caccGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGATACGCCGCTCGCCTTCTTCCAGCTCTAGCACGaa                                                       <  7:267979/149‑1 (MQ=255)
                                                                                         acaGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGATACGCCGCTCGCCTTCTTCCAGCTCTAGCACGAAg                                                      <  8:51125/149‑1 (MQ=255)
                                                                                          caGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGATACGCCGCTCGCCTTCTTCCAGCTCTAGCACGAAGc                                                     <  3:79213/149‑1 (MQ=255)
                                                                                           aGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGATACGCCGCTCGCCTTCTTCCAGCTCTAGCACGAAGcc                                                    >  1:200613/1‑149 (MQ=255)
                                                                                            ggcggcTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGATACGCCGCTCGCCTTCTTCCAGCTCTAGCACGAAGCCt                                                   <  7:146302/149‑1 (MQ=255)
                                                                                                gcTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGATACGCCGCTCGCCTTCTTCCAGCTCTAGCACGAAGCCTGaaa                                               <  1:231842/149‑1 (MQ=255)
                                                                                                  tCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGATACGCCGCTCGCCTTCTTCCAGCTCTAGCACGAAGCCTGAAAgg                                             >  7:319784/1‑149 (MQ=255)
                                                                                                  tCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGAGCAGTTGACCAGCAAATCGATACGCCGCTCGCCTTCTTCCAGCTCTAGCACGAAGCCTGAAAgg                                             >  5:226969/1‑149 (MQ=255)
                                                                                                     gggCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCAGCGCAGTTGACCAGCAAATCGATACGCCGCTCGCCTTCTTCCAGATCTAGAACCAAGcc                                                    >  5:90396/1‑139 (MQ=255)
                                                                                                          aGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGATACGCCGCTCGCCTTCTTCCAGCTCTAGCACGAAGCCTGAAAGGCGCAGCgg                                     <  5:378062/149‑1 (MQ=255)
                                                                                                           gCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGATACGCCGCTCGCCTTCTTCCAGCTCTAGCACGAAGCCTGAAAGGCGCAGCggc                                    >  2:183437/1‑149 (MQ=255)
                                                                                                                   aCTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGATACGCCGCTCGCCTTCTTCCAGCTCTAGCACGAAGCCTGAAAGGCGCAGCGGCTGGCTGAg                            <  8:25515/149‑1 (MQ=255)
                                                                                                                   aCTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGATACGCCGCTCGCCTTCTTCCAGCTCTAGCACGAAGCCTGAAAGGCGCAGCGGCTGGCTGAg                            <  1:71052/149‑1 (MQ=255)
                                                                                                                     ttCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGATACGCCGCTCGCCTTCTTCCAGCTCTAGCACGAAGCCTGAAAGgc                                            <  8:270797/131‑1 (MQ=255)
                                                                                                                         agcagATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGATACGCCGCTCGCCTTCTTCCAGCTCTAGCACGAAGCCTGAAAGGCGCAGCGGCTGGCTGAGGTCg                        >  2:47041/1‑147 (MQ=255)
                                                                                                                         agcagATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGATACGCCGCTCGCCTTCTTCCAGCTCTAGCACGAAGCCTGAAAGGCGCAGCGGCTGGCTGAGGTCGcc                      >  2:41302/1‑148 (MQ=255)
                                                                                                                              aTCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGATACGCCGCTCGCCTTCTTCCAGCTCTAGCACGAAGCCTGCAAGGCGCAGCGGCTGGCTGAGGTCGCAGgcgc                 >  4:89031/1‑149 (MQ=255)
                                                                                                                               tCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGATACGCCGCTCGCCTTCTTCCAGCTCTAGCACGAAGCCTGAAAGGCGCAGCGGCTGGCTGAGGTCGCAGGcgcg                <  1:183437/149‑1 (MQ=255)
                                                                                                                                ccTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGATACGCCGCTCGCCTTCTTCCAGCTCTAGCACGAAGCCTGAAAGGCGCAGCGGCTGGCTGAGGTCGCAGGcgcg                >  5:114959/1‑148 (MQ=255)
                                                                                                                                   gTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGATACGCCGCTCGCCTTCTTCCAGCTCTAGCACGAAGCCTGAACGGCGCAGCGGCTGGCTGAGGTCGCAGGCGCGGAAc            >  8:60170/1‑149 (MQ=255)
                                                                                                                                     ttcctggCCCATGAATGGGCCGGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGATACGCCGCTCGCCTTCTTCCAGCTCTAGCACGAAGCCTGAAAGGCGCAGCGGCTGGCTGAGGTCGCAGGCGCGGAACAg          >  6:334881/6‑149 (MQ=255)
                                                                                                                                      cgtcgGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGATACGCCGCTCGCCTTCTTCCAGCTCTAGCACGAAGCCTGAAAGGCGCAGCGGCTGGCTGAGGTCGCAGGCGCGGAACAg          <  2:127155/148‑1 (MQ=255)
                                                                                                                                        tcgGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGATACGCCGCTCGCCTTCTTCCAGCTCTAGCACGAAGCCTGAAAGGCGCAGCGGCTGGCTGAGGTCGCAGGCGCGGAACAGCAc       <  5:109032/149‑1 (MQ=255)
                                                                                                                                        ctggCCCATGAATGGGCCGACGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGATACGCCGCTCGCCTTCTTCCAGCTCTAGCACGAAGCCTGAAAGGCGCAGCGGCTGGCTGAGGTCGCAGGCGCGGAACAGCa        >  1:5432/3‑149 (MQ=255)
                                                                                                                                           gCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGATACGCCGCTCGCCTTCTTCCAGCTCTAGCACGAAGCCTGAAAGGCGCAGCGGCTGGCTGAGGTCGCATGCGCGGAACAGCACCTc    >  7:196591/1‑149 (MQ=255)
                                                                                                                                           gCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGATACGCCGCTCGCCTTCTTCCAGCTCTAGCACGAAGCCTGAAAGGCGCAGCGGCTGGCTGAGGTCGCAGGCGCGGAACAGCACCTc    <  6:360322/149‑1 (MQ=255)
                                                                                                                                            cccATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGATACGCCGCTCGCCTTCTTCCAGCTCTAGCACGAAGCCTGAAAGGCGCAGCGGCTGGCTGAGGTCGCAGGCGCGGAACAGCACCTc    >  6:235610/1‑148 (MQ=255)
                                                                                                                                             ccATTCATGGGCCAGGAACGGCCCATAGGTGCGCGAGCCCGCGCAGTTGACCAGCAAATCGATACGCCGCTCGCCTTCTTCCAGCTCTAGCACGAAGCCTGAAAGGCGCAGCGGCTGGCTGAGGTCGCAGGCGCGGAACAGCACCTCGa  >  1:188031/1‑149 (MQ=255)
                                                                                                                                             ccATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCTGCGCAGTTGACCAGCAAATCGATACGCCGCTCGCCTTCTTCCAGCTCTAGCACGAAGCCTGAAAGGCGCAGCGGCTGGCTGAGGTCGCAGGCGCGGAACAGCACCTCGa  >  4:253919/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CCATGGGCCAGGCGCCACCCCGGCCAGGCCGGCGACGTTGAGAATCTGCCCGCCACCCTGTACAGCCATCAGGTTACCGATGGCATGGCACAGGCGGCTCAGGGCCAGCACGTTGACTTCCAGCAGATCCTGTTCGTCGGCCCATTCATGGGCCAGGAACGGCCCATAGGTGCGCAAGCCCGCGCAGTTGACCAGCAAATCGATACGCCGCTCGCCTTCTTCCAGCTCTAGCACGAAGCCTGAAAGGCGCAGCGGCTGGCTGAGGTCGCAGGCGCGGAACAGCACCTCGA  >  NC_002947/516548‑516837

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: