Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,460,731 G→C 6.2% L25V (CTG→GTG)  PP_2153 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,460,7310GC 6.2% 193.0 / 3.2 64L25V (CTG→GTG) PP_2153hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (32/28);  new base C (2/2);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.87e-01

CCGCCTTCGGAAATGACCACCGGCAGCGGCTCGCCGACTTCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGAAGTGCGATCCGATCTTCGATGCGGTAGTATT  >  NC_002947/2460583‑2460877
                                                                                                                                                    |                                                                                                                                                  
ccgccTTCGGAAATGACCACCGGCAGCGGCTCGCCGACTTCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGATGCGGTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAg                                                                                                                                                    <  7:144491/149‑1 (MQ=255)
  gccTTCGGAAATGACCACCGGCAGCGGCTCGCCGACTTCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGgt                                                                                                                                                  >  8:268391/1‑149 (MQ=255)
  gccTTCGGAAATGACCACCGGCAGCGGCTCGCCGACTTCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGgt                                                                                                                                                  <  3:159602/149‑1 (MQ=255)
     ttCGGAAATGACCACCGGCAGCGGCTCGCCGACTTCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGtt                                                                                                                                               >  5:16913/1‑149 (MQ=255)
          aaaTGACCACCGGCAGCGGCTCGCCGACTTCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGAc                                                                                                                                          <  6:16913/149‑1 (MQ=255)
          aaaTGACCACCGGCAGCGGCTCGCCGACTTCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGAc                                                                                                                                          <  7:268391/149‑1 (MQ=255)
          aaaTGACCACCGGCAGCGGCTCGCCAACTTCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACACCTGTTGCgcc                                                                                                                                          <  8:219415/149‑3 (MQ=255)
                 caccGGCAGCGGCTCGCCGACTTCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCATGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGtt                                                                                                                                               <  1:91546/137‑1 (MQ=255)
                 caccGGCAGCGGCTCGCCGACTTCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGtt                                                                                                                                               >  2:91546/1‑137 (MQ=255)
                   ccGGCAGCGGCTCGCCGACTTCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCg                                                                                                                                 <  2:249781/149‑1 (MQ=255)
                     ggCAGCGGCTCGCCGACTTCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGc                                                                                                                               <  2:214694/149‑1 (MQ=255)
                      gCAGCGGCTCGCCGACTTCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGc                                                                                                                               >  3:204482/1‑148 (MQ=255)
                      gCAGCGGCTCGCCGACTTCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGGCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGc                                                                                                                               >  5:212170/1‑148 (MQ=255)
                       cAGCGGCTCGCCGACTTCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCa                                                                                                                             <  6:120581/149‑1 (MQ=255)
                         gCGGCTCGCCGACTTCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAgg                                                                                                                           <  1:95908/149‑1 (MQ=255)
                                      ttCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCt                                                                                                                >  5:41303/1‑147 (MQ=255)
                                      ttCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCt                                                                                                                >  8:121692/1‑147 (MQ=255)
                                                  ggTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACa                                                                                                  >  8:250427/1‑149 (MQ=255)
                                                       gCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAgcggc                                                                                             >  6:188387/1‑149 (MQ=255)
                                                       gCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAgcggc                                                                                             >  6:236326/1‑149 (MQ=255)
                                                        cccAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAgcggc                                                                                             <  1:171857/148‑1 (MQ=255)
                                                        cccAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAgcggc                                                                                             >  1:205052/1‑148 (MQ=255)
                                                         ccAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAACTCGCTGAGCAGGTTGAACAgcggcg                                                                                            >  7:264613/1‑148 (MQ=255)
                                                            accaccGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAatc                                                                                        <  8:227051/149‑1 (MQ=255)
                                                                 cgcgCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCtt                                                                                   <  4:204482/149‑1 (MQ=255)
                                                                    gcTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGc                                                                                 >  5:89710/1‑148 (MQ=255)
                                                                          cAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACAt                                                                          >  8:316576/1‑149 (MQ=255)
                                                                           aaGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATc                                                                         <  1:217515/149‑1 (MQ=255)
                                                                                   gcgcTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGCCTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGCCGAATCATCTTGCAACACCTCTGGTTc                                                                   >  6:91058/1‑147 (MQ=255)
                                                                                    cgcTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAg                                                                <  7:250427/149‑1 (MQ=255)
                                                                                            cTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCgg                                                        >  4:154141/1‑149 (MQ=255)
                                                                                            cTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCg                                                         <  3:190968/148‑1 (MQ=255)
                                                                                                gCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCcgc                                                    >  5:227196/1‑149 (MQ=255)
                                                                                                 cccGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAgg                                                               >  8:87966/1‑137 (MQ=255)
                                                                                                 cccGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAgg                                                               <  7:87966/137‑1 (MQ=255)
                                                                                                 cccGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCcgcg                                                   >  2:367597/1‑149 (MQ=255)
                                                                                                  ccGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCgcgc                                                  >  7:120850/1‑149 (MQ=255)
                                                                                                        gAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAggg                                            >  7:52703/1‑149 (MQ=255)
                                                                                                                 gAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTAttt                                   <  8:52703/149‑1 (MQ=255)
                                                                                                                 gAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTAttt                                   <  8:120850/149‑1 (MQ=255)
                                                                                                                   gggTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGa                                 <  7:121692/149‑1 (MQ=255)
                                                                                                                    ggTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCTAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGaa                                <  6:41303/149‑1 (MQ=255)
                                                                                                                    ggTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTg                                  >  4:74451/1‑147 (MQ=255)
                                                                                                                    ggTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGaa                                <  1:367597/149‑1 (MQ=255)
                                                                                                                    ggTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGaa                                <  2:205052/149‑1 (MQ=255)
                                                                                                                        cGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGAAGTg                             <  3:74451/148‑1 (MQ=255)
                                                                                                                           tCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGAAGTGcgat                         <  5:236326/149‑1 (MQ=255)
                                                                                                                           tCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGAAGTGcgat                         >  1:202855/1‑149 (MQ=255)
                                                                                                                           tCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGAAGTGcgat                         >  7:220098/1‑149 (MQ=255)
                                                                                                                             cTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGAAGTGcgatcc                       >  2:54156/1‑149 (MQ=255)
                                                                                                                              ttGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGAAGTGcgatcc                       >  1:126282/1‑148 (MQ=255)
                                                                                                                              ttGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGAAGTGcgatcc                       >  5:220182/1‑148 (MQ=255)
                                                                                                                                    cTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGAAGTGCGATCCGATCTTc                <  7:316576/149‑1 (MQ=255)
                                                                                                                                        aCTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGAAGTGCGATCCGATCTTCGATg            <  5:188387/149‑1 (MQ=255)
                                                                                                                                         cTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGAAGTGCGATCCGATCTTCGATGc           >  6:68880/1‑149 (MQ=255)
                                                                                                                                          tGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGAAGTGCGATCCGATCTTCGATGCg          >  8:296592/1‑149 (MQ=255)
                                                                                                                                          tGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGAAGTGCGATCCGATCTTCGATGCg          >  1:172832/1‑149 (MQ=255)
                                                                                                                                          agtcgcAACACCTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGAAGTGCGATCCGATCTTCGATGCg          >  2:32040/4‑149 (MQ=255)
                                                                                                                                           gtcgcAACACCTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGGTAAGGGGGCTTATTTGAAGTGCGATCCGATCTTCGATGCgg         >  3:227680/3‑149 (MQ=255)
                                                                                                                                           gtcgcAACACCTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGAAGTGCGATCCGATCTTCGATGCgg         <  1:32040/147‑1 (MQ=255)
                                                                                                                                            gcgcAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCACCACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGAAGTGCGATCCGATCTTCGATGCg          <  2:320509/147‑1 (MQ=255)
                                                                                                                                            gcgcAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGAAGTGCGATCCGATCTTCGATGCgg         <  8:220098/148‑1 (MQ=255)
                                                                                                                                                aaCAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTAt                                     <  6:356531/116‑1 (MQ=255)
                                                                                                                                                aaCAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTAt                                     >  5:356531/1‑116 (MQ=255)
                                                                                                                                                  cAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGAAGTGCGATCCGATCTTCGATGCGGTAGTAtt  >  5:248214/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCGCCTTCGGAAATGACCACCGGCAGCGGCTCGCCGACTTCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGAAGTGCGATCCGATCTTCGATGCGGTAGTATT  >  NC_002947/2460583‑2460877

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