Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,065,170 T→C 22.6% intergenic (+62/‑145) PP_1838 → / → PP_1839 membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,065,1700TC22.6% 187.1 / 45.8 95intergenic (+62/‑145)PP_1838/PP_1839membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (36/37);  new base C (6/15);  total (42/53)
Fisher's exact test for biased strand distribution p-value = 1.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.50e-01

TCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGA  >  NC_002947/2065023‑2065318
                                                                                                                                                    |                                                                                                                                                    
tCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTTCACAGGGttag                                                                                                                                                    >  2:138342/1‑148 (MQ=255)
tCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCg                                                                                                                                                    >  1:4505081/1‑150 (MQ=255)
tCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCg                                                                                                                                                    >  1:4444158/1‑150 (MQ=255)
tCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTACCACAGGCCTTGCACAGGGTTCg                                                                                                                                                    >  2:2093762/1‑150 (MQ=255)
 cAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCgg                                                                                                                                                   <  1:1192747/150‑1 (MQ=255)
 cAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCgg                                                                                                                                                   <  1:3950903/150‑1 (MQ=255)
  aagaagCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGa                                                                                                                                                  >  2:4490103/1‑150 (MQ=255)
      agCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCtt                                                                                                                                              <  1:3275270/150‑1 (MQ=255)
      agCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCtt                                                                                                                                              <  1:3140487/150‑1 (MQ=255)
      agCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCtt                                                                                                                                              <  1:613688/150‑1 (MQ=255)
      agCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCtt                                                                                                                                              <  1:3100260/150‑1 (MQ=255)
       gCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTtg                                                                                                                                             <  1:474917/150‑1 (MQ=255)
         tGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTtgtg                                                                                                                                           <  1:1360081/150‑1 (MQ=255)
             aGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAgg                                                                                                                                       <  2:4031745/150‑1 (MQ=255)
             aGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAgg                                                                                                                                       <  1:3258848/150‑1 (MQ=255)
             aGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGACCTGCACAAGGTCCGAACCCTGTGCaat                                                                                                                                       <  2:2191323/150‑3 (MQ=37)
              gCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCaatg                                                                                                                                      <  1:1614006/150‑4 (MQ=25)
               cGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCaatgc                                                                                                                                     <  1:3328640/150‑5 (MQ=25)
               cGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCACTGCACAAGGTCCGAACCCTGTGCaatgc                                                                                                                                     <  1:4276591/150‑5 (MQ=37)
                gTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCAATCGCGGGCTTGCCCGCGCCCCCCGGCAATGCCCAGGGTTCGGACCTTGTGCAGGtgc                                                                                                                                    >  2:1053812/1‑148 (MQ=255)
                    atCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCtgtg                                                                                                                                <  2:4695404/150‑1 (MQ=21)
                    atCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCtgtg                                                                                                                                <  2:1735466/150‑1 (MQ=21)
                    atCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCtgtg                                                                                                                                <  1:3062543/150‑1 (MQ=255)
                    atCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAAGGTCCGAACCCTGTGCAATGCCtgtg                                                                                                                                <  2:3513027/150‑1 (MQ=37)
                     tCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCGTTGCACAGGGTTCGGACCTTGTGCAGGTCCtgtgg                                                                                                                               >  2:4516830/1‑149 (MQ=255)
                      ccAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGc                                                                                                                               <  1:3446003/149‑1 (MQ=255)
                      ccAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGc                                                                                                                               >  2:3446003/1‑149 (MQ=255)
                       cAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCtgtgg                                                                                                                               >  2:1209699/1‑147 (MQ=18)
                       cAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCtgtgg                                                                                                                               <  1:1209699/148‑2 (MQ=18)
                        ggACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGGGTGCCCGCGCCCCCAGGCAGTGGGCACGGTTCGGGCCATTTGCAGGGCCGGGGCggg                                                                                                                            >  1:1247884/2‑148 (MQ=255)
                        aGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGGGTGGCCGCGCCCCCAGGCCATGGACAGGGTTGGGGCCGTGTGCCGGGCCGGGGCgcg                                                                                                                            >  1:3458291/1‑148 (MQ=255)
                         gACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGc                                                                                                                           <  1:3469921/150‑1 (MQ=18)
                         gACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGCGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGc                                                                                                                           <  2:4729617/150‑1 (MQ=18)
                           ccGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCgg                                                                                                                         <  2:3794759/150‑1 (MQ=18)
                           ccGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCgg                                                                                                                         <  2:4478153/150‑1 (MQ=18)
                           ccGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCgg                                                                                                                         <  1:1034005/150‑1 (MQ=255)
                           ccGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCgg                                                                                                                         <  2:4174862/150‑1 (MQ=255)
                               gACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGGTCGGACCGTGTGCAGGGACGGGGCGAGCGGGccc                                                                                                                     >  1:3104254/1‑148 (MQ=39)
                                   aaGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCGTTCCCGCTCCCACAGGCCTTGGACAGGGTTCGGACCTTGTGCAGGTGCAGTGCGAGCGGGCCAGccc                                                                                                                 >  2:1427903/1‑150 (MQ=255)
                                     gagCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGATTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCcgc                                                                                                               <  1:2677264/150‑1 (MQ=255)
                                      agCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCGTTGCCCCGCCCACAGGCAGTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGGGGGCGGGCCAGCCcgcg                                                                                                              >  2:941383/1‑150 (MQ=255)
                                          tgctgGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCtt                                                                                                                                              >  1:5182369/1‑114 (MQ=255)
                                          tgctgGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCtt                                                                                                                                              <  2:5182369/114‑1 (MQ=255)
                                           gctgGATAACACCGCCGGCGTGGTGCCCGGCACACCTAATTAAACCGGGGCCGCGTTGCGGCCCATTCGGGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCTGGTCGTGTGCGAGCGGGCAAGCCCGCGAATgg                                                                                                         <  2:1448560/150‑1 (MQ=38)
                                           gctgGATAACACCGCCGGCGTGGTGCAGGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACGTTGTGCAGGTCCTGTGCGAGGGGGCAAGCCCGCGAATgg                                                                                                         <  2:1448591/150‑1 (MQ=255)
                                             tgGATAACACCGCCGGCGTGGTGCATGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTAGCGGGGTTTGCCGCTCCCAAAGGGATTTGACATGGTTTGGTACTTTTTGAGGTGCTGTGGGAGCGGGGAATCCCGCGTATGGGc                                                                                                       >  2:1510131/1‑150 (MQ=255)
                                                     acCGCCGGCGTGGTGCACGGCACACCCAAGTAAACCGGGGCCGCGTTGCGGCCCATACGCGGGCTTTACCCCACCCCAAGCCCAAGCACAGCGGTCGGACCATGTGCAGGTACTGAGCGGGCGGGCAAACCCCAGCCTTGGCCCCGcaca                                                                                               >  2:2784536/1‑150 (MQ=255)
                                                                 gtgCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGGCCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGAc                                                                                   >  2:3986007/1‑150 (MQ=255)
                                                                           cacCGAAGTAAACGGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCTGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTAt                                                                         <  2:937462/150‑1 (MQ=255)
                                                                                   tAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGccc                                                                 >  2:5433581/1‑150 (MQ=255)
                                                                                       ccGGGGCCGCGTTGCGGCCCATTCGCGGGCTTTCCCGCTCCCACACGCAGTGCACAGGGTTCGGACCTTGTGCAGGTCCGGTGCGAGCGTGCAAGCCCGCGACTGGGCCCCCCACAATTTTCCCGGGCTGAACAATATACCCGccccccg                                                             >  2:3743792/1‑147 (MQ=255)
                                                                                       ccGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCCCGCCCCCAGGCCGTGCGCCGGGGTCGGGCCGTGTGCAGGGCCGGGGCGCGCGGGCCCGCCCGCGCCTGGGCCCCGCACCAGTTTACCCGCCGGAACAGTATAccccccccccg                                                             >  2:3411198/1‑147 (MQ=255)
                                                                                           ggCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACCGGCATTGCACAGGGTTCGGACCTTGTGCTGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCa                                                         <  2:3674742/150‑1 (MQ=255)
                                                                                           gcccggtttGCCGCCCTTGCGGGGGCTGCCCCCCACCCCCCGGCATGGCACGGGGTTCGGACCTGGCGCCGGTCGGGGGCCCGGGGCCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCa                                                         <  1:5154777/143‑1 (MQ=255)
                                                                                            gCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAg                                                        <  2:3949793/150‑1 (MQ=255)
                                                                                                   tggggCCCATGCGCGGGGTGCCCCGCCCCCACCGTATTTGCACTGGGTTCGGACGTTGTGCTGGTCGTGGCCGCGCGGCCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCa                                                 <  1:1427903/147‑1 (MQ=255)
                                                                                                    ggggCCAATTCGGGGGGTGGCCCGCTCCCCCTGGCTTTGCACTGGGTTCGGACCTTGGGCTGGTCGTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAg                                                <  1:383873/148‑1 (MQ=255)
                                                                                                       gCCCTTCCGCGGGCTGGCCCGCCCCCACCGCCTTTGCCCCGGTTTCGGACCTTGGCCTGGGCCGGTGCGAGCGGCCAACCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCa                                             <  1:3893690/150‑1 (MQ=255)
                                                                                                       gCCCATTCGCGGGCTTGCCCGCTCCCAACGGCCTTGCCCAGGGGTCCGACCTTGTGCAGGTCCTGTGCGAGCGGGCAATCCCGCGAAAGGGCCCCCCACAAATTTCCCCGACTGAACAATAAACCCGCCCCACCGTCACTTATTCAGGCa                                             >  2:4571565/1‑150 (MQ=25)
                                                                                                        cccATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCAt                                            <  2:3049676/150‑1 (MQ=255)
                                                                                                        cccATGGGCGGGCTGGCCCGCCCCCCCCGGCTGTGCCGGGGGTTGGGCCCTTGTGCCGGCCGTGGGCGCGGGGGCACGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCAt                                            <  1:3411198/150‑1 (MQ=255)
                                                                                                         ccATTCGCGGGCTTGCCCGCTCCCCCAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATc                                           <  2:4820741/150‑1 (MQ=255)
                                                                                                         ccATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTc                                                                                                                                     >  2:475741/1‑60 (MQ=255)
                                                                                                         ccATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTc                                                                                                                                     <  1:475741/60‑1 (MQ=255)
                                                                                                                          ccGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCg                           >  2:1694587/1‑150 (MQ=18)
                                                                                                                          ccGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCg                           >  2:4437681/1‑150 (MQ=18)
                                                                                                                          ccGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCg                           >  1:4068295/1‑150 (MQ=18)
                                                                                                                          ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCg                           >  2:1680653/1‑150 (MQ=255)
                                                                                                                            gCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGtt                         >  2:3689783/1‑150 (MQ=255)
                                                                                                                               cccACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCTAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGAc                      <  2:5103324/150‑1 (MQ=255)
                                                                                                                                ccACAGGCATTGCACAGGGTTCGTACCTTGTGCAGGACCTGTGCGAGCGGGCAAGCCCTCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACt                     <  2:427534/150‑1 (MQ=255)
                                                                                                                                ccACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACt                     >  1:885737/1‑150 (MQ=255)
                                                                                                                                ccACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACt                     >  2:5016117/1‑150 (MQ=255)
                                                                                                                                 cacaGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTg                    >  2:4868502/1‑150 (MQ=255)
                                                                                                                                  acaGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTga                   >  2:5221035/1‑150 (MQ=255)
                                                                                                                                  acaGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTga                   >  1:3335894/1‑150 (MQ=255)
                                                                                                                                  acaGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTga                   >  1:45022/1‑150 (MQ=255)
                                                                                                                                    aggacctGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTgaga                 >  2:180858/7‑150 (MQ=21)
                                                                                                                                    agcaaTTCCAAAGTGTTCTCACATTGTCCACGACCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAACGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTgaga                 <  2:1005684/146‑1 (MQ=255)
                                                                                                                                    aGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTgaga                 >  2:740484/1‑150 (MQ=255)
                                                                                                                                     ggcttggCACTGGGTTGGGACCTTGTCCCGCTCGTGTGCGCGGGGGCACGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGagag                <  1:1344939/144‑1 (MQ=255)
                                                                                                                                      gctttGCCCCGGGTCCGGCCCTGGGGCCGGTCCGGGGCGCGCGGGCACGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAgg               <  1:2549648/147‑1 (MQ=255)
                                                                                                                                      gacctGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAgg               >  1:4507160/5‑150 (MQ=25)
                                                                                                                                      gCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAgg               <  2:3390125/150‑1 (MQ=255)
                                                                                                                                      gCATGGCACTGGGGCCGGCCCTTGGGCCGGGCGGGGGCGAGCGGGCACCCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAgg               <  1:2549657/150‑1 (MQ=255)
                                                                                                                                             acaGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGtt        >  2:2187062/1‑150 (MQ=255)
                                                                                                                                             acaGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGtt        >  1:317478/1‑150 (MQ=255)
                                                                                                                                             acaGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGtt        >  2:3656386/1‑150 (MQ=255)
                                                                                                                                             acaGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGtt        >  2:2763846/1‑150 (MQ=255)
                                                                                                                                             acaGGGTACGGACCTGGTGCAGGTCCTGTGCGACCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGtt        <  1:3099403/150‑1 (MQ=255)
                                                                                                                                               cgggTTCGGCCCTTGTGCAGGACCTGTGCGCGGGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCa      <  1:138342/149‑1 (MQ=255)
                                                                                                                                               aGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCGAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCa      >  1:3062762/1‑150 (MQ=255)
                                                                                                                                                 ggTTCGGCCCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAAt    <  2:1247884/150‑1 (MQ=255)
                                                                                                                                                 ggTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAAt    >  1:351084/1‑150 (MQ=255)
                                                                                                                                                 ggTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAAt    <  2:1247716/150‑1 (MQ=255)
                                                                                                                                                 ggTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAAt    <  2:4068295/150‑1 (MQ=37)
                                                                                                                                                   ttCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGa  >  2:4528805/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGA  >  NC_002947/2065023‑2065318

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: