Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,549,653 T→A 9.6% intergenic (+18/‑154) PP_t47 → / → PP_4881 tRNA‑Leu/iron ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,549,6530TA 9.6% 313.8 / 22.1 104intergenic (+18/‑154)PP_t47/PP_4881tRNA‑Leu/iron ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (49/45);  new base A (6/4);  total (55/49)
Fisher's exact test for biased strand distribution p-value = 7.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

Only 100 of 106 total aligned reads displayed.

TGAAAGCGAACGATAGTCGCCTTCTGCAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTAGCCGTAAGCGAGGAAGTAC  >  NC_002947/5549505‑5549802
                                                                                                                                                    |                                                                                                                                                     
tgAAAGCGAACGATAGTCGCCTTCTGCAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGATGGTTTCAGGTATCAGTGACCGCAAGGTCGGGGATGTGCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCtt                                                                                                                                                      >  2:1187997/1‑150 (MQ=255)
   aaGCGAACGATAGTCGCCTTCTGCAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCt                                                                                                                                                   >  2:304620/1‑150 (MQ=255)
   aaGCGAACGATAGTCGCCTTCTGCAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCt                                                                                                                                                   >  2:1291465/1‑150 (MQ=255)
        aaCGATAGTCGCCTTCTGCAGAGGAAGTAAAGCTTTGGCCAGGTGGCGGTATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTgg                                                                                                                                              <  2:77166/150‑1 (MQ=255)
               gTCGCCTTCTGCAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGAAGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGtttttt                                                                                                                                       >  2:67728/1‑150 (MQ=255)
                tCGCCTTCTGCAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAg                                                                                                                                                 >  2:473003/1‑139 (MQ=255)
                tCGCCTTCTGCAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAg                                                                                                                                                 <  1:473003/139‑1 (MQ=255)
                tCGCCTTCTGCAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTc                                                                                                                                      >  1:1930491/1‑150 (MQ=255)
                  gCCTTCTGCAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCtt                                                                                                                                    <  1:57054/150‑1 (MQ=255)
                    cTTCTGCAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGc                                                                                                                                  <  2:1322449/150‑1 (MQ=255)
                    cTTCTGCAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGc                                                                                                                                  <  2:1930491/150‑1 (MQ=255)
                     ttCTGCAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGcc                                                                                                                                 >  1:866893/1‑150 (MQ=255)
                         gCAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGtt                                                                                                                                           <  1:72639/136‑1 (MQ=255)
                         gCAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGtt                                                                                                                                           >  2:72639/1‑136 (MQ=255)
                                 tgTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTAtgt                                                                                                                     <  2:770253/150‑1 (MQ=255)
                                  gTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCGGCTTTGTGTGACCGCAATGTCGTGGAAGTACGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTAtgtg                                                                                                                    <  1:918663/150‑1 (MQ=25)
                                   tAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGc                                                                                                                   >  2:1189155/1‑150 (MQ=255)
                                         tttGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTc                                                                                                             <  1:1167879/150‑1 (MQ=255)
                                              ccAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTCATTGGCGGGATTTGTCTGGCTTGCGATGTATGTGCCTCCTgccccc                                                                                                        >  2:318092/1‑150 (MQ=17)
                                               cAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTcc                                                                                                            <  1:960862/145‑1 (MQ=255)
                                                              gTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTg                                                                                        >  2:281747/1‑150 (MQ=255)
                                                              gTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTg                                                                                        >  1:1196371/1‑150 (MQ=255)
                                                              gTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATAACAGAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCtt                                                                                        >  2:1522629/1‑149 (MQ=255)
                                                               tAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCAAGCTAGTgg                                                                                                                                              >  1:930366/1‑95 (MQ=255)
                                                               tAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCAAGCTAGTgg                                                                                                                                              <  2:930366/95‑1 (MQ=255)
                                                               gaGACGCGCGGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGa                                                                                       <  1:411642/149‑1 (MQ=255)
                                                                aGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGAc                                                                                      <  1:617387/150‑1 (MQ=255)
                                                                aGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCAAGCTAGTGGGttttt                                                                                                                                        <  1:327761/100‑1 (MQ=255)
                                                                aGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCAAGCTAGTGGGttttt                                                                                                                                        >  2:327761/1‑100 (MQ=255)
                                                                   cgcgCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCAt                                                                                   <  2:527667/150‑1 (MQ=255)
                                                                           tttCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTtct                                                                           >  2:447275/1‑150 (MQ=255)
                                                                           tttCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCAAGCTAGTGGGttttt                                                                                                                                        >  2:1566089/1‑89 (MQ=255)
                                                                           tttCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCAAGCTAGTGGGttttt                                                                                                                                        <  1:1566089/89‑1 (MQ=255)
                                                                            ttCAGGTATCAGTGACCGCAAGGTCGTGGACGTTCGTGTCTTCTCCTTGGAACCATACAAAAACCCACTAGCTTGCTGGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTtctc                                                                          <  1:622167/150‑1 (MQ=39)
                                                                            ttCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTtctc                                                                          >  1:1352864/1‑150 (MQ=255)
                                                                                 gTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTcac                                                                     >  2:1240256/1‑150 (MQ=255)
                                                                                  tATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGt                                                                                                                                            <  1:887296/78‑1 (MQ=255)
                                                                                  tATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGt                                                                                                                                            >  2:887296/1‑78 (MQ=255)
                                                                                   aTCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAATCCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCacac                                                                   <  1:809783/150‑1 (MQ=255)
                                                                                       gTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTgg                                                                                                                                              >  1:1134818/1‑71 (MQ=255)
                                                                                       gTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTgg                                                                                                                                              <  2:1134818/71‑1 (MQ=255)
                                                                                       gTGACCGCAAGGTCGTGGAAGTTCGAGTCTTATCCTGTGCACCATACAAAAACCCACTAGCTTGCTCGTGGGTTTTTTCTTGCCTGCAATTTATGTGCCCCCTCCCCCCCTTGCACCATTGTCTGCCCATTAATTTCTCCTCACACGGcc                                                               >  2:680488/1‑150 (MQ=38)
                                                                                            cGCAAGGTCGTGGAAGTTGGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTg                                                          >  2:573279/1‑150 (MQ=255)
                                                                                            cGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGtt                                                                                                                                           <  1:628510/69‑1 (MQ=255)
                                                                                            cGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGtt                                                                                                                                           >  2:628510/1‑69 (MQ=255)
                                                                                            cGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTg                                                          >  2:1359914/1‑150 (MQ=255)
                                                                                             gCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTga                                                         <  1:1189155/150‑1 (MQ=255)
                                                                                                  gTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCtt                                                                                                     <  2:946780/101‑1 (MQ=255)
                                                                                                     gTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAAt                                                 >  1:655431/1‑150 (MQ=255)
                                                                                                      tGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAAtt                                                >  2:1646184/1‑150 (MQ=255)
                                                                                                       ggAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTGTCGTCACACGGCCACTTGAGAAACAAttt                                               >  1:713646/1‑150 (MQ=255)
                                                                                                        gAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTCGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTc                                              <  1:1187997/150‑1 (MQ=255)
                                                                                                        gAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTc                                              <  2:1580422/150‑1 (MQ=255)
                                                                                                        gAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTc                                              <  1:1359914/150‑1 (MQ=255)
                                                                                                          aGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGc                                                                                                                                  >  2:1160541/1‑64 (MQ=255)
                                                                                                          aGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGc                                                                                                                                  <  1:1160541/64‑1 (MQ=255)
                                                                                                           gTTCGTGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCCGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGtt                                           <  1:318092/150‑1 (MQ=255)
                                                                                                             tCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTa                                         <  1:1646184/150‑1 (MQ=255)
                                                                                                               gAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTAcc                                       >  1:1317828/1‑150 (MQ=255)
                                                                                                                aGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCa                                      <  2:1196371/150‑1 (MQ=255)
                                                                                                                  tcttctCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCCCCCGGCCACTTGAGAAACAATTTCGTTTACCAtt                                    >  2:743555/1‑150 (MQ=255)
                                                                                                                  tcttctCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACg                                                                  <  2:1757315/120‑1 (MQ=255)
                                                                                                                  tcttctCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACg                                                                  >  1:1757315/1‑120 (MQ=255)
                                                                                                                  tcttctCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGtt                                           >  2:9110/1‑143 (MQ=255)
                                                                                                                  tcttctCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGtt                                           <  1:9110/143‑1 (MQ=255)
                                                                                                                  tcttctCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCAtt                                    >  2:1132925/1‑150 (MQ=255)
                                                                                                                     tctcCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGtt                                 <  1:447275/150‑1 (MQ=255)
                                                                                                                      ctcCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGttt                                <  2:655431/150‑1 (MQ=255)
                                                                                                                       tcCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGttag                               <  1:67728/150‑3 (MQ=255)
                                                                                                                         cTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCa                             >  1:1641956/1‑150 (MQ=255)
                                                                                                                          tGGGCACCATACAAAAACACACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCAATTGAGAAACAATTTCGTTTACCATTGTTTCCaa                            >  2:916340/1‑150 (MQ=255)
                                                                                                                             gCACCATACAAAAACGCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAata                         <  2:1352864/150‑1 (MQ=255)
                                                                                                                               aCCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATg                       >  2:1403032/1‑150 (MQ=255)
                                                                                                                               aCCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATg                       >  1:325495/1‑150 (MQ=255)
                                                                                                                                 cATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTa                     <  1:304620/150‑1 (MQ=255)
                                                                                                                                 cATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTa                     <  1:1132925/150‑1 (MQ=255)
                                                                                                                                  aTACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATgg                                                                                            >  1:529631/1‑78 (MQ=255)
                                                                                                                                  aTACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATgg                                                                                            <  2:529631/78‑1 (MQ=255)
                                                                                                                                  aTACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTAg                    >  1:1271000/1‑150 (MQ=255)
                                                                                                                                   agaaaaaaCCCACTAGCAAGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGTCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTAGc                   >  2:546012/4‑150 (MQ=255)
                                                                                                                                    aCAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTccccccc                                                                                                       <  1:34796/65‑1 (MQ=255)
                                                                                                                                    aCAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTccccccc                                                                                                       >  2:34796/1‑65 (MQ=255)
                                                                                                                                      aaaaaCCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCg                                                                         <  2:284458/93‑1 (MQ=255)
                                                                                                                                      aaaaaCCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCg                                                                         >  1:284458/1‑93 (MQ=255)
                                                                                                                                      aaaaaCCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGtt                                           >  1:826005/1‑123 (MQ=255)
                                                                                                                                      aaaaaCCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGtt                                           <  2:826005/123‑1 (MQ=255)
                                                                                                                                       aaaaCCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAata                         <  2:1721857/140‑1 (MQ=255)
                                                                                                                                       aaaaCCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAata                         >  1:1721857/1‑140 (MQ=255)
                                                                                                                                       aaaaCCCACTAGCAAGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTAg                    <  1:331459/145‑1 (MQ=255)
                                                                                                                                       aaaaCCCACTAGCAAGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTAg                    >  2:331459/1‑145 (MQ=255)
                                                                                                                                        aaaCCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTAGCCGTaa              <  1:1291465/150‑1 (MQ=255)
                                                                                                                                        aaaCCCACTAGCAAGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTAGCCGTaa              >  1:47309/1‑150 (MQ=255)
                                                                                                                                         aaCCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACAAGGCCACTTGATAAACAATTTCGTTTACCATTGTTTCCAAATATGTAGCCGTAAg             >  2:95084/1‑150 (MQ=255)
                                                                                                                                          aCCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTAGCCGTAAGc            >  1:186460/1‑150 (MQ=255)
                                                                                                                                          aCCCACTAGCTTGCTAGTGCGCTTTTTCTTGCGGGCGATTTATGTGCCCCCTCCCCCCCTTGCACAATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCATTTACCATTGTTTCCAAATATGTAGCCGTAAGc            >  2:716015/1‑150 (MQ=255)
                                                                                                                                            ccACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTAGCCGTAAGCGa          >  1:220770/1‑150 (MQ=255)
                                                                                                                                            ccACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTAGCCGTAAGCGa          <  1:1240256/150‑1 (MQ=255)
                                                                                                                                                tAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTAGCCGTAAGCGAGGaa      <  1:281747/150‑1 (MQ=255)
                                                                                                                                                  gCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGt                                            <  1:1642106/110‑1 (MQ=255)
                                                                                                                                                    ttGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACAAGGCCACTTGATAAACAATTTCGTTTACCATTGTTTCCAAATATGTAGCCGTAAGCGAGGAAGTAc  <  1:95084/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGAAAGCGAACGATAGTCGCCTTCTGCAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTAGCCGTAAGCGAGGAAGTAC  >  NC_002947/5549505‑5549802

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