Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,048,361 C→A 5.5% S139Y (TCC→TAC)  PP_2661 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,048,3610CA 5.5% 234.8 / 2.8 73S139Y (TCC→TAC) PP_2661hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (33/36);  new base A (2/2);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-01

TCCGCCGGTGGCTGATGCTTTACGGCATTCATGTCGATGGCGTGGTCAACAGCGGTCGGCACCAGCACATCCTGGCCCAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATG  >  NC_002947/3048214‑3048509
                                                                                                                                                   |                                                                                                                                                    
tCCGCCGGTGGCTGATGCTTTACGGCATTCATGTCGATGGCGTGGTCAACAGCGGTCGGCACCAGCACATCCTGGCCCAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTcc                                                                                                                                                     >  4:142808/1‑149 (MQ=255)
tCCGCCGGTGGCTGATGCTTTACGGCATTCATGTCGATGGCGTGGTCAACAGCGGTCGGCACCAGCACATCCTGGCCCAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTcc                                                                                                                                                     >  2:371892/1‑149 (MQ=255)
      ggTGGCTGATGCTTTACGGCATTGATGTCGTTGGCGTGGTCAACAGCGGTGGGCACGAGCTCTTCCTGGGCCAGCGCGGGCTTGAAAACGCTCCCTCCTTGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGgt                                                                                                                                               <  7:70395/149‑1 (MQ=255)
            tGATGCTTTACGGCATTCATGTCGATGGCGTGGTCAACAGCGGTCGGCACCAGCACATCCTGGCCCAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTgcgc                                                                                                                                         >  5:389153/1‑149 (MQ=255)
              aTGCTTTACGGCATTCATGTCGATGGCGTGGTCAACAGCGGTCGGCACCAGCACATCCTGGCCCAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCt                                                                                                                                       <  4:211615/149‑1 (MQ=255)
              aTGCTTTACGGCATTCATGTCGATGGCGTGGTCAACAGCGGTCGGCACCAGCACATCCTGGCCCAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTCTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGcc                                                                                                                                        <  1:204102/148‑1 (MQ=255)
              aTGCTTTACGGCATTCATGTCGATGGCGTGGTCAACAGCGGTCGGCACCAGCACATCCTGGCCCAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTCTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGcc                                                                                                                                        <  1:210604/148‑1 (MQ=255)
                  tttACGGCATTCATGTCGATGGCGTGGTCAACAGCGGTCGGCACCAGCACATCCTGGCCCAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTgg                                                                                                                                     >  1:134848/1‑147 (MQ=255)
                  tttACGGCATTCATGTCGATGGCGTGGTCAACAGCGGTCGGCACCAGCACATCCTGGCCCAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGACTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCATgg                                                                                                                                     >  2:416486/1‑147 (MQ=255)
                      cGGCATTCATGTCGATGGCGTGGTCAACAGCGGTCGGCACCAGCACATCCTGGCCCAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGgtt                                                                                                                               >  8:73244/1‑148 (MQ=255)
                               tGTCGATGGCGTGGTCAACAGCGGTCGGCACCAGCACATCCTGGCCCAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACg                                                                                                                      >  2:293384/1‑149 (MQ=255)
                                   gATGGCGTGGTCAACAGAGGTCGGCACCAGCACATCCTGGCCCAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTTCATTAGGATTTTGATATTGACGTGCTCTTAGATGATTCCGAAGGTGTGGGCCTGGAAGGTGTCGACCACGGtt                                                                                                                   <  4:403178/148‑1 (MQ=255)
                                        cGTGGTCAACAGCGGTCGGCACCAGCACATCCTGGCCCAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCcgcg                                                                                                             <  6:437010/149‑1 (MQ=255)
                                        cGTGGTCAACAGCGGTCGGCACCAGCACATCCTGGCCCAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCcgcg                                                                                                             >  5:437010/1‑149 (MQ=255)
                                         gTGGTCAACAGCGGTCGGCACCAGCACATCCTGGCCCAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTcc                                                                                                                >  1:322373/1‑145 (MQ=255)
                                         gTGGTCAACAGCGGTCGGCACCAGCACATCCTGGCCCAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTcc                                                                                                                <  2:322373/145‑1 (MQ=255)
                                           ggTCAACAGCGGTCGGGACCAGCACATCCTGGCCCAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCgtgg                                                                                                          <  1:83820/149‑1 (MQ=255)
                                             tCAACAGCGGTCGGCACCAGCACATCCTCGCCCAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCgtggtg                                                                                                        >  1:245955/1‑149 (MQ=255)
                                                aCAGCGGTCGGCACCAGCACATCCTGGCCCAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCgtggtggtg                                                                                                     >  2:12840/1‑149 (MQ=255)
                                                aCAGCGGTCGGCACCAGCACATCCTGGCCCAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCgtggtggtg                                                                                                     <  7:23015/149‑1 (MQ=255)
                                                aCAGCGGTCGGCACCAGCACATCCTGGCCCAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCgtggtggtg                                                                                                     <  1:12840/149‑1 (MQ=255)
                                                  aGCGGTCGGCACCAGCACATCCTGGCCCAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCgtggtggtggt                                                                                                   >  1:69178/1‑149 (MQ=255)
                                                  aGCGGTCGGCACCAGCACATCCTGGCCCAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCgtggtggtggt                                                                                                   <  4:75720/149‑1 (MQ=255)
                                                    cGGTCGGCACCAGCACATCCTGGCCCAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGgtgt                                                                                                 >  1:403692/1‑149 (MQ=255)
                                                    cGGTCGGCACCAGCACATCCTGGCCCAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGgtgt                                                                                                 <  1:371892/149‑1 (MQ=255)
                                                          gCACCAGCACATCCTGGCCCAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCg                                                                                            >  4:278812/1‑148 (MQ=255)
                                                           cACCAGCACATCCTGGCCCAGCGCGGTCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGaa                                                                                          <  4:86496/149‑1 (MQ=255)
                                                           cACCAGCACATCCTGGCCCAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTTTGccct                                                                                            >  2:230557/1‑146 (MQ=255)
                                                           cACCAGCACATCCTGGCCCAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCg                                                                                            >  5:269890/1‑147 (MQ=255)
                                                               aGCACATCCTGGCCCAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAgg                                                                                        >  1:78008/1‑147 (MQ=255)
                                                               aGCACATCCTGGCCCAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAgg                                                                                        <  2:78008/147‑1 (MQ=255)
                                                                gCACATCCTGGCCCAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATg                                                                                     <  8:131960/149‑1 (MQ=255)
                                                                     tCCTGGCCCAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCAc                                                                                <  7:73244/149‑1 (MQ=255)
                                                                       cTGGCCCAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTg                                                                              <  3:142808/149‑1 (MQ=255)
                                                                         ggCCCAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTACGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTggg                                                                            >  8:293531/1‑149 (MQ=255)
                                                                         ggCCCAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTACGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTggg                                                                            >  8:287717/1‑149 (MQ=255)
                                                                           cccAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCgtggt                                                                                                         <  6:113670/118‑1 (MQ=255)
                                                                           cccAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCgtggt                                                                                                         >  5:113670/1‑118 (MQ=255)
                                                                           cccAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGgcg                                                                          <  5:407031/149‑1 (MQ=255)
                                                                            ccAGCGCGGGCTGGAAAACGCTCCCTCCTAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGgcg                                                                          >  1:243630/1‑148 (MQ=255)
                                                                            ccAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGgcg                                                                          >  8:183376/1‑148 (MQ=255)
                                                                              aGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGgcg                                                                          >  2:87536/1‑146 (MQ=255)
                                                                               gcgcGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAg                                                                       >  8:299357/1‑148 (MQ=255)
                                                                                     gCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGa                                                                <  6:389153/149‑1 (MQ=255)
                                                                                     gCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGa                                                                <  1:421953/149‑1 (MQ=255)
                                                                                         gAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAg                                                              <  1:293384/147‑1 (MQ=255)
                                                                                         gAAAACGCTCCCTCCAAGCTGCCTTCGGCTATTGATATTCACCTGCACGAAGATGATTCCGAAGCTGAGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGAACTGTGCGCAGAACGTGAAGGa                                                            >  5:429498/1‑149 (MQ=255)
                                                                                          aaaaCGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGAc                                                           >  2:247105/1‑149 (MQ=255)
                                                                                                 tccctccAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCggtgg                                                    >  1:321339/1‑149 (MQ=255)
                                                                                                   cctccAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCggtgg                                                    <  2:403692/147‑1 (MQ=255)
                                                                                                   cctccAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGgtggtg                                                  <  6:269890/149‑1 (MQ=255)
                                                                                                    ctccAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGgtggtg                                                  >  2:341154/1‑148 (MQ=255)
                                                                                                       cAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGAtg                                              >  5:336412/1‑149 (MQ=255)
                                                                                                       cAAGCTGCCTTCGGCTTGTGATATTGACCTGCTCGTAGTTGATTCCGAAGGTGTGCGCCTGGAAGGTGTGGACCACGGTTTCCGCGTGGTGGTGGTGTTCCCGAAGGATGAGCACTGTGCGCAGTAGGTGAAGGACGCGGTGGTGGAtg                                              <  6:429498/149‑1 (MQ=255)
                                                                                                         aGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTACGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGAtgtg                                            <  7:287717/149‑1 (MQ=255)
                                                                                                         aGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTACGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGAtgtg                                            <  7:293531/149‑1 (MQ=255)
                                                                                                                      ttttGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAg                                                                    >  4:172975/1‑112 (MQ=255)
                                                                                                                      ttttGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAg                                                                    <  3:172975/112‑1 (MQ=255)
                                                                                                                        ttGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGctggcc                             >  6:2900/1‑149 (MQ=255)
                                                                                                                         tGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGGCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGctggcc                             <  1:230557/148‑1 (MQ=255)
                                                                                                                           atatTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGctggcctgg                          >  2:385566/1‑149 (MQ=255)
                                                                                                                              ttGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGgcag                       <  1:385566/149‑1 (MQ=255)
                                                                                                                              ttGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGgca                        <  7:299357/148‑1 (MQ=255)
                                                                                                                              ttGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGgca                        <  2:134848/148‑1 (MQ=255)
                                                                                                                                   cTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAAcc                  >  2:340162/1‑149 (MQ=255)
                                                                                                                                     gCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGctggcct                            <  3:372371/137‑1 (MQ=255)
                                                                                                                                     gCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGctggcct                            >  4:372371/1‑137 (MQ=255)
                                                                                                                                      cACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCa               <  7:183376/149‑1 (MQ=255)
                                                                                                                                        cGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCAcc             >  8:132841/1‑149 (MQ=255)
                                                                                                                                            gatgatTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTa         <  6:149979/149‑1 (MQ=255)
                                                                                                                                             atgatTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTAt        <  2:245955/149‑1 (MQ=255)
                                                                                                                                              tgatTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTAt        <  2:69178/148‑1 (MQ=255)
                                                                                                                                                 ttCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGaaa    >  4:127939/1‑149 (MQ=255)
                                                                                                                                                 ttCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGaaa    <  5:357269/149‑1 (MQ=255)
                                                                                                                                                  tCCGAAGGTGTGCGCCTGGAAGGTGTAGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAAt   >  6:31628/1‑149 (MQ=255)
                                                                                                                                                   ccGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGACGTGTGCCCGAAGGATGAGCAATCGGAGCAGAAGGTGAATGACGCCGTGGTGGATGTGAAGCCACGGATGGCCTGCAAGCAACCGCACCGTTATGAAATg  >  3:19102/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCCGCCGGTGGCTGATGCTTTACGGCATTCATGTCGATGGCGTGGTCAACAGCGGTCGGCACCAGCACATCCTGGCCCAGCGCGGGCTGGAAAACGCTCCCTCCAAGCTGCCTTCGGCTTTTGATATTGACCTGCACGTAGATGATTCCGAAGGTGTGCGCCTGGAAGGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATG  >  NC_002947/3048214‑3048509

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑