Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 337,839 C→T 20.9% intergenic (‑86/+43) PP_0279 ← / ← PP_0280 hypothetical protein/amino acid ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947337,8390CT20.9% 211.3 / 51.9 107intergenic (‑86/+43)PP_0279/PP_0280hypothetical protein/amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (50/33);  new base T (11/11);  total (63/44)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.04e-01

Only 100 of 112 total aligned reads displayed.

TCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAA  >  NC_002947/337690‑337986
                                                                                                                                                      |                                                                                                                                                   
tCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTac                                                                                                                                                     >  2:1827021/1‑150 (MQ=255)
tCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTGAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTac                                                                                                                                                     >  2:4233932/1‑150 (MQ=255)
 cAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTaca                                                                                                                                                    <  1:5312508/150‑1 (MQ=255)
    gAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCt                                                                                                                                                 <  2:2058843/150‑1 (MQ=25)
    gAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCt                                                                                                                                                 <  1:1642154/150‑1 (MQ=25)
      aCTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCTCAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCtgt                                                                                                                                               <  1:1827021/150‑1 (MQ=255)
      aCTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCtgt                                                                                                                                               <  1:4691836/150‑1 (MQ=255)
       cTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGGCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCtgtg                                                                                                                                              <  2:2046765/150‑1 (MQ=255)
       cTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCAACCTTGAATGTGGCGCTATACCtgtg                                                                                                                                              <  2:4169876/150‑1 (MQ=25)
       cTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCAACCTTGAATGTGGCGCTATACCtgtg                                                                                                                                              <  1:852961/150‑1 (MQ=25)
       cTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCtgtg                                                                                                                                              <  1:5408898/150‑1 (MQ=255)
       cTGGTCGAGGGCCTAGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATAGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCtgtg                                                                                                                                              <  1:1941001/150‑1 (MQ=255)
       cTGATCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCtgtg                                                                                                                                              <  1:907753/150‑1 (MQ=255)
        tGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTgg                                                                                                                                             <  2:702849/150‑1 (MQ=25)
        tGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTgg                                                                                                                                             <  2:2759456/150‑1 (MQ=25)
        tGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTATACCTGTgg                                                                                                                                             <  1:4745017/150‑1 (MQ=255)
        tGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTgg                                                                                                                                             <  2:1966255/150‑1 (MQ=255)
        tGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTgg                                                                                                                                             <  1:3179612/150‑1 (MQ=255)
            cGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGc                                                                                                                                         <  1:3117337/150‑1 (MQ=255)
             gAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCg                                                                                                                                        >  2:2695796/1‑150 (MQ=255)
              aGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCgg                                                                                                                                       <  1:4313747/150‑1 (MQ=21)
              aGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCgg                                                                                                                                       <  2:3316510/150‑1 (MQ=255)
               gggCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCtgt                                                                                                                                               >  2:5279322/1‑141 (MQ=255)
               gggCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCtgt                                                                                                                                               <  1:5279322/141‑1 (MQ=255)
                 gCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTgg                                                                                                                                             >  1:1352024/1‑141 (MQ=255)
                 gCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTgg                                                                                                                                             <  2:1352024/141‑1 (MQ=255)
                 gCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACAGTCAAGGTTGGCGCAACACCAGTGGGAGCGGGCg                                                                                                                                    >  2:415089/1‑150 (MQ=255)
                    tGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCtgtg                                                                                                                                              <  2:4326799/137‑1 (MQ=21)
                    tGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCtgtg                                                                                                                                              >  1:4326799/1‑137 (MQ=21)
                    tGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATATCGCCACATTCAAGGTTGGCGCTAAACCTGTGGGAGCGGGCGTGc                                                                                                                                 >  2:874197/1‑150 (MQ=255)
                     ggAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCGGGCGTGcc                                                                                                                                <  2:3504355/150‑1 (MQ=21)
                              ggCTTCGGTGGGGGCTGTGGTGGTGGTTGGGCATGGCACGGATGTGTGGTGGTTTTGCGGGCCTGTCCGCGGGCACGCCCGCGCCCACAGGTATCGCGCCCCATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGca                                                                                                                       <  2:3932622/149‑1 (MQ=37)
                               gCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTac                                                                                                                                                     <  1:1067177/119‑1 (MQ=255)
                               gCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTac                                                                                                                                                     >  2:1067177/1‑119 (MQ=255)
                                 ttAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGATACCGCCACAATCAAGGTTGGCGCAACACCGGTGGGAGCCGGCGGGCCCCCGCAGcactc                                                                                                                    >  2:206568/1‑150 (MQ=38)
                                 ttAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCGCCCACAGGGATACCCCCACAATCAAAGTTGGCGCTACACCAGGGGGAGCGGGCGTGCCCCCCAAGcaggc                                                                                                                    >  2:206346/1‑148 (MQ=255)
                                     gTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCGGTGGGAGCGGGCGTGCCCGCGAAGCAGTCCGTg                                                                                                                >  2:3112410/1‑150 (MQ=255)
                                               tgGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGt                                                                                                                                   >  1:4522491/1‑121 (MQ=255)
                                               tgGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGt                                                                                                                                   <  2:4522491/121‑1 (MQ=255)
                                                  aGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCGCCCACAGGTATAGCGCCCCATTCAAGGTTGGCGCTACACCGGTGGGAGCGGGCGTGCCCGCGGAGCCGGCAGGGAGGTGTCAGCGGc                                                                                                   >  2:5203406/1‑150 (MQ=38)
                                                              aTGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTACCCCAGGTATAACGCCCCATTCAAGGGTGGCCCAACACCAGTGGGAGCGGGCGGGCCCCCGAAGCAGTCAGGGGGGTCTCAGCGGCCGTggccggcaa                                                                                       >  2:3226459/1‑149 (MQ=255)
                                                               tGGCAAGGAGGTGGGGTGGTTTTGCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATc                                                                                      <  2:2213043/150‑1 (MQ=255)
                                                                    tggttgtgtgGTGTTTTTGAGGGCCGCTTCGCGGGCCCGCCGGCTCCCCCCGTTATCGCCCCACATTCACGGGTGGCGCTACCCCGGTGGGCGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCAcgc                                                                                 <  2:590768/146‑1 (MQ=255)
                                                                    cggctgtgggGGGGTTGGGCGGGCCGCGGCGCGGGCACGCCCGCGCCCCCCGGTATCGCGCCACATTCCCGGTGGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCAcgc                                                                                 <  2:4497403/146‑1 (MQ=25)
                                                                               tggtTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATTGCGCCACAATGAAGGTTGGCGCTACACCAGTGGGAGCGGGCGTGCCCGCGAAGCAGTCCGTGAGGTCAGAGCGGCCTTGGCAGGCATACACGCGCAACCAccgc                                                                      >  2:3732144/1‑150 (MQ=37)
                                                                                                 cgcgGCCCCCCCCGCACCCACTGTTACACCCCCACATTAAAGCTTGGCGATACACCTGTGGCAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCt                                                    <  2:2534083/150‑1 (MQ=255)
                                                                                                    gggCACGCCCGCTCCCACTGGTATCGCGCCACATTCAAGGTGGGCGCTACCCCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGc                                                 <  2:3199899/150‑1 (MQ=255)
                                                                                                              gCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGa                                       >  2:1853134/1‑150 (MQ=255)
                                                                                                                tCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCCGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                                     >  1:5018094/1‑150 (MQ=255)
                                                                                                                tCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                                     >  1:344233/1‑150 (MQ=255)
                                                                                                                 cccACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                                     >  2:4933948/1‑149 (MQ=25)
                                                                                                                 cccACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                                     <  1:4933948/149‑1 (MQ=25)
                                                                                                                  ccACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGc                                                 <  2:4509781/136‑1 (MQ=255)
                                                                                                                  ccACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGc                                                 >  1:4509781/1‑136 (MQ=255)
                                                                                                                  ccACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACga                                   >  2:4409295/1‑150 (MQ=255)
                                                                                                                  ccACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACga                                   >  2:2133483/1‑150 (MQ=255)
                                                                                                                  ccACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACga                                   >  2:2609708/1‑150 (MQ=255)
                                                                                                                  ccACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGGGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACga                                   >  1:5303490/1‑150 (MQ=255)
                                                                                                                   cacaGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACgag                                  >  1:5211265/1‑150 (MQ=25)
                                                                                                                   cacaGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCGGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACgag                                  >  2:3383189/1‑150 (MQ=255)
                                                                                                                    acaGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACgaga                                 >  2:1276230/1‑150 (MQ=25)
                                                                                                                    acaGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACgaga                                 >  2:587017/1‑150 (MQ=25)
                                                                                                                    acaGGTGTAGCGCCAACCTTGAATGTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACga                                   <  1:2894729/148‑1 (MQ=25)
                                                                                                                    acaGGTGTAGCGCCAACCTTGAATGTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACga                                   >  2:2894729/1‑148 (MQ=25)
                                                                                                                    acaGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACgaga                                 >  2:4308889/1‑150 (MQ=255)
                                                                                                                     caGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGAt                                >  2:2582751/1‑150 (MQ=255)
                                                                                                                      aGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATc                               >  2:805506/1‑150 (MQ=255)
                                                                                                                       ggTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCa                              >  1:5356224/1‑150 (MQ=255)
                                                                                                                       ggTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCa                              >  1:5356306/1‑150 (MQ=255)
                                                                                                                        gTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGaa                                                                                                                          <  1:398270/57‑1 (MQ=255)
                                                                                                                        gTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGaa                                                                                                                          >  2:398270/1‑57 (MQ=255)
                                                                                                                        gTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCAc                             >  1:3550532/1‑150 (MQ=255)
                                                                                                                        gTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCAc                             >  2:3287802/1‑150 (MQ=255)
                                                                                                                         tataGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCAcc                            >  2:3059673/1‑150 (MQ=255)
                                                                                                                          ataGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCa                           >  1:2747970/1‑150 (MQ=255)
                                                                                                                           tcgcgcCACTTTAAAGGTTGGCGCTACCCCTGTGGGTGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAg                          <  1:3569542/148‑1 (MQ=255)
                                                                                                                           taGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAg                          >  1:109408/1‑150 (MQ=255)
                                                                                                                           taGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAg                          >  1:4421590/1‑150 (MQ=255)
                                                                                                                              cgcCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTa                       >  2:1966920/1‑150 (MQ=255)
                                                                                                                                ccACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGcagt                                                                                                                     >  2:1956324/1‑54 (MQ=255)
                                                                                                                                ccACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGcagt                                                                                                                     <  1:1956324/54‑1 (MQ=255)
                                                                                                                                ccACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGa                     >  2:2194223/1‑150 (MQ=255)
                                                                                                                                      ttcttggTGGGCGCTACCCCTGTGGGCGGGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACg                <  1:2695796/145‑1 (MQ=255)
                                                                                                                                       tCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGc               >  1:3423109/1‑150 (MQ=255)
                                                                                                                                         aaGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCACCCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCt             >  2:5388602/1‑150 (MQ=255)
                                                                                                                                         aaGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGt                        >  1:1306292/1‑139 (MQ=255)
                                                                                                                                         aaGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGt                        <  2:1306292/139‑1 (MQ=255)
                                                                                                                                         aaGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGAACCAGCATGAACGAGATCACCAGGTAGAACACGCCt             >  2:5373492/1‑150 (MQ=255)
                                                                                                                                           gtttGCCGCTACCTCGGTGTGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTgc           <  1:2042674/148‑1 (MQ=255)
                                                                                                                                            gTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTgcg          >  1:3020325/1‑150 (MQ=255)
                                                                                                                                               ggCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTgcggcg       >  2:775223/1‑150 (MQ=255)
                                                                                                                                                 cgcTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCgaa     <  1:2609708/150‑1 (MQ=255)
                                                                                                                                                  gcTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCgaag    >  2:1730773/1‑150 (MQ=255)
                                                                                                                                                   cTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCgaaga   >  1:1344588/1‑150 (MQ=255)
                                                                                                                                                   cTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCgaaga   >  1:1600202/1‑150 (MQ=255)
                                                                                                                                                   cTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCgaaga   >  1:1600211/1‑150 (MQ=255)
                                                                                                                                                   cTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCgaaga   >  2:4044513/1‑150 (MQ=255)
                                                                                                                                                   cTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCgaaga   >  1:1283780/1‑150 (MQ=255)
                                                                                                                                                   cTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCgaaga   >  2:1541625/1‑150 (MQ=255)
                                                                                                                                                    tACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCgaagaa  >  1:2166721/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
TCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAA  >  NC_002947/337690‑337986

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