Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 668,110 T→A 11.7% intergenic (+71/+81) yfhM → / ← PP_0571 lipoprotein/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947668,1100TA11.7% 417.1 / 5.8 155intergenic (+71/+81)yfhM/PP_0571lipoprotein/membrane protein
Reads supporting (aligned to +/- strand):  ref base T (63/67);  new base A (19/3);  total (84/71)
Fisher's exact test for biased strand distribution p-value = 9.48e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 156 total aligned reads displayed.

GCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCG  >  NC_002947/667961‑668255
                                                                                                                                                     |                                                                                                                                                 
gcAGGTCGATTCGATGTACGGGCCGAACCTGCTGGCCGTGGGCGATGGGCATGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGTGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAAt                                                                                                                                                   <  2:2582339/150‑1 (MQ=39)
      cGAATCGATGTACCGGCCGAACCTGCAGGCGGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATCGGGCCGCAAAGCGGCCCCGACAATTCAAgc                                                                                                                                             <  1:377049/150‑1 (MQ=255)
       gAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCCACAAATCAAgcg                                                                                                                                            >  2:2052144/1‑150 (MQ=255)
       gAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCGATATGGCCGCCAAGCGGCCCCGCCCATTCAAgcg                                                                                                                                            >  1:4114361/1‑150 (MQ=255)
        aaTCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAgcgg                                                                                                                                           >  2:4356142/1‑150 (MQ=255)
        aaTCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAgcgg                                                                                                                                           <  2:5384063/150‑1 (MQ=255)
        aaTCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGGGGCGATAGGGCCGCCAAGCGGGCCCGCCCAATCAAgcgg                                                                                                                                           >  1:150530/1‑150 (MQ=255)
            gATGTACCGGCCGAACCTGCAGGCGGTGGGCGTTGGGCACGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGaa                                                                                                                                       <  2:2013388/150‑1 (MQ=255)
             aTGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAg                                                                                                                                      <  2:1263274/150‑1 (MQ=255)
                tACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCAAAGCGGCCCCCACAATTAAAGCGGCGAAGCTg                                                                                                                                   >  2:3002802/1‑150 (MQ=255)
                  ccGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAACGCGGCCCCGACAATTCAAGCGGCGAAGCTGaa                                                                                                                                 <  2:5468874/150‑1 (MQ=255)
                  ccGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGCCAAATCAAGCGGCGAAGCTGaa                                                                                                                                 >  2:508111/1‑150 (MQ=255)
                   cGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGaaa                                                                                                                                <  2:3928585/150‑1 (MQ=255)
                     gCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCCGCGCGCTAGGGCCGCCAAGCGGCCCCGCCCAATCAAGCGGCGCAGCGGAAAGc                                                                                                                              >  2:4425715/1‑150 (MQ=255)
                       cGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCt                                                                                                                            <  2:2733046/150‑1 (MQ=255)
                        gAACCTGCAGGCCGTGGGCGATGGGCATGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAATGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTg                                                                                                                           <  1:1280900/150‑1 (MQ=255)
                        gAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGTAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTg                                                                                                                           <  2:4009135/150‑1 (MQ=255)
                        gAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAgcggc                                                                                                                                          >  1:5354457/1‑135 (MQ=255)
                        gAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTg                                                                                                                           <  2:605854/150‑1 (MQ=255)
                        gAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCGGAAAGCCTg                                                                                                                           >  2:1309580/1‑150 (MQ=255)
                         ccccTGCCGGCCGGGGGCGCTGGGCACGGGGAGAGGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGGGGGCGCGGCCTTGCGTCGCGATCGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTgg                                                                                                                          <  2:1502895/148‑1 (MQ=255)
                         acccTGCAGGCCGTGGGCGTTGGGCATGGGGAGATGACGGTTAAAGCCCGCTAGAGACCCGCTCTACGCGTAATGTGTGGGAGCGGCCTTGCGTCGCTATTGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTgg                                                                                                                          <  1:4858314/148‑1 (MQ=38)
                         aaCCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTgg                                                                                                                          <  1:2232763/150‑1 (MQ=255)
                          aCCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGCGTGCGTCGCGATAGGGCCCCAAAGCGGCCCCGACAAATAAACCCGCGAACCAGAAAGACAggg                                                                                                                         >  1:5266940/1‑150 (MQ=255)
                          aCCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTggg                                                                                                                         <  2:2736370/150‑1 (MQ=255)
                            cTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCCATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCCAAGCGGAAAGCCAGGGgc                                                                                                                       >  2:4821848/1‑150 (MQ=39)
                            cTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGGCGCGATAGGGCCGCACAGCGGCCCCGACCATTCAAGCGGCGAAGCTGAAAGCCTggggg                                                                                                                       >  2:3904106/1‑149 (MQ=255)
                             tGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGGGGGCGGGGCGCTGGGGCCGCACGGCGGGCCCCCCCCATCAAGCGGCGCCGCGGGAGGCCGGGGgcg                                                                                                                      >  1:5425327/1‑149 (MQ=255)
                                aGGGCGTGGGCGATGGGCAAGGGTTGTTGTGGGTGAAAGCCCGCTTGATCCCCGGTATACGGTGAATGTGTGGGTGCGGGCTTGGGTCGCGTTGGGGCCGCAAAGGGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGgctgct                                                                                                                   <  1:397866/150‑1 (MQ=255)
                                aGGCCGTGGGCGATGGGCATGGGTAGATGACGGTTAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCTATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGgctgct                                                                                                                   <  1:2606763/150‑1 (MQ=255)
                                   ccGTGGGGTTTGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAACGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCtgtg                                                                                                                <  1:3479100/150‑1 (MQ=255)
                                   ccGTGGGCTATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCtgtg                                                                                                                <  1:5569447/150‑1 (MQ=255)
                                   ccGTGGGCGTTGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGGGGAATGTGTGGGAGCGGCCTTGCGCCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCtgtg                                                                                                                <  1:3870283/150‑1 (MQ=255)
                                     gTGGGCTATGGGCAAGGGTAGATGACGTTGAAAGCCCGCTAGATACCCGCTCTACGCGTAAAGTTTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCa                                                                                                              <  1:2233514/150‑1 (MQ=38)
                                     gTGGGCGTTGGGCAAGGGGATATGACGGTGAACGCCCGCTAGAGACCCGCTCTCCGGGGAAAGTGTGGGCGCGGCCTTGGGTCGCGATAGGGCCGCAACGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCa                                                                                                              <  1:1046523/150‑1 (MQ=255)
                                     gTGGGCGTGGGGCATGGGGAGATGCGGGTTAAAGCCCGCTAGAGACCCGCTCTAGGCGGAAAGTGTGGGCGCGGCCTGGCGTCGCTATCGGGCCGCAACGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCa                                                                                                              <  1:477941/150‑1 (MQ=255)
                                     gTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGTAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCa                                                                                                              <  1:3165011/150‑1 (MQ=255)
                                     gTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCa                                                                                                              <  2:4358450/150‑1 (MQ=255)
                                     gTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCa                                                                                                              >  2:3304674/1‑150 (MQ=255)
                                          cgttGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCt                                                                                                         <  2:4166679/147‑1 (MQ=255)
                                          cGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGGTTGGGGGGGGGTTGGGGCGGGTAGGGGGGCCGGCAATTCAAGCGGGGAAGCTGATAGGCTGGGGGTGGTGTGGGGCCCt                                                                                                         >  1:4671530/1‑150 (MQ=255)
                                          cGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCt                                                                                                         <  2:1608281/150‑1 (MQ=255)
                                          cGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGGCGCGCTAGGGCCGCCAAGCGGCCCCGACAAATCAAGCGGCGAAGCGGCAAGCCTGGGGCTGCTGTGCAGcccg                                                                                                         >  2:4847493/1‑149 (MQ=255)
                                           gATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCCAAGCGGCCCCGACAAATCAAGCGGCGCAGCTGAAAACCAGGGGCAGCGGTGCAGcccag                                                                                                        >  1:4617969/1‑148 (MQ=38)
                                              gggCATGGGGAGATGACGGTGAAAGCCCGCTAGAGCCCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCACAGCAGccccaggc                                                                                                     <  2:4329782/150‑6 (MQ=37)
                                                gCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGGGGGTGGCGATAGGGCCGCAAAGCGGGCCCGACACATCAAGCGGCGCAGCGGAAAGCCTGGGGCTGCTGTGCGGCCCAGtcgcg                                                                                                   <  1:5040422‑M2/150‑6 (MQ=255)
                                                gCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCATGCGTCGCGCTAGGGCCCCAAATCGGCCCCCACAAAAAAATCGGCGAAGCAGAAAGCCAGGGGCGGCTGTGCAGCCCAGgcgcg                                                                                                   <  2:1885696‑M2/150‑6 (MQ=255)
                                                 catggggAGATGACGGTGAACGCCCGCTAGAGACCCGCTCTACGCGGAACGTGTGGGAGCGGCCTTGCGTCGCGATCGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGa                                                                                                  <  1:1609651/147‑1 (MQ=255)
                                                 cAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcga                                                                                                  >  2:2193562‑M2/1‑144 (MQ=255)
                                                    ggggAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgc                                                                                               <  2:5531117‑M2/150‑10 (MQ=255)
                                                      ggtgaTGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCaa                                                                                             <  1:4522785/147‑1 (MQ=255)
                                                        agaTGACGGTGAAAGCCCGCTAGAGACCCTCTCTACGCGGAATGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgactcaagg                                                                                           >  2:4561160‑M2/1‑137 (MQ=255)
                                                        agaTGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAACGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaagg                                                                                           >  1:2171185‑M2/1‑137 (MQ=255)
                                                         gaTGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaaggc                                                                                          >  2:2781507‑M2/1‑136 (MQ=255)
                                                                     cccccGCTCGACACCCGCTCTACGGGTAATGTGGGGGGGCGGCGTTGGGACGCGATGGGGCCGCAACGCGCCCCCGACAATGCCCGCGGCGAAGGTTAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCACGGCCGCTCCCACaaa                                                                              <  2:3625191/148‑1 (MQ=255)
                                                                        ccGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAAATCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCACGGCCGCTCCCACAAAGAc                                                                           >  2:3621885/1‑150 (MQ=255)
                                                                         cGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCCATAGGGCCCCAAAGCGGCCCCCACAAAACAAACGGCCGAGCTGCAACCCGGGGGGCGCTGTGCCCCCCGGTCGCGCCGCAACGCCGCTCCCCCAAAGAcc                                                                          >  2:747575/1‑150 (MQ=255)
                                                                          gCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCCAAGCCGCCCCGACAAATCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCGGTCGCGACGCCAGGCCGCGCCCCCAAAGACcg                                                                         >  2:148562/1‑150 (MQ=37)
                                                                          gCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCa                                                                                                              >  1:1643540/1‑113 (MQ=255)
                                                                          gCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCa                                                                                                              <  2:1643540/113‑1 (MQ=255)
                                                                                         tACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAAATCAAGAGGCGAAGCAGAAACCCAGGGGCTGCTGTGCAGCCCAGTCGCGACGCAAGGCCGCTCCCACAAATACCGCGCAGATCAGTGGAt                                                          >  2:4206578/1‑150 (MQ=25)
                                                                                         tACGCGGAAAGTGTGGGAGCGGCCGTTCGTGGCGGTATGGCCCCCAAGCGGCCCCCACCAATCAACCGGCGAAGCTGTAAGCCTGGGGCTGGGGGGCAGCCCTGGCGCGACGCAAGGCCGCTCCCCCCAAGCACGCGCCGCTAAGTGGAt                                                          >  2:5172560/1‑150 (MQ=255)
                                                                                         tACGCGGAAAGTGTGGGAGCGGCCGTGCCTCGCGATAGGGCCGCAAAACGGCCCCGACAAAGCAAGCGGCGCAGCGGGAAGCCGGGGGCGGCGGGGCAGCCCGGGGGCGACGCAACGCCGCGCCCCCCAAGAACCCCCAGAGCAGGGGAt                                                          >  2:577723/1‑150 (MQ=255)
                                                                                         tACGCGGAAAGTGTGGGAGCGGCCATGCCTCGCGATAGGGCCGCCAAGCGGCCCCGACAAATAAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGGGCAGCCCAGTCGCGACGCAAGGCCGCTCCCCCCAAGACCGCGCCCATAACTGGAt                                                          >  2:5172604/1‑150 (MQ=25)
                                                                                            gcgGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCAc                                                       >  2:4632490/1‑150 (MQ=255)
                                                                                              ggAAAGTGTGGGCGCGGCGTTGCGTCGCGATAGGGCCGCAACGCGGCCCCGACAATTCCAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGATGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCAccc                                                     <  2:5097144/150‑1 (MQ=37)
                                                                                                        ggAGCGGCCGTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTAAAGCGGCGAACCTGAAGGCCAGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACCAAGACCGCGCAGATCAGTGGATAACCCACCTCAacac                                           >  2:1228058/1‑150 (MQ=37)
                                                                                                              gCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAg                                     >  2:4414966/1‑150 (MQ=255)
                                                                                                              gCCTTGCGTCGCGATACGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGTCCTCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCAACCACAg                                     >  2:3841560/1‑150 (MQ=255)
                                                                                                               ccTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGc                                    >  2:2310893/1‑150 (MQ=255)
                                                                                                                 ttGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGTGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGccc                                  >  1:8139/1‑150 (MQ=255)
                                                                                                                     gTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatctcccctt                              <  1:4408169‑M2/150‑75 (MQ=255)
                                                                                                                      tCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCa                             >  1:5606731/1‑150 (MQ=255)
                                                                                                                      tCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCa                             >  1:2797237/1‑150 (MQ=255)
                                                                                                                       ccgtttcgggCCGAAACGCGCCCCCCAAAATACATGCGGCGAAGCTGAACGCCGGGGGCTGCTGAGCCGCCCGGCCCCCACCCACGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCAc                            <  1:5568022/143‑1 (MQ=255)
                                                                                                                         cgATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatctccccttgccc                          <  1:930167‑M2/150‑79 (MQ=255)
                                                                                                                         cgATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCAcca                          >  2:1224046/1‑150 (MQ=255)
                                                                                                                          gATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatctccccttgccca                         <  1:300473‑M2/150‑80 (MQ=255)
                                                                                                                          gATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCAccag                         >  2:1837656/1‑150 (MQ=255)
                                                                                                                           aTAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatctccccttgcccat                        <  1:2834861‑M2/150‑81 (MQ=255)
                                                                                                                           aTAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCAccagc                        >  1:3648305/1‑150 (MQ=255)
                                                                                                                           aTAGGGCCGCAAAACGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcagctccccgtgcccag                        <  1:2834835‑M2/150‑81 (MQ=255)
                                                                                                                             ggggCCGCTTTGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGAGCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCAccagcca                      >  2:3117263/2‑150 (MQ=255)
                                                                                                                              gggCCGCTTTGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCAccagccag                     >  1:3224188/1‑150 (MQ=255)
                                                                                                                              gggCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCAccagccag                     >  1:986668/1‑150 (MQ=255)
                                                                                                                               ggCCGCAACGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGCCGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAga                    <  2:3177926/150‑1 (MQ=255)
                                                                                                                               ggCCGCAAAGCGGCCCCGACAATTCAAGGGGCTAAGATGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCACGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAga                    <  2:3177994/150‑1 (MQ=255)
                                                                                                                               ggCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGAc                                                                           <  1:4987189/95‑1 (MQ=255)
                                                                                                                               ggCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGAc                                                                           >  2:4987189/1‑95 (MQ=255)
                                                                                                                               ggCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAga                    >  1:1022238/1‑150 (MQ=255)
                                                                                                                                  cGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTAtcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatctccccttgcccatcgcccac                 <  1:1556057‑M2/150‑88 (MQ=255)
                                                                                                                                   gCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAgatgat                >  2:5302724/1‑150 (MQ=255)
                                                                                                                                    cAAATCGGCCCAGACAATTCAAGCGGCTAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTTTCGCGAAGCAAGGACGCTCCAACAAATACCGCTCAGATCAGTGGATAACCCAGCTCATCACAAACAGCACCAGCACCAGCCAGATGAtt               >  1:1203352/1‑150 (MQ=255)
                                                                                                                                       aGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCg            >  1:4150568/1‑150 (MQ=255)
                                                                                                                                           gCCCCTACAATTCAAGCGGCGACGCTGAACGCGTGGGGGTGGTGTGCCGCCGTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGga        <  2:5425327/150‑1 (MQ=255)
                                                                                                                                              ccGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGgatga     >  2:4958227/1‑150 (MQ=255)
                                                                                                                                              ccGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGgatga     >  2:1662244/1‑150 (MQ=255)
                                                                                                                                                gACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCATCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCAACCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATc   >  2:448165/1‑150 (MQ=255)
                                                                                                                                                gACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGt                                                              >  1:633463/1‑91 (MQ=255)
                                                                                                                                                gACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATc   >  2:2928783/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCG  >  NC_002947/667961‑668255

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