Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,645,315 A→G 6.6% N223S (AAC→AGC)  PP_4954 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,645,3150AG 6.6% 224.2 / 4.1 76N223S (AAC→AGC) PP_4954membrane protein
Reads supporting (aligned to +/- strand):  ref base A (34/37);  new base G (3/2);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

GCAGAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCACTGACCGAG  >  NC_002947/5645172‑5645460
                                                                                                                                               |                                                                                                                                                  
gcagAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAg                                                                                                                                               >  6:235824/1‑149 (MQ=255)
 cagAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGc                                                                                                                                              >  1:15698/1‑149 (MQ=255)
  agagCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGc                                                                                                                                              >  3:89293/1‑148 (MQ=255)
              gACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGatcatc                                                                                                                                <  1:368755/149‑1 (MQ=255)
                        cGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCggg                                                                                                                      >  4:24917/1‑149 (MQ=255)
                          aGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCTAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTg                                                                                                                    <  2:231393/149‑1 (MQ=255)
                               ggCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGCTCATCCTGCACCGGGTGGGCaa                                                                                                               >  7:228837/1‑149 (MQ=255)
                               ggCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGc                                                                                                                 <  1:42501/147‑1 (MQ=255)
                                gCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTgg                                                                                                                   >  8:297353/1‑144 (MQ=255)
                                gCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTgg                                                                                                                   <  7:297353/144‑1 (MQ=255)
                                  aGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAg                                                                                                            <  3:24917/149‑1 (MQ=255)
                                       ccTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGccaa                                                                                                        <  8:182973/148‑1 (MQ=255)
                                       ccTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGccaa                                                                                                        >  7:182973/1‑148 (MQ=255)
                                          gAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTgg                                                                                                                   >  3:226748/1‑134 (MQ=255)
                                          gAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTgg                                                                                                                   <  4:226748/134‑1 (MQ=255)
                                             gCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACt                                                                                                 >  5:91198/1‑149 (MQ=255)
                                               tCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTc                                                                                                >  1:260435/1‑148 (MQ=255)
                                                   ggCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAAc                                                                                                              >  2:318095/1‑130 (MQ=255)
                                                   ggCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAAc                                                                                                              <  1:318095/130‑1 (MQ=255)
                                                   ggCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGGAgg                                                                                           >  3:116677/1‑149 (MQ=255)
                                                      cAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAg                                                                                                                                               >  2:174604/1‑95 (MQ=255)
                                                      cAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAg                                                                                                                                               <  1:174604/95‑1 (MQ=255)
                                                      cAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGGAGGTaa                                                                                        <  4:116677/149‑1 (MQ=255)
                                                      cAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGt                                                                                          >  3:107633/1‑147 (MQ=255)
                                                      cAGCCTGTTTGAACCGGGCCGCAAGCTCGACGTCGAGGCACTGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAATGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTaa                                                                                        <  5:256106/149‑1 (MQ=255)
                                                        gCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGc                                                                                      <  6:91198/149‑1 (MQ=255)
                                                         ccTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGcc                                                                                     <  4:169908/149‑1 (MQ=255)
                                                         ccTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGcc                                                                                     <  2:260435/149‑1 (MQ=255)
                                                             tttGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTa                                                                                 >  3:216163/1‑149 (MQ=255)
                                                              ttGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTAc                                                                                >  7:198562/1‑149 (MQ=255)
                                                               tGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTa                                                                                 >  3:78791/1‑147 (MQ=255)
                                                               tGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTa                                                                                 <  4:78791/147‑1 (MQ=255)
                                                               tGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATAATCCTGCACAGGGTGCGCAACAGCCCACCAAACTCGAAGGTAAGCCGCTAc                                                                                >  2:363851/1‑148 (MQ=255)
                                                                gAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTAcc                                                                               <  5:235824/148‑1 (MQ=255)
                                                                    cGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAg                                                                          >  2:232227/1‑149 (MQ=255)
                                                                    cGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAg                                                                          <  4:216163/149‑1 (MQ=255)
                                                                         cGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGt                                                                      <  1:379564/148‑1 (MQ=255)
                                                                                tCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGt                                                                                                                     >  6:332020/1‑94 (MQ=255)
                                                                                tCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGt                                                                                                                     <  5:332020/94‑1 (MQ=255)
                                                                                      tCGAGGCACGGCGCCTGGAACTGGCCCAGCATAACGGCAAGGTGGTGCCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTTGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAg                                                        <  8:261315/149‑1 (MQ=255)
                                                                                       cGAGGCACTGCGCCTGTAACTTGCCCAGCAGAACGGCAAGGTGGTGTCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAg                                                        <  8:266410/148‑1 (MQ=255)
                                                                                        gAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAg                                                                                                            >  5:167623/1‑95 (MQ=255)
                                                                                        gAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAgg                                                       <  3:167595/148‑1 (MQ=255)
                                                                                        gAGGCACGGCGACTGGAACTGGCCCAGCAGAACTGCAAGGTGGTGGACGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAg                                                                                                            <  6:167623/95‑1 (MQ=255)
                                                                                         aGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGat                                                     <  2:15698/149‑1 (MQ=255)
                                                                                             aCGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATcc                                                 <  1:38201/149‑1 (MQ=255)
                                                                                             aCGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATc                                                  >  3:93842/1‑148 (MQ=255)
                                                                                                  gcCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAg                                            >  3:65129/1‑149 (MQ=255)
                                                                                                      ggAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAg                                                                          <  1:220502/115‑1 (MQ=255)
                                                                                                       gAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCgt                                         <  4:93842/147‑1 (MQ=255)
                                                                                                       gAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCgt                                         <  3:137150/147‑1 (MQ=255)
                                                                                                             gCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAATGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTg                                  <  1:363851/148‑1 (MQ=255)
                                                                                                               ccAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGc                                           <  6:7573/137‑1 (MQ=255)
                                                                                                               ccAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGc                                           >  5:7573/1‑137 (MQ=255)
                                                                                                               ccAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCt                               >  2:312623/1‑149 (MQ=255)
                                                                                                                 agcagAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCg                             >  2:333596/1‑149 (MQ=255)
                                                                                                                   cagAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCa                           <  4:107633/149‑1 (MQ=255)
                                                                                                                    agAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCAc                          >  5:284525/1‑149 (MQ=255)
                                                                                                                         ggCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTAccc                     <  8:228837/149‑1 (MQ=255)
                                                                                                                           cAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGc                                           <  4:155740/125‑1 (MQ=255)
                                                                                                                           cAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGc                                           >  3:155740/1‑125 (MQ=255)
                                                                                                                           cAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTAcccc                    >  4:253779/1‑148 (MQ=255)
                                                                                                                               gtggtgGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACaa               >  5:256054/1‑149 (MQ=255)
                                                                                                                                tggtggCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAAc              <  8:198562/149‑1 (MQ=255)
                                                                                                                                ccttggCTGCGTTCAGCGCGGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAAc              >  1:39673/4‑149 (MQ=255)
                                                                                                                                  ttggCTGCGTTCAGCGCGGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTAccc                     <  8:99996/139‑1 (MQ=255)
                                                                                                                                  ttggCTGCGTTCAGCGCGGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTAccc                     >  7:99996/2‑140 (MQ=255)
                                                                                                                                  ttggCTGCGTTCAGCGCGGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAAc              >  7:54876/2‑147 (MQ=255)
                                                                                                                                  ttggCTGCGTTCAGCGCGGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTCTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAAc              <  8:54876/146‑1 (MQ=255)
                                                                                                                                  gtggCCGCGCTGAACGCAGCCAAGGAGATGCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACg             >  1:28958/1‑149 (MQ=255)
                                                                                                                                  gtggCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTAc                       >  1:81777/1‑138 (MQ=255)
                                                                                                                                  gtggCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTAc                       <  2:81777/138‑1 (MQ=255)
                                                                                                                                       cgcgCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCa                                                <  5:104679/108‑1 (MQ=255)
                                                                                                                                       cgcgCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCa                                                >  6:104679/1‑108 (MQ=255)
                                                                                                                                       cgcgCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCACTGa       >  8:253375/1‑149 (MQ=255)
                                                                                                                                           cTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCACTGACCGa   <  1:333596/149‑1 (MQ=255)
                                                                                                                                              aaCGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAg                                                        >  2:141246/1‑93 (MQ=255)
                                                                                                                                              aaCGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAg                                                        <  1:141246/93‑1 (MQ=255)
                                                                                                                                              aaCGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCACTGACCGAg  <  2:28958/147‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GCAGAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCACTGACCGAG  >  NC_002947/5645172‑5645460

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: