Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,520,944 A→T 8.0% intergenic (+17/+19) 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,9440AT 8.0% 330.3 / 12.9 108intergenic (+17/+19)mcbR/yncDDNA‑binding transcriptional dual regulator McbR/putative TonB‑dependent outer membrane receptor
Reads supporting (aligned to +/- strand):  ref base A (56/43);  new base T (6/3);  total (62/46)
Fisher's exact test for biased strand distribution p-value = 7.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 111 total aligned reads displayed.

GGAGCATATAGAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACAC  >  NC_000913/1520798‑1521093
                                                                                                                                                  |                                                                                                                                                     
ggAGCATATAGAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATAc                                                                                                                                                    >  2:600463/1‑150 (MQ=255)
      tataGAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCa                                                                                                                                              >  2:1009330/1‑150 (MQ=255)
          gAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGTATTTGCCGGGCAttta                                                                                                                                          <  2:1969604/150‑4 (MQ=255)
          gAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCaa                                                                                                                                          >  1:1970661/1‑150 (MQ=255)
          gAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCaa                                                                                                                                          >  2:1125113/1‑150 (MQ=255)
           aaaaCTATCATCTGTTACTCGCTGCGTTAAAACCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGTATTTGCCGGGCAtttac                                                                                                                                         <  2:1471863/150‑5 (MQ=255)
               cTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTaa                                                                                                                                     >  1:1608486/1‑150 (MQ=255)
                 atcatcTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAAtt                                                                                                                                   <  1:612613/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                                                                                                                      >  2:1103571/1‑150 (MQ=255)
                              cGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCa                                                                                                                      >  1:652092/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:153369/1‑150 (MQ=255)
                                          aGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCAt                                                                                                          >  1:147330/1‑150 (MQ=255)
                                          aGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCAt                                                                                                          >  1:570049/1‑150 (MQ=255)
                                            cAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGc                                                                                                        >  2:1635826/1‑150 (MQ=255)
                                            cAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGc                                                                                                        >  2:1131165/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)
                                                                 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)
                                                                        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)
                                                                                     aTGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGc                                                               >  1:22023/1‑150 (MQ=255)
                                                                                          acaaAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGtt                                                          >  2:820111/1‑150 (MQ=255)
                                                                                          acaaAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGtt                                                          >  2:908282/1‑150 (MQ=255)
                                                                                          acaaAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGtt                                                          <  1:600463/150‑1 (MQ=255)
                                                                                           caaAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGttt                                                         >  2:1048937/1‑150 (MQ=255)
                                                                                           caaAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGttt                                                         <  2:848986/150‑1 (MQ=255)
                                                                                           caaAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGttt                                                         >  2:56902/1‑150 (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)
                                                                                                       aTTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACa                                             >  2:150906/1‑150 (MQ=255)
                                                                                                        ttttATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACaa                                            <  1:1135893/150‑1 (MQ=255)
                                                                                                         tttATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAAttt                                                                                          <  2:1987445/103‑1 (MQ=255)
                                                                                                          ttATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATg                                          <  1:1076300/150‑1 (MQ=255)
                                                                                                          ttATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATg                                          <  2:156400/150‑1 (MQ=255)
                                                                                                                 aGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAAc                                   <  1:1185702/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:1728024/150‑1 (MQ=255)
                                                                                                                      tACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACg                              <  2:652092/150‑1 (MQ=255)
                                                                                                                      tACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACg                              >  1:1515924/1‑150 (MQ=255)
                                                                                                                        cAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGAc                                                      <  1:1209456/124‑1 (MQ=255)
                                                                                                                        cAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGAc                                                      >  2:1209456/1‑124 (MQ=255)
                                                                                                                        cAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAAc                               >  1:415312/1‑147 (MQ=255)
                                                                                                                        cAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAAc                               <  2:415312/147‑1 (MQ=255)
                                                                                                                           tCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTc                         <  1:1107800/150‑1 (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)
                                                                                                                                    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)
                                                                                                                                                  aTACCGTGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGacac  <  2:1843688/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GGAGCATATAGAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACAC  >  NC_000913/1520798‑1521093

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: