Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,238,938 T→G 14.6% intergenic (+68/+60) uvrB → / ← PP_1975 uvrABC excinuclease protein B/EmrB/QacA family drug resistance transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,238,9380TG14.6% 214.3 / 26.0 87intergenic (+68/+60)uvrB/PP_1975uvrABC excinuclease protein B/EmrB/QacA family drug resistance transporter
Reads supporting (aligned to +/- strand):  ref base T (39/35);  new base G (7/6);  total (46/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

CCGCGACCTGGAGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATCCCGACAAA  >  NC_002947/2238792‑2239087
                                                                                                                                                  |                                                                                                                                                     
ccGCGACCTGGAGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTgcg                                                                                                                                                    >  2:4720103/1‑150 (MQ=255)
ccGCGACCTGGAGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCAGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCCCCGGCTgcg                                                                                                                                                    >  2:4719960/1‑150 (MQ=255)
  gcgAACTGGAGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGcc                                                                                                                                                  >  2:3543568/1‑150 (MQ=255)
  gcgAACTGGAGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGcc                                                                                                                                                  >  2:3544046/1‑150 (MQ=255)
  gcgAACTGGAGTTAGAACCTGCGGCGCAGATCCGGCCCGAGATTGCGCAGCTCCGGGAGCGGTAGATCAGCAGCTGAAGCCCGGCGCGGGCCCTGTGGGAGCGGGCGAGCCCGCGAAGAAGTCAACGCGGTGGCAGGCACCGCCtccccc                                                                                                                                                  >  2:3543798/1‑145 (MQ=255)
     aCCTGGAGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCGGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGgt                                                                                                                                               <  1:5451280/150‑1 (MQ=255)
     aCCTGGAGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGgt                                                                                                                                               <  1:87494/150‑1 (MQ=255)
     aCCTGGAGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGgt                                                                                                                                               <  2:3599867/150‑1 (MQ=255)
       cTGGAGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCGAACACCGGCGCAGCCGgtgc                                                                                                                                             <  1:4226984/150‑2 (MQ=39)
       cTGGAGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCGAACACCGGCGCAGCCGgtgc                                                                                                                                             <  1:2665754/150‑2 (MQ=39)
        tGGAGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGtt                                                                                                                                            <  2:683498/150‑1 (MQ=255)
        tGGAGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCGGTGGGAGCGGGCGGGCCCGCGCAGAAATAAACGCCGGGGCGGGCCCCCGCGGCGCCCGTGtt                                                                                                                                            >  1:2316241/1‑150 (MQ=255)
         ggAGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCGCAGCCGgtgcca                                                                                                                                           <  2:1820732/150‑4 (MQ=255)
         ggAGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCGCAGCCGgtgcca                                                                                                                                           <  1:1356704/150‑4 (MQ=255)
           aGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGGGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCGAACACCGGCGCAGCCGgtgccagc                                                                                                                                         <  2:3743986/150‑6 (MQ=38)
           aGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTcgc                                                                                                                                         <  2:1614917/150‑1 (MQ=255)
             ttCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCgg                                                                                                                                       <  1:4807708/150‑1 (MQ=255)
             ttCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCgg                                                                                                                                       <  1:4103815/150‑1 (MQ=255)
                gAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCCACGCGGGGGCTGGCCCCGGCTGCGCCGGTGGTCGCGGGcc                                                                                                                                    >  2:1498892/1‑149 (MQ=255)
                gAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCa                                                                                                                                    >  2:4826052/1‑150 (MQ=255)
                   gCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCg                                                                                                                                                 >  2:4305940/1‑134 (MQ=255)
                   gCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCg                                                                                                                                                 <  1:4305940/134‑1 (MQ=255)
                   gCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGc                                                                                                                                 <  2:2894490/150‑1 (MQ=255)
                     tGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGccc                                                                                                                               <  1:3720251/150‑1 (MQ=255)
                     tGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGccc                                                                                                                               <  1:4704856/150‑1 (MQ=255)
                     tGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGccc                                                                                                                               <  1:3720300/150‑1 (MQ=255)
                            cAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCca                                                                                                                        <  2:840283/150‑1 (MQ=255)
                                tGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAgg                                                                                                                    <  1:3390566/150‑1 (MQ=255)
                                tGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAgg                                                                                                                    <  1:765799/150‑1 (MQ=255)
                                   gTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGAc                                                                                                                 >  2:2774361/1‑150 (MQ=255)
                                   gTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGGGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGAc                                                                                                                 <  2:4366374/150‑1 (MQ=255)
                                     gACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCGGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCGGGCCCCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCg                                                                                                               <  2:2358380/150‑1 (MQ=255)
                                        gagaTTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGCAGAATTCAACGCGGTGGCTGGCCCCGGCTGCGCCGGTGTTCGCGGGCAGGCCCGCTCCCCCCGGGACCGccg                                                                                                            >  1:5366612/1‑148 (MQ=37)
                                           aTTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGGGCCCGCGAAGAATTCAAAGCGGTGGCTGGCACCCGCTGCGCCGGTGTTAGCGGGCATTCCCCCTCCCACACGGACCGCaactg                                                                                                         >  1:4194058/1‑150 (MQ=255)
                                            ttGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGGGCGGGCGGGCCCCCGAAGAAATAAAAGCCGTGGCGGGCACCGGCGGCGCCCGGGTTGGCGGGGAGGCCCCCGCCCACACGGGCCCCagcaga                                                                                                        >  1:253805/1‑145 (MQ=255)
                                             tGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCatctgat                                                                                                       <  2:3557232/150‑1 (MQ=255)
                                                gcAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGgt                                                                                                                                               <  1:2822127/107‑1 (MQ=255)
                                                gcAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGgt                                                                                                                                               >  2:2822127/1‑107 (MQ=255)
                                                  aGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGTATTATTGAAAGCGGTTGCTGGCACCGGGTGGGCCGGTTTTGGGGGGGATTGCCGCTGCCACAGGGACCGCatatgatattta                                                                                                  >  1:779955/1‑145 (MQ=255)
                                                  aGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTa                                                                                                  >  2:4913656/1‑150 (MQ=255)
                                                       cGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCGGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGc                                                                                             <  1:5340736/150‑1 (MQ=255)
                                                          gAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGGGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTc                                                                                          <  2:722890/150‑1 (MQ=255)
                                                          gAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTc                                                                                          >  2:1052297/1‑150 (MQ=255)
                                                                    tCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTc                                                                                <  1:5448136/150‑1 (MQ=255)
                                                                    tCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCAACTGAGAGGTAATAGCCTCAAAGGGCTTc                                                                                >  1:425564/1‑150 (MQ=25)
                                                                        cagcTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGACACCGGCGCAGCCGgtgc                                                                                                                                             >  1:1861066/1‑84 (MQ=21)
                                                                        cagcTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGACACCGGCGCAGCCGgtgc                                                                                                                                             <  2:1861066/85‑2 (MQ=21)
                                                                           cTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGGGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGcc                                                                         <  2:645629/150‑1 (MQ=255)
                                                                                gCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGCAGAAGTCAAAGCGGTGGCTGGCACCGGCTGCGCCGGTGGTCGCGGGCCTGCCCCCTCCCACAGGGACCGCCTCAGATCGGTAATAGCCACAAAGGGCGTACCCAGccaaaca                                                                    >  2:11180/1‑145 (MQ=18)
                                                                                       gcgGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTcc                                                             <  1:1763045/150‑1 (MQ=255)
                                                                                       gcgGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTcc                                                             <  2:4518973/150‑1 (MQ=255)
                                                                                             ccTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAAc                                                       >  2:5253402/1‑150 (MQ=255)
                                                                                               tgtgGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACtt                                                     >  1:2570131/1‑150 (MQ=255)
                                                                                                 tgGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGc                                                   <  2:3635897/150‑1 (MQ=255)
                                                                                                  gggAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCg                                                  <  1:2169486/150‑1 (MQ=255)
                                                                                                    gAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCggg                                                >  2:1052418/1‑150 (MQ=255)
                                                                                                      gCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCatctgat                                                                                                       <  1:3453741/93‑1 (MQ=255)
                                                                                                      gCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCatctgat                                                                                                       >  2:3453741/1‑93 (MQ=255)
                                                                                                         ggCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTcagc                                           <  2:3173286/150‑1 (MQ=255)
                                                                                                               cccGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCAcc                                     >  1:3775483/1‑150 (MQ=255)
                                                                                                               cccGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCAcc                                     >  1:272730/1‑150 (MQ=255)
                                                                                                               cccGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCAcc                                     >  1:272690/1‑150 (MQ=255)
                                                                                                                ccGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCa                                                                                                             <  2:343589/77‑1 (MQ=255)
                                                                                                                ccGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCa                                                                                                             >  1:343589/1‑77 (MQ=255)
                                                                                                                  gcgAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTcc                                                                                                                          >  1:4117433/1‑62 (MQ=255)
                                                                                                                  gcgAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTcc                                                                                                                          <  2:4117433/62‑1 (MQ=255)
                                                                                                                        aaTTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATg                            >  2:1510213/1‑150 (MQ=255)
                                                                                                                          ttCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCt                          >  1:1154740/1‑150 (MQ=255)
                                                                                                                           tCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTc                         >  2:3008493/1‑150 (MQ=255)
                                                                                                                             aaCGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCAc                       >  1:4920317/1‑150 (MQ=255)
                                                                                                                             aaCGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCAc                       >  2:1632782/1‑150 (MQ=255)
                                                                                                                              aCGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCAcc                      <  1:4002414/150‑1 (MQ=255)
                                                                                                                                gcgGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGc                    <  2:620391/150‑1 (MQ=255)
                                                                                                                                gcgGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCGTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCAAGCTCACCGc                    >  1:70543/1‑150 (MQ=255)
                                                                                                                                    tggctggcACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAAc                >  2:3546340/1‑150 (MQ=255)
                                                                                                                                     cgcgaacACCGGCGCAGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAAcc               >  2:4455606/7‑150 (MQ=39)
                                                                                                                                        tggcACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGa            >  1:268876/1‑150 (MQ=255)
                                                                                                                                        tggcACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGa            >  1:3878148/1‑150 (MQ=255)
                                                                                                                                        tggcACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGa            >  2:4705851/1‑150 (MQ=255)
                                                                                                                                        tggcACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGa            >  1:4918641/1‑150 (MQ=255)
                                                                                                                                        gaacACCGGCGCAGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGa            >  2:2120558/4‑150 (MQ=255)
                                                                                                                                         aacACCGGCGCAGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGAt           >  2:1417884/3‑150 (MQ=255)
                                                                                                                                          gcACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATc          >  1:2890819/1‑150 (MQ=255)
                                                                                                                                           cACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATcc         >  1:143129/1‑150 (MQ=255)
                                                                                                                                           cACCGGCGCAGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATcc         >  2:3684503/1‑150 (MQ=255)
                                                                                                                                              cGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGTCAACCCGATCCCGa      >  2:1039364/1‑150 (MQ=255)
                                                                                                                                              cGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATCCCGa      <  1:4047830/150‑1 (MQ=255)
                                                                                                                                              cGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATCCCGa      <  1:5408219/150‑1 (MQ=255)
                                                                                                                                                  tGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATCCCGACaaa  >  2:267199/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CCGCGACCTGGAGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATCCCGACAAA  >  NC_002947/2238792‑2239087

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: