Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,520,941 C→T 7.3% intergenic (+14/+22) mcbR → / ← yncD DNA‑binding transcriptional dual regulator McbR/putative TonB‑dependent outer membrane receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,520,9410CT 7.3% 331.0 / 10.9 110intergenic (+14/+22)mcbR/yncDDNA‑binding transcriptional dual regulator McbR/putative TonB‑dependent outer membrane receptor
Reads supporting (aligned to +/- strand):  ref base C (54/46);  new base T (6/2);  total (62/48)
Fisher's exact test for biased strand distribution p-value = 2.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.40e-01

Only 100 of 110 total aligned reads displayed.

TTACGGGAGCATATAGAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATT  >  NC_000913/1520793‑1521085
                                                                                                                                                    |                                                                                                                                                
ttACGGGAGCATATAGAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCa                                                                                                                                                 <  1:969156/150‑1 (MQ=255)
ttACGGGAGCATATAGAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCa                                                                                                                                                 <  2:1488660/150‑1 (MQ=255)
     ggAGCATATAGAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATAc                                                                                                                                            >  2:600463/1‑150 (MQ=255)
               gAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGTATTTGCCGGGCAttta                                                                                                                                  <  2:1969604/150‑4 (MQ=255)
               gAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCaa                                                                                                                                  >  2:1125113/1‑150 (MQ=255)
               gAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCaa                                                                                                                                  >  1:1970661/1‑150 (MQ=255)
                aaaaCTATCATCTGTTACTCGCTGCGTTAAAACCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGTATTTGCCGGGCAtttac                                                                                                                                 <  2:1471863/150‑5 (MQ=255)
                   aCTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTa                                                                                                                              >  1:810978/1‑150 (MQ=255)
                    cTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTaa                                                                                                                             >  1:1608486/1‑150 (MQ=255)
                      atcatcTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAAtt                                                                                                                           <  2:217323/150‑1 (MQ=255)
                      atcatcTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATAGTTAAAGTAAATGCCCGTCAATTACCGGTCATGAAGGTAAtt                                                                                                                           <  2:1262344/150‑1 (MQ=255)
                        catcTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTAc                                                                                                                         >  1:848986/1‑150 (MQ=255)
                         atcTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACt                                                                                                                        <  2:1970661/150‑1 (MQ=255)
                           cTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCa                                                                                                                      >  2:290215/1‑150 (MQ=255)
                           cTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCa                                                                                                                      <  1:1689851/150‑1 (MQ=255)
                            tGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATGCCGGGCATCAAGGTAATTACTCaa                                                                                                                     >  2:134595/1‑150 (MQ=255)
                                 cTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAAtctc                                                                                                                >  1:990146/1‑150 (MQ=255)
                                   cGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCa                                                                                                              >  2:1992141/1‑150 (MQ=255)
                                   cGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCa                                                                                                              >  1:652092/1‑150 (MQ=255)
                                   cGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCa                                                                                                              >  2:1103571/1‑150 (MQ=255)
                                         gTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAat                                                                                                        <  1:1125113/150‑1 (MQ=255)
                                           tAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAatat                                                                                                      >  1:685026/1‑150 (MQ=255)
                                           tAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATAACTCAAATCTCCACGCAatat                                                                                                      >  1:1373174/1‑150 (MQ=255)
                                           tAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGAAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAatat                                                                                                      >  1:1843688/1‑150 (MQ=255)
                                            aaaaACAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATAtt                                                                                                     >  2:1107800/1‑150 (MQ=255)
                                               aGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCAt                                                                                                  >  1:153369/1‑150 (MQ=255)
                                               aGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCAt                                                                                                  >  1:147330/1‑150 (MQ=255)
                                               aGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCAt                                                                                                  >  1:570049/1‑150 (MQ=255)
                                                 cAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGc                                                                                                >  2:1131165/1‑150 (MQ=255)
                                                 cAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGc                                                                                                >  2:1635826/1‑150 (MQ=255)
                                                  aaaaGACAAAGAGGGATGCGGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGcc                                                                                               >  2:1076300/1‑150 (MQ=255)
                                                       aCAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAAAAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGacac                                                                                          <  2:175925/150‑1 (MQ=255)
                                                          aaGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCAt                                                                                       <  1:1635826/150‑1 (MQ=255)
                                                          aaGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCAt                                                                                       >  2:1185702/1‑150 (MQ=255)
                                                               ggATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAAttt                                                                                  <  2:810978/150‑1 (MQ=255)
                                                                 aTGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATGCCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCg                                                                                <  1:134595/150‑1 (MQ=255)
                                                                  tGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGt                                                                               >  2:1135893/1‑150 (MQ=255)
                                                                    cAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTcc                                                                             <  2:153369/150‑1 (MQ=255)
                                                                    cAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTcc                                                                             >  2:1087783/1‑150 (MQ=255)
                                                                    cAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTcc                                                                             <  2:570049/150‑1 (MQ=255)
                                                                      taCATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCgg                                                                           <  1:1009330/149‑1 (MQ=255)
                                                                        cATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGc                                                                         <  1:1707585/150‑1 (MQ=255)
                                                                        cATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGc                                                                         >  1:1382904/1‑150 (MQ=255)
                                                                           tGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGaa                                                                      <  1:1103571/150‑1 (MQ=255)
                                                                            gTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAg                                                                     >  2:1621042/1‑150 (MQ=255)
                                                                             tCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAgg                                                                    <  2:1382904/150‑1 (MQ=255)
                                                                                  cTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGt                                                               >  1:339437/1‑150 (MQ=255)
                                                                                  cTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGt                                                               <  2:268187/150‑1 (MQ=255)
                                                                                   tGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTa                                                              <  1:641684/150‑1 (MQ=255)
                                                                                        ttATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATc                                                         >  2:1712052/1‑150 (MQ=255)
                                                                                        ttATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATc                                                         >  2:885942/1‑150 (MQ=255)
                                                                                          aTGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGc                                                       >  1:22023/1‑150 (MQ=255)
                                                                                             caacaaAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAACGTTCGTAATATCGCCCg                                                    >  2:1309327/1‑150 (MQ=255)
                                                                                               acaaAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGtt                                                  <  1:600463/150‑1 (MQ=255)
                                                                                               acaaAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGtt                                                  >  2:820111/1‑150 (MQ=255)
                                                                                               acaaAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGtt                                                  >  2:908282/1‑150 (MQ=255)
                                                                                                caaAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGttt                                                 >  2:1048937/1‑150 (MQ=255)
                                                                                                caaAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGttt                                                 <  2:848986/150‑1 (MQ=255)
                                                                                                   aaTATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGAc                                              <  2:990146/150‑1 (MQ=255)
                                                                                                    atatTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACt                                             <  1:1131165/150‑1 (MQ=255)
                                                                                                     tatTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACta                                            >  1:906587/1‑149 (MQ=255)
                                                                                                         gCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTg                                        <  1:290215/150‑1 (MQ=255)
                                                                                                          cTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGa                                       >  2:525482/1‑150 (MQ=255)
                                                                                                            aTTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACa                                     >  1:45190/1‑150 (MQ=255)
                                                                                                             ttttATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACaa                                    <  1:1135893/150‑1 (MQ=255)
                                                                                                              tttATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAAttt                                                                                  >  1:1987445/1‑103 (MQ=255)
                                                                                                               ttATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATg                                  <  1:1076300/150‑1 (MQ=255)
                                                                                                               ttATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATg                                  <  2:156400/150‑1 (MQ=255)
                                                                                                                      aGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAAc                           <  1:1185702/150‑1 (MQ=255)
                                                                                                                       gCAATACAATCGTTAAAGTAAATGCCCGTCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAAcc                          <  2:1115673/150‑1 (MQ=255)
                                                                                                                       gCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAAcc                          >  1:1043691/1‑150 (MQ=255)
                                                                                                                       gCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAAcc                          <  1:1350051/150‑1 (MQ=255)
                                                                                                                        cAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCa                         <  1:398452/150‑1 (MQ=255)
                                                                                                                           tACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACg                      <  2:1728024/150‑1 (MQ=255)
                                                                                                                           tACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACg                      <  2:147330/150‑1 (MQ=255)
                                                                                                                           tACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACg                      >  1:1515924/1‑150 (MQ=255)
                                                                                                                           tACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACg                      <  2:652092/150‑1 (MQ=255)
                                                                                                                             cAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGAc                                              <  1:1209456/124‑1 (MQ=255)
                                                                                                                             cAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAAc                       >  1:415312/1‑147 (MQ=255)
                                                                                                                             cAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAAc                       <  2:415312/147‑1 (MQ=255)
                                                                                                                                tCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTc                 <  1:1107800/150‑1 (MQ=255)
                                                                                                                                 cGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCt                >  1:49064/1‑150 (MQ=255)
                                                                                                                                  gTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAata                                                           <  1:273101/106‑1 (MQ=255)
                                                                                                                                  gTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAata                                                           >  2:273101/1‑106 (MQ=255)
                                                                                                                                  gTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCtt               >  1:1904749/1‑150 (MQ=255)
                                                                                                                                  gTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGGAATATCGCCCGTTTGACTCATTGACAATGACAGAACCACCGTATTCtt               >  1:1897539/1‑150 (MQ=255)
                                                                                                                                     taccttgaTGCCCGGTATTTGCCGGGCATCAAGGTAATTACTCAAATCCCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTAt            >  2:1276323/8‑150 (MQ=255)
                                                                                                                                     taccttgaTGCCCGGTATTTGCCGGGCATCAAGGTAATTACTCAAATCCCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTAt            >  2:1283097/8‑150 (MQ=255)
                                                                                                                                     aaaGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTAt            <  2:45190/150‑1 (MQ=255)
                                                                                                                                      aaGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATc           >  2:289502/1‑150 (MQ=255)
                                                                                                                                      aaGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATc           <  2:1608486/150‑1 (MQ=255)
                                                                                                                                       ccttgaTGCCCGGTATTTGCCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCg          >  1:730609/6‑150 (MQ=255)
                                                                                                                                       aGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCGAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCg          <  1:1389558/150‑1 (MQ=255)
                                                                                                                                       aGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCg          >  2:1397285/1‑150 (MQ=255)
                                                                                                                                        cttgaTGCCCGGTATTTGCCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCTTAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGa         >  1:1563931/5‑150 (MQ=255)
                                                                                                                                         ttgaTGCCCGGTATTTGCCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGaa        >  2:777876/4‑150 (MQ=255)
                                                                                                                                          aaaTGCCCGGCAAATACCGGGCATCAAGTTAATTACTCAAATCTCCACTCAATATTCATTCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAAt       <  1:208455/150‑1 (MQ=255)
                                                                                                                                            aTGCCCGGTATTTGCCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATaa     >  1:801347/1‑150 (MQ=255)
                                                                                                                                            aTGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATaa     >  2:371139/1‑150 (MQ=255)
                                                                                                                                               cccGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAAtt  <  1:525482/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TTACGGGAGCATATAGAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATT  >  NC_000913/1520793‑1521085

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