Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,591,619 T→A 5.4% L258M (TTG→ATG) ‡ PP_3169 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,591,6190TA 5.4% 463.3 / 12.9 149L258M (TTG→ATG) ‡PP_3169membrane protein
Reads supporting (aligned to +/- strand):  ref base T (65/76);  new base A (4/4);  total (69/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.02e-01

Only 100 of 152 total aligned reads displayed.

GGTCTGGCGGGTGCTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATGGTCGGT  >  NC_002947/3591471‑3591765
                                                                                                                                                    |                                                                                                                                                  
ggTCTGGCGGGTGCTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCtt                                                                                                                                                   <  2:324580/150‑1 (MQ=255)
 gTCTGGCGGGTGCTGGCCAAGACCCCGAATGGCGATTACTGCGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTg                                                                                                                                                  >  1:3047640/1‑150 (MQ=255)
    tGGCGGGTGCTGGTCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTgcgc                                                                                                                                               <  1:1263275/150‑1 (MQ=255)
      gCGGGTGCTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTg                                                                                                                                             <  2:1147179/150‑1 (MQ=255)
      gCGGGTGCTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTg                                                                                                                                             <  2:2643778/150‑1 (MQ=255)
       cGGGTGCTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTgg                                                                                                                                            <  1:1272441/150‑1 (MQ=255)
           tGCTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCa                                                                                                                                        <  2:586085/150‑1 (MQ=255)
               ggCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGa                                                                                                                                    >  2:2449708/1‑150 (MQ=255)
               agCCAAGACCCCGAATGGTCAGTACTCCGAACACATCAGCAGCGTGCTCCACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGAACCAGCGCAAGCGCTGGTGCAACgg                                                                                                                                    <  1:793732/149‑2 (MQ=255)
                gCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGAt                                                                                                                                   <  1:372033/150‑1 (MQ=255)
                gCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGAt                                                                                                                                   >  1:1603836/1‑150 (MQ=255)
                gCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGAt                                                                                                                                   <  2:772914/150‑1 (MQ=255)
                        cccGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCt                                                                                                                           >  1:2804619/1‑150 (MQ=255)
                        cccGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCt                                                                                                                           >  2:1878816/1‑150 (MQ=255)
                         ccGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTg                                                                                                                          >  1:737169/1‑150 (MQ=255)
                                gCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATg                                                                                                                   <  1:717736/150‑1 (MQ=255)
                                gCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATg                                                                                                                   <  2:2229697/150‑1 (MQ=255)
                                  gATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGAc                                                                                                                 <  2:2355039/150‑1 (MQ=255)
                                      actacGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGa                                                                                                             >  1:998719/1‑150 (MQ=255)
                                       ctacGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCTAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGAACAGAt                                                                                                            <  1:645757/150‑1 (MQ=255)
                                        tacGCAGACATCAGCGACGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGCACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGAtt                                                                                                           >  2:1062149/1‑150 (MQ=255)
                                        tacGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGGACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGAtt                                                                                                           >  1:1115390/1‑150 (MQ=255)
                                        tacGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTgg                                                                                                                                            >  1:2975317/1‑117 (MQ=255)
                                        tacGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTgg                                                                                                                                            <  2:2975317/117‑1 (MQ=255)
                                        tacGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCAAGCGCTGGt                                                                                                                                           <  1:2731775/118‑1 (MQ=255)
                                         acGAAGGCATCAGCAGCGTTCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTg                                                                                                          <  1:157008/150‑1 (MQ=255)
                                          cGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTgg                                                                                                         <  1:2449708/150‑1 (MQ=255)
                                           gAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGc                                                                                                        >  1:944374/1‑150 (MQ=255)
                                           gAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGc                                                                                                        <  2:890663/150‑1 (MQ=255)
                                              ggCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGAc                                                                                                     <  2:1853165/150‑1 (MQ=255)
                                                 aTCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGc                                                                                                        <  2:1716476/144‑1 (MQ=255)
                                                 aTCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCa                                                                                                  >  1:259949/1‑150 (MQ=255)
                                                 aTCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCa                                                                                                  >  2:1049195/1‑150 (MQ=255)
                                                  tCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGa                                                                                                                  <  2:1031425/133‑1 (MQ=255)
                                                  tCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGa                                                                                                                  >  1:1031425/1‑133 (MQ=255)
                                                    agcagcGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACtg                                                                                               <  1:1036235/150‑1 (MQ=255)
                                                     gcagcGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACtgg                                                                                              >  1:1078137/1‑150 (MQ=255)
                                                         cGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTggtggt                                                                                          <  1:494876/150‑1 (MQ=255)
                                                         cGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTggtggt                                                                                          >  2:2690162/1‑150 (MQ=255)
                                                            gCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCAc                                                                                       <  2:2156143/150‑1 (MQ=255)
                                                             cTCGATCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCAcc                                                                                      >  2:340977/1‑150 (MQ=255)
                                                              tCGACCAGCAACCACCCGAAATGCAATGTTGGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCg                                                                                     >  1:859187/1‑150 (MQ=255)
                                                              tCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGc                                                                                                        <  1:1811011/131‑1 (MQ=255)
                                                              tCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGc                                                                                                        >  2:1811011/1‑131 (MQ=255)
                                                                gACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGAc                                                                                   >  1:2533723/1‑150 (MQ=255)
                                                                gACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGAc                                                                                   >  1:1449960/1‑150 (MQ=255)
                                                                   cAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTAt                                                                                                                    >  2:1694673/1‑114 (MQ=255)
                                                                   cAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTAt                                                                                                                    <  1:1694673/114‑1 (MQ=255)
                                                                       aaCCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCa                                                                            >  1:450379/1‑150 (MQ=255)
                                                                         ccaccCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTggtggt                                                                                          <  1:1116271/134‑1 (MQ=255)
                                                                         ccaccCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATg                                                                          <  2:3047640/150‑1 (MQ=255)
                                                                          caccCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATgg                                                                         <  2:1081591/150‑1 (MQ=255)
                                                                          caccCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATgg                                                                         <  2:2533723/150‑1 (MQ=255)
                                                                            cccGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCAc                                                                                                 <  1:441416/124‑1 (MQ=255)
                                                                              cGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACtggtg                                                                                            >  1:550483/1‑127 (MQ=255)
                                                                                   tGCAACGCTTGTCGCGCAATCCGGAACTGGCGCAAGCAGTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGAGCATGGGGATCCCgg                                                                <  1:2204846/150‑1 (MQ=255)
                                                                                    gCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGc                                                                                                                         <  1:2043811/92‑1 (MQ=255)
                                                                                    gCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGc                                                                                                                         >  2:2043811/1‑92 (MQ=255)
                                                                                    gCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGa                                                                                                                  <  2:570454/99‑1 (MQ=255)
                                                                                      aaCGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGa                                                             >  2:2389792/1‑150 (MQ=255)
                                                                                        cGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATgg                                                                         <  1:1745928/136‑1 (MQ=255)
                                                                                         gCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCAc                                                                                                                                       >  2:1589155/1‑73 (MQ=255)
                                                                                          cTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTa                                                         <  1:2389792/150‑1 (MQ=255)
                                                                                            tGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTaca                                                       <  2:1539300/150‑1 (MQ=255)
                                                                                            tGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGATGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTaca                                                       <  2:1398974/150‑1 (MQ=255)
                                                                                             gTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACtggt                                                                                             <  2:2475605/111‑1 (MQ=255)
                                                                                             gTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACtggt                                                                                             >  1:2475605/1‑111 (MQ=255)
                                                                                             gTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTacac                                                      <  2:669481/150‑1 (MQ=255)
                                                                                              tCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACAcc                                                     >  1:2377058/1‑150 (MQ=255)
                                                                                                    aaCCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTgg                                                                                                         >  1:2037726/1‑92 (MQ=255)
                                                                                                    aaCCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTgg                                                                                                         <  2:2037726/92‑1 (MQ=255)
                                                                                                      ccTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTg                                                                                                          >  1:64554/1‑89 (MQ=255)
                                                                                                      ccTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGAATCCCGGGGAATTACACCTTTCGCtt                                             <  2:453354/150‑1 (MQ=255)
                                                                                                         ggAACTGGCGCAAGCACTGAACGATGACCCGTTGAACCAGCGCAAGCGCTgg                                                                                                                                            >  1:843636/1‑52 (MQ=25)
                                                                                                         ggAACTGGCGCAAGCACTGAACGATGACCCGTTGAACCAGCGCAAGCGCTgg                                                                                                                                            <  2:843636/52‑1 (MQ=25)
                                                                                                           aaCTGGCGCAAGCACTTAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGCCTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGTTCACCGACCTGCGCATGGGGAGCCCGGGGAATTACACCTTTCGCTTCAACa                                        <  2:1115390/150‑1 (MQ=255)
                                                                                                            aCTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACAt                                       >  2:159671/1‑150 (MQ=255)
                                                                                                                   cAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAg                                >  2:1137346/1‑150 (MQ=255)
                                                                                                                        aCTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGa                           <  2:1925232/150‑1 (MQ=255)
                                                                                                                        aCTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGa                           <  2:2804619/150‑1 (MQ=255)
                                                                                                                        aCTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGa                           >  2:2558836/1‑150 (MQ=255)
                                                                                                                         cTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAgcgcaac                          >  1:1978355/1‑147 (MQ=255)
                                                                                                                           gAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCtt                                             <  1:2889055/129‑1 (MQ=255)
                                                                                                                           gAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCtt                                             >  2:2889055/1‑129 (MQ=255)
                                                                                                                            aaCGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTc                                            >  1:225553/1‑129 (MQ=255)
                                                                                                                               gATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTaca                                                       <  1:1976964/115‑1 (MQ=255)
                                                                                                                                  gACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAggg                 <  2:958230/150‑1 (MQ=255)
                                                                                                                                    cccGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCa               >  1:2443430/1‑150 (MQ=255)
                                                                                                                                        gtgaaCCAGCGCAAGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTgg           >  2:873575/5‑150 (MQ=255)
                                                                                                                                         tGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCTCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGa          <  1:2708089/150‑1 (MQ=255)
                                                                                                                                           aaCCAGCGCAAGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTaca                                                       <  2:2482205/102‑1 (MQ=255)
                                                                                                                                            aCCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCTTGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATgg       <  2:450379/150‑1 (MQ=255)
                                                                                                                                             ccAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATGGt      <  2:1537469/150‑1 (MQ=255)
                                                                                                                                              cAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATGGTc     >  2:1431286/1‑150 (MQ=255)
                                                                                                                                              cAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATGGTc     >  1:346780/1‑150 (MQ=255)
                                                                                                                                               aGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATGGTCg    <  1:2558836/150‑1 (MQ=255)
                                                                                                                                               aGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTATCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATGGTCg    >  1:2636167/1‑150 (MQ=255)
                                                                                                                                                 cgcTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATGGTCGgt  <  1:1062149/150‑1 (MQ=255)
                                                                                                                                                    ttGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCggg                                                               <  1:706043/86‑1 (MQ=255)
                                                                                                                                                    ttGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCggg                                                               >  2:706043/1‑86 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GGTCTGGCGGGTGCTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATGGTCGGT  >  NC_002947/3591471‑3591765

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