Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,645,319 A→G 6.3% A224A (GCA→GCG PP_4954 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,645,3190AG 6.3% 235.9 / 3.5 79A224A (GCA→GCGPP_4954membrane protein
Reads supporting (aligned to +/- strand):  ref base A (35/39);  new base G (3/2);  total (38/41)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

GCAGAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCACTGACCGAGGCGTTTT  >  NC_002947/5645172‑5645467
                                                                                                                                                   |                                                                                                                                                     
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                                                                                                                          <  4:226748/134‑1 (MQ=255)
                                          gAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTgg                                                                                                                          >  3:226748/1‑134 (MQ=255)
                                             gCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACt                                                                                                        >  5:91198/1‑149 (MQ=255)
                                               tCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTc                                                                                                       >  1:260435/1‑148 (MQ=255)
                                                   ggCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAAc                                                                                                                     <  1:318095/130‑1 (MQ=255)
                                                   ggCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAAc                                                                                                                     >  2:318095/1‑130 (MQ=255)
                                                   ggCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGGAgg                                                                                                  >  3:116677/1‑149 (MQ=255)
                                                      cAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAg                                                                                                                                                      <  1:174604/95‑1 (MQ=255)
                                                      cAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAg                                                                                                                                                      >  2:174604/1‑95 (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                                                                                        <  4:78791/147‑1 (MQ=255)
                                                               tGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTa                                                                                        >  3:78791/1‑147 (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                                                                                                                            <  5:332020/94‑1 (MQ=255)
                                                                                tCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGt                                                                                                                            >  6:332020/1‑94 (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                            >  7:99996/2‑140 (MQ=255)
                                                                                                                                  ttggCTGCGTTCAGCGCGGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTAccc                            <  8:99996/139‑1 (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)
                                                                                                                                                   aGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCACTGACCGAGGCGtttt  >  2:17949/1‑149 (MQ=255)
                                                                                                                                                   aGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCACTGACCGAGGCGtttt  >  3:295006/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCAGAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCACTGACCGAGGCGTTTT  >  NC_002947/5645172‑5645467

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: