Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,061,390 C→A 5.9% L107M (CTG→ATG) ‡ PP_3579 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,061,3900CA 5.9% 242.4 / 2.7 76L107M (CTG→ATG) ‡PP_3579membrane protein
Reads supporting (aligned to +/- strand):  ref base C (37/34);  new base A (3/2);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGCCCGGAAAAGGTCCGCAAGGTCGCCGGCGTCGAAGACGAGAACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCGGGCAAAACCACGAACCCCCACTGCGCTTCGCTGATGGCGAGCG  >  NC_002947/4061244‑4061532
                                                                                                                                                  |                                                                                                                                              
ccagccCGGAAAAGGTCCGCAAGGTCGCCGGCGTCGAAGACGAGAACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTg                                                                                                                                              >  3:65893/1‑149 (MQ=255)
ccagccCGGAAAAGGTCCGCAAGGTCGCCGGCGTCGAAGACGAGAACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTg                                                                                                                                              >  2:250756/1‑149 (MQ=255)
   gccCGGAAAAGGTCCGCAAGGTCGCCGGCGTCGAAGACGAGAACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGcc                                                                                                                                           <  8:18413/149‑1 (MQ=255)
   gccAGGAAAAGGTCCGCAAGGTCGCCGGCTTCGAAGACGAGAACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTGTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGcc                                                                                                                                           <  2:339393/149‑1 (MQ=255)
            aGGTCCGCAAGGTCGCCGGCGTCGAAGACGAGAACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTAt                                                                                                                                  <  6:296927/149‑1 (MQ=255)
             ggTCCGCAAGGTCGCCGGCGTCGAAGACGAGAACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCAGGGCCGTGccttgca                                                                                                                                 <  4:346637/149‑7 (MQ=255)
                ccGCAAGGTCGCCGGCGTCGAAGACGAGAACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCAGGGCCGTGcc                                                                                                                                      >  4:106251/1‑140 (MQ=255)
                ccGCAAGGTCGCCGGCGTCGAAGACGAGAACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCAGGGCCGTGcc                                                                                                                                      <  3:106251/141‑2 (MQ=255)
                 cGCAAGGTCGCCGGCGTCGAAGACGAGAACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTa                                                                                                                             <  1:102147/149‑1 (MQ=255)
                 cGCAAGGTCGCCGGCGTCGAAGACGAGAACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTa                                                                                                                             >  4:271029/1‑149 (MQ=255)
                 cGCAAGGTCGCCGGCGTCGAAGACGAGAACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTa                                                                                                                             >  2:54390/1‑149 (MQ=255)
                  gCAAGGTCGCCGGCGTCGAAGACGAGAACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCCGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTa                                                                                                                             >  7:149353/1‑148 (MQ=255)
                  gCAAGGTCGCCGGCGTCGAAGACGAGAACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGACTCCCGCCCCGGCGTGCAAGGCAGGACCATGGCGGTGCCCTAtctcac                                                                                                                             >  1:84877/1‑146 (MQ=255)
                            cGGCGTCGAAGACGAGAACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGa                                                                                                                  >  1:123245/1‑149 (MQ=255)
                             ggCGTCGAAGACGAGAACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGAc                                                                                                                 <  3:278406/149‑1 (MQ=255)
                                           gaACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTgg                                                                                                                       >  1:6545/1‑129 (MQ=255)
                                           gaACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTgg                                                                                                                       <  2:6545/129‑1 (MQ=255)
                                           gaACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGTCCGTCGCCGGTTCCt                                                                                                   <  6:390494/149‑1 (MQ=255)
                                           gaACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCt                                                                                                   >  4:367277/1‑149 (MQ=255)
                                           gaACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCt                                                                                                   >  2:50474/1‑149 (MQ=255)
                                              cGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGc                                                                                                                                       <  3:5813/110‑1 (MQ=255)
                                              cGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGc                                                                                                                                       >  4:5813/1‑110 (MQ=255)
                                               gCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGgct                                                                                               >  5:135316/1‑149 (MQ=255)
                                                ccGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGgctg                                                                                              <  7:189525/149‑1 (MQ=255)
                                                 cGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGgctgc                                                                                             <  6:135316/149‑1 (MQ=255)
                                                    ccTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGa                                                                                          >  3:285889/1‑149 (MQ=255)
                                                    ccTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGa                                                                                          >  7:343907/1‑149 (MQ=255)
                                                        gTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATc                                                                                        <  1:50474/147‑1 (MQ=255)
                                                        gTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATc                                                                                        <  1:250756/147‑1 (MQ=255)
                                                        gTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCgg                                                                                      <  1:54390/149‑1 (MQ=255)
                                                             gctgACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGTCTGCa                                                                                 >  2:362506/1‑149 (MQ=255)
                                                              ctgACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGc                                                                                  <  4:65893/147‑1 (MQ=255)
                                                                    gtgtGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCAc                                                                          >  6:367668/1‑149 (MQ=255)
                                                                      gtgCATCGCCTCCATCGCCAGCCTTGCCGCGGTCACCCTGCAAATGGTCTCAAGCCGCGGCCTGAAAGTCACGGCCCTGGCCGTGCACTATCTCTACACTGGCATGACCTTAGCAGTTTCCTGGATGCTGATCGGCTGCATAttaaacc                                                                        >  2:268065/1‑144 (MQ=255)
                                                                      gtgCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCTGCTGAATATTCAcca                                                                        >  2:262138/1‑148 (MQ=255)
                                                                         cATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTg                                                                      <  4:181330/148‑1 (MQ=255)
                                                                               cgccATCGCCAGCCTTCCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTg                                                                                                                                        <  4:113702/76‑1 (MQ=255)
                                                                                gccATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTAcgc                                                              <  3:367277/149‑1 (MQ=255)
                                                                                  cATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTAcgccc                                                            >  3:306676/1‑149 (MQ=255)
                                                                                    tCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTAcgcccgc                                                          <  3:271029/149‑1 (MQ=255)
                                                                                      gCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCt                                                        >  7:345769/1‑149 (MQ=255)
                                                                                       ccAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACGCTGCACTACGCCCGCCTg                                                       >  5:209911/1‑149 (MQ=255)
                                                                                         aGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGtt                                                     >  3:184389/1‑149 (MQ=255)
                                                                                            cTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTa                                                  <  8:343907/149‑1 (MQ=255)
                                                                                             ttGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGCCCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTac                                                 >  5:302747/1‑149 (MQ=255)
                                                                                                ccGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACAcc                                              >  3:118778/1‑149 (MQ=255)
                                                                                                   cgGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCggg                                           <  4:118778/149‑1 (MQ=255)
                                                                                                          ccTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATc                                                                                        >  1:37143/1‑97 (MQ=255)
                                                                                                          ccTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATc                                                                                        <  2:37143/97‑1 (MQ=255)
                                                                                                          ccTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATAGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCGGGCAAAAc                                     >  8:180053/1‑148 (MQ=255)
                                                                                                          ccTGCAACTGGGCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCgg                                                                                      <  3:244753/99‑1 (MQ=255)
                                                                                                           cTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGtt                                                     >  7:15316/1‑131 (MQ=255)
                                                                                                           cTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGtt                                                     <  8:15316/131‑1 (MQ=255)
                                                                                                           cTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCGGGCAAAACCa                                   <  5:224492/149‑1 (MQ=255)
                                                                                                           cTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCGGGCAAAACCa                                   <  2:123245/149‑1 (MQ=255)
                                                                                                              cAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCGGGCAAAACCACGa                                <  1:268065/149‑1 (MQ=255)
                                                                                                               aaCTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTTCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCGGGCAAAACCACGa                                <  1:262138/148‑1 (MQ=255)
                                                                                                                   ggTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCCCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTacac                                               <  3:389813/129‑1 (MQ=255)
                                                                                                                   ggTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCCCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTacac                                               >  4:389813/1‑129 (MQ=255)
                                                                                                                     tctcCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTg                                                                                   <  3:13935/91‑1 (MQ=255)
                                                                                                                     tctcCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTg                                                                                   >  4:13935/1‑91 (MQ=255)
                                                                                                                         cagcagCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTATAAACTGTCCAGACCGTCGCCGGTTCCTCGCTGCTGATCTGCTGCATAATAAACCTGCACTAAGCCCGCCTGTTATACAAAGGGCAAAACAACGAAACACCACTgc                     >  5:414818/1‑149 (MQ=255)
                                                                                                                         cAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCGGGCAAAACCACGAACCCCCACTg                      <  1:362506/148‑1 (MQ=255)
                                                                                                                         cAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCATCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCGGGCAAAACCACGAACCCCCAc                        <  7:253014/146‑1 (MQ=255)
                                                                                                                            ccGCGGTCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCGGGCAAAACCACGAACCCCCACTGCGCt                  >  2:20227/1‑149 (MQ=255)
                                                                                                                               cgGTCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCGGGCAAAACCACGAACCCCCACTGCGCTTCg               <  1:20227/149‑1 (MQ=255)
                                                                                                                                   cTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCGGGCAAAACCACGAACCCCCACTGCGCTTCGCTGa           >  5:13260/1‑149 (MQ=255)
                                                                                                                                    tGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTCCACCGGGCAAAACCACGAACCCCCACTGCGCTTCGCTGAt          >  4:275654/1‑149 (MQ=255)
                                                                                                                                     gCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCGGGCAAAACCACGAACCCCCACTGCGCTTCGCTGATg         >  5:224928/1‑149 (MQ=255)
                                                                                                                                       aaGGCACGGCCCTGGCCGTGCACTATCTCTATACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCGGGCAAAACCACGAACCCCCACTGCGCTTCGCTGATGGc       >  5:414139/1‑149 (MQ=255)
                                                                                                                                       aaGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACGCTGCACTACGCCCGCCTGTTCTACACCGGGCAAAACCACGAACCCCCACTGCGCTTCGCTGATGGc       <  6:209911/149‑1 (MQ=255)
                                                                                                                                       aaGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCACCTGATCGGCTGCATATTCACACTCCACTACGCCCGCCTCTTCTACACCGGGCAAAACCACGAAACCCCACAGCGCTTCGCtaaagga       >  2:248623/1‑143 (MQ=255)
                                                                                                                                        gtgCACGGCCAGGGCCGTGC‑CTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCGGGCAAAACCACGAACCCCCACTGCGCTTCGCTGATGGc       >  1:157393/3‑147 (MQ=255)
                                                                                                                                         ggCACGGGCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCgg                                                                                      <  6:230853/68‑1 (MQ=255)
                                                                                                                                         ggCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCgg                                                                                      >  5:230853/1‑68 (MQ=255)
                                                                                                                                            aCGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCGGGCAACACCACGAACCCCCACTGCGCTTCGCTGATGGCGAgcg  >  6:387229/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
CCAGCCCGGAAAAGGTCCGCAAGGTCGCCGGCGTCGAAGACGAGAACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCGGGCAAAACCACGAACCCCCACTGCGCTTCGCTGATGGCGAGCG  >  NC_002947/4061244‑4061532

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: ‑