Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,520,943 A→T 7.4% intergenic (+16/+20) 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,9430AT 7.4% 341.6 / 12.5 108intergenic (+16/+20)mcbR/yncDDNA‑binding transcriptional dual regulator McbR/putative TonB‑dependent outer membrane receptor
Reads supporting (aligned to +/- strand):  ref base A (56/44);  new base T (6/2);  total (62/46)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 110 total aligned reads displayed.

GGAGCATATAGAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACA  >  NC_000913/1520798‑1521092
                                                                                                                                                 |                                                                                                                                                     
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)
                 atcatcTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAAtt                                                                                                                                  <  1:612613/150‑1 (MQ=255)
                 atcatcTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAAtt                                                                                                                                  <  2:217323/150‑1 (MQ=255)
                 atcatcTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATAGTTAAAGTAAATGCCCGTCAATTACCGGTCATGAAGGTAAtt                                                                                                                                  <  2:1262344/150‑1 (MQ=255)
                   catcTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTAc                                                                                                                                >  1:848986/1‑150 (MQ=255)
                   catcTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTAc                                                                                                                                <  1:1626205/150‑1 (MQ=255)
                    atcTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACt                                                                                                                               <  2:1970661/150‑1 (MQ=255)
                      cTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCa                                                                                                                             <  1:1689851/150‑1 (MQ=255)
                      cTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCa                                                                                                                             >  2:290215/1‑150 (MQ=255)
                       tGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATGCCGGGCATCAAGGTAATTACTCaa                                                                                                                            >  2:134595/1‑150 (MQ=255)
                            cTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAAtctc                                                                                                                       >  1:990146/1‑150 (MQ=255)
                              cGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCa                                                                                                                     >  1:652092/1‑150 (MQ=255)
                              cGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCa                                                                                                                     >  2:1103571/1‑150 (MQ=255)
                              cGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCa                                                                                                                     >  2:1992141/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:147330/1‑150 (MQ=255)
                                          aGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCAt                                                                                                         >  1:570049/1‑150 (MQ=255)
                                          aGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCAt                                                                                                         >  1:153369/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                                                                                              >  2:1185702/1‑150 (MQ=255)
                                                     aaGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCAt                                                                                              <  1:1635826/150‑1 (MQ=255)
                                                          ggATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAAttt                                                                                         <  2:810978/150‑1 (MQ=255)
                                                           gATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTc                                                                                        >  1:232326/1‑150 (MQ=255)
                                                            aTGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATGCCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCg                                                                                       <  1:134595/150‑1 (MQ=255)
                                                             tGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGt                                                                                      >  2:1135893/1‑150 (MQ=255)
                                                               cAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTcc                                                                                    >  2:1087783/1‑150 (MQ=255)
                                                               cAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTcc                                                                                    <  2:153369/150‑1 (MQ=255)
                                                               cAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTcc                                                                                    <  2:570049/150‑1 (MQ=255)
                                                                 taCATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCgg                                                                                  <  1:1009330/149‑1 (MQ=255)
                                                                   cATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGc                                                                                >  1:1382904/1‑150 (MQ=255)
                                                                   cATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGc                                                                                <  1:1707585/150‑1 (MQ=255)
                                                                      tGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGaa                                                                             <  1:1103571/150‑1 (MQ=255)
                                                                       gTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAg                                                                            >  2:1621042/1‑150 (MQ=255)
                                                                        tCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAgg                                                                           <  2:1382904/150‑1 (MQ=255)
                                                                             cTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGt                                                                      <  2:268187/150‑1 (MQ=255)
                                                                             cTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGt                                                                      >  1:339437/1‑150 (MQ=255)
                                                                              tGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTa                                                                     <  1:641684/150‑1 (MQ=255)
                                                                                   ttATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATc                                                                >  2:1712052/1‑150 (MQ=255)
                                                                                   ttATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATc                                                                >  2:885942/1‑150 (MQ=255)
                                                                                        caacaaAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAACGTTCGTAATATCGCCCg                                                           >  2:1309327/1‑150 (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:147330/150‑1 (MQ=255)
                                                                                                                      tACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACg                             <  2:652092/150‑1 (MQ=255)
                                                                                                                      tACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACg                             >  1:1515924/1‑150 (MQ=255)
                                                                                                                      tACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACg                             <  2:1728024/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)
                                                                                                                                                aaaTACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGac   <  1:150906/150‑1 (MQ=255)
                                                                                                                                                 aaTACCGTGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCTACACGATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTTTTCTTTATCGAATAAATTATCGaca  <  2:1270226/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GGAGCATATAGAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACA  >  NC_000913/1520798‑1521092

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: