Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,843,106 G→T 5.7% intergenic (+62/+42) PP_1642 → / ← PP_1643 hypothetical protein/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,843,1060GT 5.7% 886.0 / 28.3 263intergenic (+62/+42)PP_1642/PP_1643hypothetical protein/membrane protein
Reads supporting (aligned to +/- strand):  ref base G (152/96);  new base T (8/7);  total (160/103)
Fisher's exact test for biased strand distribution p-value = 5.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 263 total aligned reads displayed.

GCAGGCGTTCAACCAGTGGCATGCAGAGAGCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCC  >  NC_002947/1842960‑1843252
                                                                                                                                                  |                                                                                                                                                  
  agGCGTTCAACCAGTGGCATGCAGAGAGCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCtgt                                                                                                                                               <  2:466154/150‑1 (MQ=255)
     cGTTCAACCAGTGGCATGCAGAGAGCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCa                                                                                                                                            <  1:547482/150‑1 (MQ=255)
     cGTTCAACCAGTGGCATGCAGAGAGCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCa                                                                                                                                            <  2:5687266/150‑1 (MQ=255)
     cGTTCAACCAGTGGCATGCAGAGAGCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCa                                                                                                                                            <  2:558478/150‑1 (MQ=255)
     cGTTCAACCAGTGGCATGCAGAGAGCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCa                                                                                                                                            <  2:856850/150‑1 (MQ=255)
       ttCAACCAGTGGCATGCAGAGAGCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGt                                                                                                                                          >  2:3467341/1‑150 (MQ=255)
          aaCCAGTGGCATGCAGAGAGCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTAcc                                                                                                                                       >  2:6795892/1‑150 (MQ=255)
          aaCCAGTGGCATGCAGAGAGCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTAcc                                                                                                                                       <  1:3453471/150‑1 (MQ=255)
          aaCCAGTGGCATGCAGAGAGCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGGTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCAGTGCAGTAcc                                                                                                                                       >  2:2427307/1‑150 (MQ=255)
          aaCCAGTGGCATGCAGAGAGCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCGTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCCCTGAAGGCCGGTGCCGTAcc                                                                                                                                       >  1:258807/1‑150 (MQ=255)
          aaCCAGTGACATGCAGAGAGCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTAcc                                                                                                                                       >  2:3694214/1‑150 (MQ=255)
           aCCAGTGGCATGCAGAGAGTGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCt                                                                                                                                      >  1:3093769/1‑150 (MQ=255)
           aCCAGTGGCATGCAGAGAGCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCt                                                                                                                                      <  2:735634/150‑1 (MQ=255)
            ccAGTGGCATGCAGAGAGCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCCGTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACGTAAGGCCTGTGCAGTACCtg                                                                                                                                     >  2:4263575/1‑150 (MQ=255)
             cAGTGGCATGCAGAGAGCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTAAAGGCCTGTGCAGTACCtgt                                                                                                                                    >  1:3631436/1‑150 (MQ=255)
               gTGGCATGCAGAGAGCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATACCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTgg                                                                                                                                  <  2:841119/150‑1 (MQ=255)
                   cATGCAGAGAGCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGGTACTGCACAGGCCTTAAGTGCGGTGCGATTACTGTGGGAGc                                                                                                                              <  1:4593136/150‑1 (MQ=18)
                    aTGCAGAGAGCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGGCCTTAAGTGCGGTGCGATTACTGTGGGAGCg                                                                                                                             <  1:3797707/150‑1 (MQ=18)
                     tGCAGAGAGCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCgg                                                                                                                            >  1:5610205/1‑150 (MQ=255)
                         gagagCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGt                                                                                                                        >  2:4569810/1‑150 (MQ=255)
                           gagCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGGCCTTAAGTGCGGTGCGATTACTGTGGGAGCGGGCATGc                                                                                                                      <  1:5201123/150‑1 (MQ=17)
                            agCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGcc                                                                                                                     >  1:5807959/1‑150 (MQ=255)
                             gcgggcggCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACCGGTACTGCACAGGCCTTAAGTGCGGTGCGATTACTGTGGGAGCGGGCATGccc                                                                                                                    <  2:5512140/150‑1 (MQ=17)
                               gggcggCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCcgc                                                                                                                  >  1:3161977/1‑150 (MQ=255)
                                     ccGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGacaag                                                                                                            >  1:146628/1‑149 (MQ=255)
                                     ccGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGaaaaa                                                                                                            >  2:1897345/1‑150 (MQ=255)
                                       gCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATc                                                                                                          >  1:2405263/1‑150 (MQ=255)
                                          gtgCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTaa                                                                                                       <  2:1550721/150‑1 (MQ=255)
                                                gTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCt                                                                                                 >  2:6524866/1‑150 (MQ=255)
                                                    tCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTcgc                                                                                             >  2:586149/1‑150 (MQ=255)
                                                      ggTGATGCTGGCGCACCCCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTcgccc                                                                                           <  2:6295506/150‑1 (MQ=255)
                                                       gTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTcgcccg                                                                                          >  2:1010637/1‑150 (MQ=255)
                                                       gTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCGGGGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTcgcccg                                                                                          >  1:103001/1‑150 (MQ=255)
                                                         gATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTcgcccgcc                                                                                        <  2:6573572/150‑1 (MQ=255)
                                                            gCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCTTGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACAATCTAACTGCCTGCGCCCGCCagg                                                                                     >  2:4544102/1‑148 (MQ=255)
                                                                  gcACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTcgcc                                                                                            <  1:3390061/137‑1 (MQ=255)
                                                                     cacCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCa                                                                            >  1:2152774/1‑150 (MQ=255)
                                                                     cacCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCGGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCa                                                                            >  2:6128930/1‑150 (MQ=255)
                                                                        cGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGc                                                                         >  1:3165143/1‑150 (MQ=255)
                                                                        cGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGc                                                                         >  2:2651873/1‑150 (MQ=255)
                                                                         gTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCt                                                                        <  1:4056091/150‑1 (MQ=255)
                                                                          ccaccccttGCGCCTTGTGGTGCTCCCCTGCGCGGGCAGGCCCGCTCCCACAGTAAACGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTg                                                                       <  1:5529289/143‑1 (MQ=38)
                                                                            aaCCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGa                                                                     >  1:3509115/1‑150 (MQ=255)
                                                                               cATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACt                                                                  <  1:5106676/150‑1 (MQ=255)
                                                                                 ttGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTc                                                                <  1:2971387/150‑1 (MQ=255)
                                                                                   gcgCCTTGTGTTGCCCCCATGCGCGGGCAGGCCCGCGCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTCCCTGTGGGCGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCa                                                              <  1:3148111/148‑1 (MQ=255)
                                                                                    aGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAg                                                             <  1:6007813/150‑1 (MQ=255)
                                                                                        ttGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCa                                                         >  1:6674387/1‑150 (MQ=255)
                                                                                         tgtgTTGCTCCCATGCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCAc                                                        <  1:4861196/150‑1 (MQ=255)
                                                                                           tgtTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGt                                                      >  1:3870245/1‑150 (MQ=255)
                                                                                               gCTCCCATTCGCGGGCATGCCCGCTCCCCCAGTAATCGCACCGCACTTAAGGCCGGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAg                                                  <  2:258807/150‑1 (MQ=255)
                                                                                               gCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCGGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACAGCCTTCGCCCGCCAAGCGCCGCGCCTGCTGGAACTTCCAGCGCACGTACAg                                                  >  1:200320/1‑150 (MQ=38)
                                                                                               gCTCCCATTCGCGGGCATGCCCGCTCCCACAGGAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTCCAg                                                  >  2:1168079/1‑150 (MQ=255)
                                                                                                 tCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTa                                                >  2:334137/1‑150 (MQ=255)
                                                                                                 tCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTa                                                >  1:6209431/1‑150 (MQ=255)
                                                                                                  cccATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAg                                               >  2:4304355/1‑150 (MQ=255)
                                                                                                  cccATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAg                                               >  1:1103657/1‑150 (MQ=255)
                                                                                                   ccATGCGCGGGCTTGCCCGCTCCCCCAGTACTCGCACCGCACTTACGGCCTGGGCAGGACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGc                                              <  2:5187018/150‑1 (MQ=255)
                                                                                                     aTTCGCGGGCATGCCCGCTCCCCCAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGccc                                            <  1:5305150/150‑1 (MQ=255)
                                                                                                     aTTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGccc                                            <  1:1557928/150‑1 (MQ=255)
                                                                                                     aTTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGccc                                            <  2:5013551/150‑1 (MQ=255)
                                                                                                     aTTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACGTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCGAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACGTCCAGCGCACGTACAGTAGccc                                            >  2:1074/1‑150 (MQ=255)
                                                                                                        cgcgGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGc                                         <  1:4570005/150‑1 (MQ=255)
                                                                                                          cgGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCaa                                       >  2:5120098/1‑150 (MQ=255)
                                                                                                          cgGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCaa                                       <  2:4965018/150‑1 (MQ=255)
                                                                                                          cgGGCACGCCCGCTCCCACAGGTACTGCACAGGCCTTAAGTGCGGTGCGATTACTGTGGGAGCGGGCATGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCaa                                       >  1:6739310/1‑150 (MQ=17)
                                                                                                           gggCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAAt                                      >  1:5066331/1‑150 (MQ=255)
                                                                                                           gggCACGCCCGCTCCCACAGGTACTGCACAGGCCTTAAGTGCGGTGCGATTACTGTGGGAGCGGGCATGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAAt                                      >  1:3340870/1‑150 (MQ=17)
                                                                                                            ggcgcgCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCACTg                                     >  1:822025/6‑150 (MQ=255)
                                                                                                              cATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGaa                                   <  2:1854634/150‑1 (MQ=255)
                                                                                                               aggCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAAc                                  <  1:4129934/148‑1 (MQ=255)
                                                                                                               acgCCCGCTCCCACAGGTACTGCACAGGCCTTAAGTGCGGTGCGATTACTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCACTGAAc                                  >  1:6831371/3‑150 (MQ=17)
                                                                                                                 gCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAg                                <  2:2238470/150‑1 (MQ=255)
                                                                                                                 gCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAg                                >  1:1992931/1‑150 (MQ=255)
                                                                                                                 gCCCGCTCCCACAGGTACTGCACAGGCCTTAAGTGCGGTGCGATTACTGTGGGAGCGGGCATGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAg                                >  2:1134971/1‑150 (MQ=17)
                                                                                                                   ccGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGcc                              >  1:1916275/1‑150 (MQ=255)
                                                                                                                   ccGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGcc                              >  2:3953705/1‑150 (MQ=255)
                                                                                                                   ccGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGcc                              >  1:4167962/1‑150 (MQ=255)
                                                                                                                       tCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAg                          >  1:3871034/1‑150 (MQ=255)
                                                                                                                       tCCCACAGTAATCGCACCGCACTTAAGGCCTGAGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCAAGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAg                          >  2:6209740/1‑150 (MQ=255)
                                                                                                                         ccACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATc                                                                                                          <  1:5018858/68‑1 (MQ=255)
                                                                                                                         ccACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATc                                                                                                          <  1:5018663/68‑1 (MQ=255)
                                                                                                                         ccACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATc                                                                                                          >  2:5018858/1‑68 (MQ=255)
                                                                                                                          cacaGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCt                       >  1:1289128/1‑150 (MQ=255)
                                                                                                                          cacaGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCt                       >  1:182388/1‑150 (MQ=255)
                                                                                                                          cacaGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCt                       >  1:1632792/1‑150 (MQ=255)
                                                                                                                          cacaGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCt                       >  2:3269854/1‑150 (MQ=255)
                                                                                                                          cacaGTAATCGCACCGCACTGAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCt                       >  2:2894971/1‑150 (MQ=255)
                                                                                                                           acaGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCACTGAACAGCCCCAGGCTc                      >  2:2066697/1‑150 (MQ=255)
                                                                                                                           acaGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTc                      >  2:2770952/1‑150 (MQ=255)
                                                                                                                           acaGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTc                      >  1:844017/1‑150 (MQ=255)
                                                                                                                              gTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCAcc                   <  1:3694214/150‑1 (MQ=255)
                                                                                                                                      accgcacTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTc           >  1:2718019/1‑150 (MQ=255)
                                                                                                                                       ccgcacTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCa          >  2:2364411/1‑150 (MQ=255)
                                                                                                                                       ccgcacTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCa          >  2:2158063/1‑150 (MQ=255)
                                                                                                                                        cgcacTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCaa         >  2:3800345/1‑150 (MQ=255)
                                                                                                                                        cgcacTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCaa         <  1:2216429/150‑1 (MQ=255)
                                                                                                                                        cgcacTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCaa         >  2:1101100/1‑150 (MQ=255)
                                                                                                                                          cacTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAagc       >  2:5806367/1‑150 (MQ=255)
                                                                                                                                             ttAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAagccag    >  1:563344/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCAGGCGTTCAACCAGTGGCATGCAGAGAGCGGGCGGCCGCTGTGCCAGTTGTCGGTGATGCTGGCGCACACCGTCAACCATTGAGCCTTGTGTTGCTCCCATTCGCGGGCATGCCCGCTCCCACAGTAATCGCACCGCACTTAAGGCCTGTGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAATCTAACTGCCTTCGCCCGCCATGCGCCGCTCATGCTGGAACTTCCAGCGCACGTACAGTAGCCCGGCAATGAACAGCCCCAGGCTCACCAGCACTTCAAGCCAGCC  >  NC_002947/1842960‑1843252

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: