Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,520,941 C→T 6.8% intergenic (+14/+22) mcbR → / ← yncD colanic acid and biofilm gene transcriptional regulator, MqsR‑controlled/putative iron outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,520,9410CT 6.8% 536.9 / 19.3 162intergenic (+14/+22)mcbR/yncDcolanic acid and biofilm gene transcriptional regulator, MqsR‑controlled/putative iron outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (82/69);  new base T (5/6);  total (87/75)
Fisher's exact test for biased strand distribution p-value = 7.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 167 total aligned reads displayed.

GGAGCATATAGAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGA  >  NC_000913/1520798‑1521090
                                                                                                                                               |                                                                                                                                                     
ggAGCATATAGAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATAc                                                                                                                                                 <  2:2001043/150‑1 (MQ=255)
 gAGCATATAGAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATAcc                                                                                                                                                <  1:911168/150‑1 (MQ=255)
  aGCATATAGAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCg                                                                                                                                               <  2:3582045/150‑1 (MQ=255)
   gCATATAGAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCgg                                                                                                                                              <  1:3235170/150‑1 (MQ=255)
   gCATATAGAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCgg                                                                                                                                              <  2:200118/150‑1 (MQ=255)
     atatAGAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGTATTTGCCGGGc                                                                                                                                            <  2:1100757/150‑1 (MQ=255)
     atatAGAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGc                                                                                                                                            >  1:222516/1‑150 (MQ=255)
      tataGAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCa                                                                                                                                           >  1:1239729/1‑150 (MQ=255)
        taGAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGTATTTGCCGGGCAtt                                                                                                                                         <  1:3214443/150‑2 (MQ=255)
        taGAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGTATTTGCCGGGCAtt                                                                                                                                         <  1:2738484/150‑2 (MQ=255)
         aGAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATAGAATCGTTAAAGTAAATGCCCGGCAATTACAGGGCATct                                                                                                                                        >  2:844725/1‑149 (MQ=255)
          gAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCaa                                                                                                                                       <  2:2061721/150‑1 (MQ=255)
           aaaaCTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAg                                                                                                                                      <  1:1092242/150‑1 (MQ=255)
               cTATCATCTGTGACTCGCTGCGTTAAAAGCAACAGACAACGAGGGATGCAGACATTGTCTTGCTGAAATTCTGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTaa                                                                                                                                  <  2:2793792/150‑1 (MQ=255)
                tATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAAt                                                                                                                                 <  1:2652837/150‑1 (MQ=255)
                    atcTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACt                                                                                                                             >  2:1409986/1‑150 (MQ=255)
                    atcTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACt                                                                                                                             >  2:367596/1‑150 (MQ=255)
                      cTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCa                                                                                                                           >  1:3575937/1‑150 (MQ=255)
                        gTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCaaa                                                                                                                         >  1:353788/1‑150 (MQ=255)
                         ttACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAAt                                                                                                                        <  1:3361528/150‑1 (MQ=255)
                           aCTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCGGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAAtct                                                                                                                      >  2:2899999/1‑150 (MQ=255)
                           aCTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAAtct                                                                                                                      <  1:2839853/150‑1 (MQ=255)
                           aCTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAAtct                                                                                                                      <  1:3381093/150‑1 (MQ=255)
                            cTCGCTGCGTTAAAAGCTAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAAtctc                                                                                                                     <  2:3566571/150‑1 (MQ=255)
                              cGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCa                                                                                                                   <  1:3695022/150‑1 (MQ=255)
                                      tAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAatat                                                                                                           <  1:1970946/150‑1 (MQ=255)
                                        aaaGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTc                                                                                                         >  2:2104231/1‑150 (MQ=255)
                                        aaaGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTc                                                                                                         <  2:197703/150‑1 (MQ=255)
                                        aaaGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTc                                                                                                         <  2:1994660/150‑1 (MQ=255)
                                            cAAAAGACACAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGc                                                                                                     <  1:1047975/150‑1 (MQ=255)
                                            cAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGc                                                                                                     <  1:2964247/150‑1 (MQ=255)
                                            cAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGc                                                                                                     >  1:552372/1‑150 (MQ=255)
                                            cAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGc                                                                                                     >  2:1622535/1‑150 (MQ=255)
                                            cAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGc                                                                                                     >  2:2922354/1‑150 (MQ=255)
                                            cAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGc                                                                                                     <  2:3033815/150‑1 (MQ=255)
                                            cAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGc                                                                                                     >  2:3611052/1‑150 (MQ=255)
                                            cAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGc                                                                                                     >  2:399451/1‑150 (MQ=255)
                                             aaaaGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGcc                                                                                                    >  2:2189884/1‑150 (MQ=255)
                                             aaaaGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGcc                                                                                                    >  1:1848834/1‑150 (MQ=255)
                                             aaaaGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGcc                                                                                                    >  1:2259841/1‑150 (MQ=255)
                                                      agagGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATa                                                                                           >  1:2655424/1‑150 (MQ=255)
                                                          ggATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAAttt                                                                                       >  2:167739/1‑150 (MQ=255)
                                                                aGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCg                                                                                 <  1:2689062/150‑1 (MQ=255)
                                                                  aCATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCggg                                                                               >  2:1418186/1‑150 (MQ=255)
                                                                     ttGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGa                                                                            >  1:1408916/1‑150 (MQ=255)
                                                                      tGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCTCGCAATATGCATGCCGACACCATAATTTCGTCGGGGCGaa                                                                           >  1:2295663/1‑150 (MQ=255)
                                                                      tGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGaa                                                                           <  1:1014336/150‑1 (MQ=255)
                                                                      tGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGaa                                                                           <  1:412756/150‑1 (MQ=255)
                                                                         cTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGt                                                                        <  2:519153/150‑1 (MQ=255)
                                                                         cTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGt                                                                        <  1:3885280/150‑1 (MQ=255)
                                                                             cttaaaTTATGCAACAAAATATTGCTATTATATACCAGCAATACAATCGTTAAACTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGt                                                                    <  2:3101615/147‑1 (MQ=255)
                                                                              tGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCCAGGTAATAACTCAAATCTCCACGCAATATTCATGCCGACACCATAATATCGTCCGGGCGAAGGTTCGTa                                                                   >  1:2237846/1‑150 (MQ=255)
                                                                              tGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTa                                                                   >  2:705457/1‑150 (MQ=255)
                                                                               gAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGACGGTTCGTaa                                                                  >  2:636896/1‑150 (MQ=255)
                                                                               gAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTaa                                                                  >  2:3769479/1‑150 (MQ=255)
                                                                                aaaTTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAat                                                                 >  1:2993905/1‑150 (MQ=255)
                                                                                aaaTTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAat                                                                 >  2:2618635/1‑150 (MQ=255)
                                                                                  aTTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAatat                                                               >  2:3084700/1‑150 (MQ=255)
                                                                                  aTTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAatat                                                               >  2:2407446/1‑150 (MQ=255)
                                                                                   ttATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATc                                                              <  2:2090643/150‑1 (MQ=255)
                                                                                     aTGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGc                                                            <  1:1293705/150‑1 (MQ=255)
                                                                                       gCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGccc                                                          >  2:3930093/1‑150 (MQ=255)
                                                                                        caacaaAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCg                                                         <  2:1445742/150‑1 (MQ=255)
                                                                                        aaacaaAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCg                                                         <  1:1694404/149‑1 (MQ=255)
                                                                                         aacaaAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGt                                                        >  1:3336833/1‑150 (MQ=255)
                                                                                         aacaaAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGt                                                        >  1:707611/1‑150 (MQ=255)
                                                                                          acaaAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGtt                                                       >  1:582972/1‑150 (MQ=255)
                                                                                          acaaAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGtt                                                       >  2:2495992/1‑150 (MQ=255)
                                                                                           caaAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGttt                                                      >  1:2441213/1‑150 (MQ=255)
                                                                                           caaAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGttt                                                      >  1:2441735/1‑150 (MQ=255)
                                                                                            aaaaTATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATGTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTg                                                     >  1:3496828/1‑150 (MQ=255)
                                                                                             aaaTATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGa                                                    <  1:2573281/150‑1 (MQ=255)
                                                                                                     cTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGa                                            >  2:2255533/1‑150 (MQ=255)
                                                                                                     cTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGa                                            <  1:2104231/150‑1 (MQ=255)
                                                                                                      tATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGAc                                           >  2:1796739/1‑150 (MQ=255)
                                                                                                       aTTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACa                                          <  1:1301048/150‑1 (MQ=255)
                                                                                                           tataCCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGa                                      <  2:3848877/150‑1 (MQ=255)
                                                                                                           tataCCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGa                                      <  2:2822411/150‑1 (MQ=255)
                                                                                                           tataCCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGa                                      <  2:129303/150‑1 (MQ=255)
                                                                                                               ccAGCATTACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGa                                  <  2:1935671/150‑1 (MQ=255)
                                                                                                                 agaaaTACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATAAACGTAATTACTCAAATCTCCACGCAATATTCATGCCGAAACCATAATAACGTCAGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCAATGACAATGACAGAAc                                >  2:225286/4‑150 (MQ=255)
                                                                                                                 aGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAAc                                <  1:1541841/150‑1 (MQ=255)
                                                                                                                 aGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAAc                                >  1:1805743/1‑150 (MQ=255)
                                                                                                                 aGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAAc                                >  1:3854473/1‑150 (MQ=255)
                                                                                                                  gCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAAcc                               >  1:1127685/1‑150 (MQ=255)
                                                                                                                            cGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCt                     <  2:339263/150‑1 (MQ=255)
                                                                                                                               tAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTa                  >  1:802655/1‑150 (MQ=255)
                                                                                                                                 aaGTAAATGCCCGGCAAATACCGTGCATCAACGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATc                <  2:1288642/150‑1 (MQ=255)
                                                                                                                                   gTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGa              >  2:2335498/1‑150 (MQ=255)
                                                                                                                                    ttgaTGCCCGGTATTTGCCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTa                  >  2:1761971/4‑145 (MQ=255)
                                                                                                                                     aaaTGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAAt            <  1:208404/150‑1 (MQ=255)
                                                                                                                                        tGCCCGGTATTTGCCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATaaa         >  2:646549/1‑150 (MQ=255)
                                                                                                                                         gCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAAt        >  1:1178624/1‑150 (MQ=255)
                                                                                                                                            cGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTAt     >  1:1671594/1‑150 (MQ=255)
                                                                                                                                             ggCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATc    <  1:2899999/150‑1 (MQ=255)
                                                                                                                                              gCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCg   >  1:3157136/1‑150 (MQ=255)
                                                                                                                                               cAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGa  >  1:382547/1‑150 (MQ=255)
                                                                                                                                               cAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGa  <  2:2867650/150‑1 (MQ=255)
                                                                                                                                               cAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGa  <  2:1894679/150‑1 (MQ=255)
                                                                                                                                               cAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGa  >  1:784361/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
GGAGCATATAGAAAACTATCATCTGTTACTCGCTGCGTTAAAAGCAAAAGACAAAGAGGGATGCAGACATTGTCTTGCTGAAATTATGCAACAAAATATTGCTATTTTATACCAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGA  >  NC_000913/1520798‑1521090

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