Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,493,464 T→C 5.1% intergenic (‑30/+22) PP_5715 ← / ← PP_5716 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,493,4640TC 5.1% 643.4 / 14.7 194intergenic (‑30/+22)PP_5715/PP_5716hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (81/103);  new base C (5/5);  total (86/108)
Fisher's exact test for biased strand distribution p-value = 7.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-01

Only 100 of 205 total aligned reads displayed.

GCGGTCAGGTAGGCTGCCAGCTGAACATCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGAAATTCAACACCGCACCACAGGCAAAATTGGAGGTGC  >  NC_002947/5493316‑5493613
                                                                                                                                                    |                                                                                                                                                     
gCGGTCAGGTAGGCTGCCAGCTGAACATCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCgtg                                                                                                                                                      >  2:1380061/1‑150 (MQ=255)
 cGGTCAGGTAGGCTGCCAGCTGAACATCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCgtgt                                                                                                                                                     <  2:1911591/150‑1 (MQ=255)
 cGGTCAGGTAGGCTGCCAGCTGAACATCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCgtgt                                                                                                                                                     <  2:266327/150‑1 (MQ=255)
 cGGTCAGGTAGGCTGCCAGCTGAACATCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCgtgt                                                                                                                                                     >  2:3148117/1‑150 (MQ=255)
       ggTAGGCTGCCAGCTGAACATCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAg                                                                                                                                               <  2:564604/150‑1 (MQ=255)
       ggTAGGCTGCCAGCTGAACATCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAg                                                                                                                                               <  2:1489469/150‑1 (MQ=255)
       ggTAGGCTGCCAGCTGAACATCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAg                                                                                                                                               <  2:4052468/150‑1 (MQ=255)
       ggTAGGCTGCCAGCTGAACATCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAg                                                                                                                                               <  1:1766035/150‑1 (MQ=255)
       ggTAGGCTGCCAGCTGAACATCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATCGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAg                                                                                                                                               <  2:3909694/150‑1 (MQ=255)
       ggTAGGCTGCCAGCTGAACATCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAACACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAg                                                                                                                                               <  2:2162065/150‑1 (MQ=255)
        gTAGGCTGCCAGCTGAACATCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAgg                                                                                                                                              <  1:2794721/150‑1 (MQ=255)
         tAGGCTGCCAGCTGAACATCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGc                                                                                                                                             >  1:191225/1‑150 (MQ=255)
             cTGCCAGCTGAACATCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAAc                                                                                                                                         <  1:1951205/150‑1 (MQ=255)
              tGCCAGCTGAACATCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCTGTCGACACGGCGTCGACGGGcctgga                                                                                                                                        <  1:3868515/150‑5 (MQ=255)
              tGCCAGCTGAACATCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACt                                                                                                                                        >  1:839867/1‑150 (MQ=255)
               gCCAGCTGAACATCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACtt                                                                                                                                       >  1:2415408/1‑150 (MQ=255)
                 cAGCTGAACATCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCa                                                                                                                                     <  2:4101427/150‑1 (MQ=255)
                  aGCTGAACATCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACACAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCAt                                                                                                                                    <  2:1080245/150‑1 (MQ=255)
                   gCTGAACATCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATg                                                                                                                                   <  1:1975715/150‑1 (MQ=255)
                   gCTGAACATCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATg                                                                                                                                   <  1:3768924/150‑1 (MQ=255)
                     tGAACATCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGt                                                                                                                                 <  2:2586900/150‑1 (MQ=255)
                      gAACATCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTc                                                                                                                                >  2:153475/1‑150 (MQ=255)
                       aaCATCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCa                                                                                                                               <  2:1993034/150‑1 (MQ=255)
                       aaCATCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCa                                                                                                                               <  1:2973745/150‑1 (MQ=255)
                       aaCATCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCa                                                                                                                               <  1:703997/150‑1 (MQ=255)
                         cATCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCaga                                                                                                                             <  1:773836/150‑1 (MQ=255)
                           tCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGt                                                                                                                           >  2:3054211/1‑150 (MQ=255)
                             cAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTc                                                                                                                         <  2:1764013/150‑1 (MQ=255)
                              aaCAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTcc                                                                                                                        <  1:2787346/150‑1 (MQ=255)
                                     cTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCca                                                                                                                 >  2:2920190/1‑150 (MQ=255)
                                     cTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCca                                                                                                                 <  1:2246899/150‑1 (MQ=255)
                                      ttCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCcag                                                                                                                >  2:1345647/1‑150 (MQ=255)
                                      ttCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCcag                                                                                                                >  1:3448200/1‑150 (MQ=255)
                                      ttCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCcag                                                                                                                >  2:2644241/1‑150 (MQ=255)
                                      ttCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCcag                                                                                                                >  2:3112995/1‑150 (MQ=255)
                                      ttCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCcag                                                                                                                >  1:2834745/1‑150 (MQ=255)
                                         acgagcagAAGCACTGACATTGTTGTCGATAAATACCAGGCTATTGCCGTTCAACCACAAAGAGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACATCATGGTCAGAGTTCCACTGCcagcag                                                                                                             >  1:1623307/8‑150 (MQ=255)
                                            tCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAggtc                                                                                                          <  2:3559470/150‑3 (MQ=255)
                                            tCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGAc                                                                                                          <  1:3072968/150‑1 (MQ=255)
                                            tCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGAc                                                                                                          >  2:3472520/1‑150 (MQ=255)
                                               tGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTg                                                                                                       <  2:4281293/150‑1 (MQ=255)
                                               tGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTg                                                                                                       <  2:2012503/150‑1 (MQ=255)
                                               tGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTg                                                                                                       <  1:3375616/150‑1 (MQ=255)
                                                  aGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAg                                                                                                    <  2:1459497/150‑1 (MQ=255)
                                                   gCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAgg                                                                                                   <  1:2873499/150‑1 (MQ=255)
                                                       tGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCAt                                                                                               <  1:576577/150‑1 (MQ=255)
                                                       tGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCAt                                                                                               >  1:2641087/1‑150 (MQ=255)
                                                       tGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCAt                                                                                               >  2:1012645/1‑150 (MQ=255)
                                                        gACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGTTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCAtt                                                                                              >  1:3807238/1‑150 (MQ=255)
                                                        gACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCAtt                                                                                              >  2:4352242/1‑150 (MQ=255)
                                                         aCATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTg                                                                                             <  2:1306197/150‑1 (MQ=255)
                                                            ttgttgTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTgtcg                                                                                          >  2:2689932/1‑150 (MQ=255)
                                                            ttgttgTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTgtcg                                                                                          >  2:841499/1‑150 (MQ=255)
                                                              gttgtCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTgtcgtc                                                                                        <  2:4178210/150‑1 (MQ=255)
                                                                tgtCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCaa                                                                                      <  2:4400671/150‑1 (MQ=255)
                                                                     aTCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTg                                                                                 >  1:3564295/1‑150 (MQ=255)
                                                                     aTCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTg                                                                                 <  2:387149/150‑1 (MQ=255)
                                                                     aTCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTg                                                                                 <  2:2212949/150‑1 (MQ=255)
                                                                        aaTACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGtt                                                                              >  2:477901/1‑150 (MQ=255)
                                                                          tACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTc                                                                            >  1:4179863/1‑150 (MQ=255)
                                                                          tACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCATCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGACTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTc                                                                            <  1:2934018/150‑1 (MQ=255)
                                                                           aCCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCt                                                                           <  2:647885/150‑1 (MQ=255)
                                                                                       ccGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGa                                                               >  1:2662643/1‑150 (MQ=255)
                                                                                              aCCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGt                                                        >  2:2304061/1‑150 (MQ=255)
                                                                                               ccACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTc                                                       <  2:3433739/150‑1 (MQ=255)
                                                                                                cacaAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCg                                                      >  1:4429663/1‑150 (MQ=255)
                                                                                                 acaAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGt                                                     >  2:2067837/1‑150 (MQ=255)
                                                                                                 acaAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGt                                                     >  2:2544672/1‑150 (MQ=255)
                                                                                                  caAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTa                                                    >  1:4325805/1‑150 (MQ=255)
                                                                                                  caAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTa                                                    >  2:2754058/1‑150 (MQ=255)
                                                                                                    aaGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATc                                                  <  2:2604743/150‑1 (MQ=255)
                                                                                                        gtAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATcggcg                                              <  1:4053827/150‑1 (MQ=255)
                                                                                                         tAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGt                                                        >  1:4514744/1‑139 (MQ=255)
                                                                                                          aaTTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCggcggc                                            <  1:1507157/150‑1 (MQ=255)
                                                                                                           aTTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCt                                           >  1:291325/1‑150 (MQ=255)
                                                                                                           aTTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCt                                           >  1:3981687/1‑150 (MQ=255)
                                                                                                            ttGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTa                                          <  1:4099596/150‑1 (MQ=255)
                                                                                                            ttGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTa                                          <  1:1053273/150‑1 (MQ=255)
                                                                                                                  tGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGaaa                                    <  1:3618411/150‑1 (MQ=255)
                                                                                                                  tGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGaaa                                    >  1:392556/1‑150 (MQ=255)
                                                                                                                  tGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGaaa                                    <  2:3344287/150‑1 (MQ=255)
                                                                                                                    cATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGAAAtt                                  <  1:1498090/150‑1 (MQ=255)
                                                                                                                         gTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGAAATTCAaca                             <  1:2644241/150‑1 (MQ=255)
                                                                                                                         gTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGAAATTCAaca                             >  1:1979318/1‑150 (MQ=255)
                                                                                                                          tATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGAAATTCAacac                            <  2:2385363/150‑1 (MQ=255)
                                                                                                                            tGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGAAATTCACCACCg                          <  1:1865620/150‑1 (MQ=255)
                                                                                                                              cccTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGAAATTCAACACCGca                        <  1:1321234/150‑1 (MQ=255)
                                                                                                                               ccTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGAAATTCAACACCGcac                       >  1:2331863/1‑150 (MQ=255)
                                                                                                                                cTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGAAATTCAACACCGcacc                      >  2:2457607/1‑150 (MQ=255)
                                                                                                                                cTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGAAATTCAACACCGcacc                      >  1:640085/1‑150 (MQ=255)
                                                                                                                                     ggCCTGTCGACACGGCGTCGACGGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGgc                                                                   >  1:1996007/1‑100 (MQ=38)
                                                                                                                                      gCCTGTCGACACGGCGTCGACGGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGAAATTCAACACCGCACCACAGGc                >  2:1432542/1‑150 (MQ=255)
                                                                                                                                        ctgTCGACACGGCGTCGACGGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGAAATTCAACACCGCACCACAGGCaa              >  1:2152378/3‑150 (MQ=255)
                                                                                                                                          gTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGAAATTCAACACCGCACCACAGGCaaaa            >  2:3782184/1‑150 (MQ=255)
                                                                                                                                          gTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGAAATTCAACACCGCACCACAGGCaaaa            >  1:2733197/1‑150 (MQ=255)
                                                                                                                                          gTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGAAATTCAACACCGCACCACAGGCaaaa            >  2:1621876/1‑150 (MQ=255)
                                                                                                                                               cGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTTCAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGAAATTCAACACCGCACCACAGGCAAAATTGGa       <  1:1227772/150‑1 (MQ=255)
                                                                                                                                               cGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGAAATTCAACACCGCACCACAGGCAAAATTGGa       <  2:739173/150‑1 (MQ=255)
                                                                                                                                                   gtgtCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGAAATTCAACACCGCACCACAGGCAAAATTGGAGGTg   <  2:3003668/150‑1 (MQ=255)
                                                                                                                                                   gtgtCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGAAATTCAACACCGCACCACAGGCAAAATTGGAGGTg   <  2:761481/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCGGTCAGGTAGGCTGCCAGCTGAACATCCAACAGTGCTTCTGCTCGTGAAGCACTGACATTGTTGTCGATCAATACCAGGCTATTGCCGTTCAACCACAAAGTGTAATTGGCTTGCATCAGTATGCCCTCCAGGCCCGTCGACGCCGTGTCGACAGGCCAACTTCATGGTCAGAGTTCCAGTGCCAGCAGGACCTGCAGGGCATTGTCGTCAAGTTTGGTTTCTACAATGGCGCGAACTTTGTCGTATCGGCGGCTATCGAAATTCAACACCGCACCACAGGCAAAATTGGAGGTGC  >  NC_002947/5493316‑5493613

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