Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,238,939 G→C 12.8% intergenic (+69/+59) 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,9390GC12.8% 246.7 / 25.1 86intergenic (+69/+59)uvrB/PP_1975uvrABC excinuclease protein B/EmrB/QacA family drug resistance transporter
Reads supporting (aligned to +/- strand):  ref base G (40/35);  new base C (5/6);  total (45/41)
Fisher's exact test for biased strand distribution p-value = 7.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

CCGCGACCTGGAGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATCCCGACAAAG  >  NC_002947/2238792‑2239088
                                                                                                                                                   |                                                                                                                                                     
ccGCGACCTGGAGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTgcg                                                                                                                                                     >  2:4720103/1‑150 (MQ=255)
ccGCGACCTGGAGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCAGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCCCCGGCTgcg                                                                                                                                                     >  2:4719960/1‑150 (MQ=255)
  gcgAACTGGAGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGcc                                                                                                                                                   >  2:3544046/1‑150 (MQ=255)
  gcgAACTGGAGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGcc                                                                                                                                                   >  2:3543568/1‑150 (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                                                                                                                                            <  1:1356704/150‑4 (MQ=255)
         ggAGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCGCAGCCGgtgcca                                                                                                                                            <  2:1820732/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:3720300/150‑1 (MQ=255)
                     tGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGccc                                                                                                                                <  1:4704856/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                                                                                                                                                >  2:2822127/1‑107 (MQ=255)
                                                gcAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGgt                                                                                                                                                <  1:2822127/107‑1 (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                                                                                                                                              <  2:1861066/85‑2 (MQ=21)
                                                                        cagcTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGACACCGGCGCAGCCGgtgc                                                                                                                                              >  1:1861066/1‑84 (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:272690/1‑150 (MQ=255)
                                                                                                               cccGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCAcc                                      >  1:272730/1‑150 (MQ=255)
                                                                                                               cccGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCAcc                                      >  1:3775483/1‑150 (MQ=255)
                                                                                                                ccGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCa                                                                                                              >  1:343589/1‑77 (MQ=255)
                                                                                                                ccGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCa                                                                                                              <  2:343589/77‑1 (MQ=255)
                                                                                                                  gcgAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTcc                                                                                                                           <  2:4117433/62‑1 (MQ=255)
                                                                                                                  gcgAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTcc                                                                                                                           >  1:4117433/1‑62 (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             >  2:4705851/1‑150 (MQ=255)
                                                                                                                                        tggcACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGa             >  1:4918641/1‑150 (MQ=255)
                                                                                                                                        tggcACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGa             >  1:3878148/1‑150 (MQ=255)
                                                                                                                                        tggcACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGa             >  1:268876/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:5408219/150‑1 (MQ=255)
                                                                                                                                              cGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATCCCGa       <  1:4047830/150‑1 (MQ=255)
                                                                                                                                                  tGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATCCCGACaaa   >  2:267199/1‑150 (MQ=255)
                                                                                                                                                   gcgcGGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATCCCGACAAAg  <  1:4719960/150‑1 (MQ=255)
                                                                                                                                                   gcgcCGGTGTTCGGGGGCAGGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATCCCGACAAAg  <  1:4720103/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCGCGACCTGGAGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATCCCGACAAAG  >  NC_002947/2238792‑2239088

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: