Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,762,947 A→T 7.3% I64F (ATC→TTC)  PP_3326 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,762,9470AT 7.3% 187.4 / 5.4 69I64F (ATC→TTC) PP_3326hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (34/30);  new base T (2/3);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.71e-01

TCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGG  >  NC_002947/3762801‑3763094
                                                                                                                                                    |                                                                                                                                                   
tCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATc                                                                                                                                                   <  1:144877/149‑1 (MQ=255)
    cGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGTCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCgg                                                                                                                                               >  1:114452/1‑149 (MQ=255)
    cGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCgg                                                                                                                                               <  3:366755/149‑1 (MQ=255)
     gTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGt                                                                                                                                              >  8:216081/1‑149 (MQ=255)
     gTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGt                                                                                                                                              >  2:104083/1‑149 (MQ=255)
      tCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGt                                                                                                                                              >  7:8823/1‑148 (MQ=255)
        aaCGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCAc                                                                                                                                           >  8:317348/1‑149 (MQ=255)
            cAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTc                                                                                                                                       <  8:194063/149‑1 (MQ=255)
               gcgcCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGATCATGTTCCCGGTCACgtcgaac                                                                                                                                    <  8:194223/149‑7 (MQ=255)
                cgcCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAgc                                                                                                                                   <  4:176200/149‑1 (MQ=255)
                  cccACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAgcgc                                                                                                                                 <  5:38572/149‑1 (MQ=255)
                   ccACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGcc                                                                                                                                >  6:69527/1‑149 (MQ=255)
                        gCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCg                                                                                                                           <  3:114736/149‑1 (MQ=255)
                          cTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGa                                                                                                                         <  8:8823/149‑1 (MQ=255)
                             tgGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGGTCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATc                                                                                                                       <  5:34426/148‑1 (MQ=255)
                              ggACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCAt                                                                                                                     <  7:216081/149‑1 (MQ=255)
                               gACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATg                                                                                                                    >  1:185049/1‑149 (MQ=255)
                                           aCGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAAc                                                                                                        >  7:280277/1‑149 (MQ=255)
                                           aCGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAAc                                                                                                        >  2:299392/1‑149 (MQ=255)
                                               ggCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGCTCGACCGTGACCGGGATCATGTTCCCGGTCa                                                                                                                                            <  6:312517/110‑1 (MQ=38)
                                               ggCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGCTCGACCGTGACCGGGATCATGTTCCCGGTCa                                                                                                                                            >  5:312517/1‑110 (MQ=38)
                                                gCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACt                                                                                                    >  5:20558/1‑148 (MQ=255)
                                                 cAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGcc                                                                                                                  <  3:197743/133‑1 (MQ=255)
                                                 cAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGcc                                                                                                                  >  4:197743/1‑133 (MQ=255)
                                                 cAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTa                                                                                                  >  4:292052/1‑149 (MQ=255)
                                                   atatCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACCTTGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTAc                                                                                                 <  1:104083/149‑1 (MQ=255)
                                                    tatCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTAccc                                                                                               >  7:268578/1‑149 (MQ=255)
                                                       cGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGc                                                                                            >  8:217266/1‑149 (MQ=255)
                                                         aGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGATCGCCTCGCGGATCATGCCGTGGTGGCACCACTTACCCCGCCt                                                                                          >  5:51578/1‑149 (MQ=255)
                                                                 tACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCAc                                                                                   <  7:317348/148‑1 (MQ=255)
                                                                   cGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTg                                                                                <  8:280277/149‑1 (MQ=255)
                                                                      ccACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCg                                                                             <  2:114452/149‑1 (MQ=255)
                                                                      ccACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCg                                                                             >  6:331465/1‑149 (MQ=255)
                                                                       cACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGt                                                                            >  3:264809/1‑149 (MQ=255)
                                                                                     tAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACg                                                              >  7:18358/1‑149 (MQ=255)
                                                                                       aCGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACgg                                                             <  6:19554/148‑1 (MQ=255)
                                                                                         gCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCAc                                                                                                     <  2:154198/106‑1 (MQ=255)
                                                                                         gCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCAc                                                                                                     >  1:154198/1‑106 (MQ=255)
                                                                                         gCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGtgctgc                                                                       <  7:320744/136‑1 (MQ=255)
                                                                                         gCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGtgctgc                                                                       >  8:320744/1‑136 (MQ=255)
                                                                                           cTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCtt                                                        <  7:324354/149‑1 (MQ=255)
                                                                                             aCAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTg                                                       <  3:350040/148‑1 (MQ=255)
                                                                                             aCAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAAAGGGCCTTg                                                       >  4:350040/1‑148 (MQ=255)
                                                                                               aaGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGc                                                                                        >  3:342109/1‑113 (MQ=255)
                                                                                               aaGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGc                                                                                        <  4:342109/113‑1 (MQ=255)
                                                                                               aaGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCCGGAACATGATACCGGTAACGGGCCGGCGCCGGGCGGGTCATGCCGTGCTGGCACCACTTAACCCCCCTGCGGCACATGCCGTGCTGCGGGGCAACCGGCCTTGGt                                                     >  6:50160/1‑148 (MQ=255)
                                                                                                 ggCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTca                                                   >  2:35134/1‑148 (MQ=255)
                                                                                                    ggAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAgc                                               >  3:269936/1‑149 (MQ=255)
                                                                                                     gAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAgcg                                              >  1:37519/1‑149 (MQ=255)
                                                                                                     gAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAgcg                                              >  2:50827/1‑149 (MQ=255)
                                                                                                          gATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCCCGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCg                                          >  7:38630/1‑148 (MQ=255)
                                                                                                              aaGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGTCAACGGGCCTTGGTCACAGCGCCCGGACaa                                     <  1:35134/149‑1 (MQ=255)
                                                                                                                gcgcgcCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAg                                    >  5:32372/1‑148 (MQ=255)
                                                                                                                    gcCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATg                               <  6:20558/149‑1 (MQ=255)
                                                                                                                    cgccaggCGCTCGACCGTGACCGGGATCATGTTCCCGGTCACGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATg                               >  4:45687/6‑149 (MQ=37)
                                                                                                                     ccTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATgg                              <  8:268578/149‑1 (MQ=255)
                                                                                                                      cTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATggg                             <  6:32372/149‑1 (MQ=255)
                                                                                                                        gggCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGAGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGt                            >  3:126177/1‑148 (MQ=255)
                                                                                                                            gTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTcaca                                                 >  7:175760/1‑123 (MQ=255)
                                                                                                                            gTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTcaca                                                 <  8:175760/123‑1 (MQ=255)
                                                                                                                              tCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCGGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAg                     >  5:196107/1‑149 (MQ=255)
                                                                                                                               cGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAg                     <  1:299392/148‑1 (MQ=255)
                                                                                                                                 aCGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGAt                  >  8:3871/1‑149 (MQ=255)
                                                                                                                                 aCGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGAt                  <  4:269936/149‑1 (MQ=255)
                                                                                                                                 aCGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGAt                  <  8:221480/149‑1 (MQ=255)
                                                                                                                                     tGACCGGGATCATGTTCCCGGTCACGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCa              <  3:45687/149‑1 (MQ=255)
                                                                                                                                           ggAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCAt         <  2:37519/149‑1 (MQ=255)
                                                                                                                                                 aTGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTc      >  7:180650/1‑147 (MQ=255)
                                                                                                                                                 aTGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCaa    <  5:331465/149‑1 (MQ=255)
                                                                                                                                                    aTCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAgg  <  4:264809/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGG  >  NC_002947/3762801‑3763094

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: