Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 838,921 T→A 5.9% H96L (CAT→CTT)  ybiX ← Fe(II)‑dependent oxygenase superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913838,9210TA 5.9% 2241.3 / 34.9 661H96L (CAT→CTT) ybiXFe(II)‑dependent oxygenase superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (268/338);  new base A (53/2);  total (321/340)
Fisher's exact test for biased strand distribution p-value = 1.84e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 670 total aligned reads displayed.

CCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCG  >  NC_000913/838772‑839070
                                                                                                                                                     |                                                                                                                                                     
ccGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGTAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGGGTAATGCGCCATCCACAt                                                                                                                                                       >  1:5294368/1‑150 (MQ=255)
   ggAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGGTTTGTGGATGGCTGCGTACGGCGCCATCCACATGGa                                                                                                                                                    >  1:10840849/1‑150 (MQ=255)
     aGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGaaa                                                                                                                                                  <  1:10370568/150‑1 (MQ=255)
       ttttACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAAcc                                                                                                                                                >  1:1416258/1‑150 (MQ=255)
        tttACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCa                                                                                                                                               >  1:6334881/1‑150 (MQ=255)
         ttACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACAGCGCCATCCACAAGGAAACCAt                                                                                                                                              >  1:9438763/1‑150 (MQ=255)
          tACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATa                                                                                                                                             <  1:8895940/150‑1 (MQ=255)
           aCCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAg                                                                                                                                            <  2:133364/150‑1 (MQ=255)
           aCCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAg                                                                                                                                            <  1:4449161/150‑1 (MQ=255)
            cccGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAgg                                                                                                                                           >  2:2072020/1‑150 (MQ=255)
            cccGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAgg                                                                                                                                           >  1:2558247/1‑150 (MQ=255)
              cGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGtt                                                                                                                                         >  2:2138549/1‑150 (MQ=255)
               gATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGttt                                                                                                                                        >  2:9840654/1‑150 (MQ=255)
               gATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGttt                                                                                                                                        >  2:8899978/1‑150 (MQ=255)
               gATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCAAAGGttt                                                                                                                                        >  1:7539891/1‑150 (MQ=255)
               gATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAACCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGttt                                                                                                                                        <  2:6986220/150‑1 (MQ=255)
                 tgttgtCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCTATGGTTTcc                                                                                                                                      <  2:2481555/150‑2 (MQ=255)
                  gttgtCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCAt                                                                                                                                     <  1:7076329/150‑1 (MQ=255)
                      tCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGtt                                                                                                                                 >  2:7363334/1‑150 (MQ=255)
                      tCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCAACCACAAGGAAACCATAGGTTTCATTGtt                                                                                                                                 >  1:8653514/1‑150 (MQ=255)
                       ccGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTc                                                                                                                                <  2:10498141/150‑1 (MQ=255)
                       ccGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTc                                                                                                                                >  2:460806/1‑150 (MQ=255)
                          aaGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGa                                                                                                                             <  1:11021566/150‑1 (MQ=255)
                           aGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGAt                                                                                                                            <  1:5079532/150‑1 (MQ=255)
                              tgtCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGc                                                                                                                         <  2:1466006/150‑1 (MQ=255)
                              tgtCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGc                                                                                                                         <  2:3864398/150‑1 (MQ=255)
                                tCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGa                                                                                                                       <  2:10178360/150‑1 (MQ=255)
                                 cATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGAt                                                                                                                      >  2:3721570/1‑150 (MQ=255)
                                   ttAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTa                                                                                                                    <  2:6686802/150‑1 (MQ=255)
                                   ttAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCAACCACAAGGAAACCATAGGTTTCATTGTTATGATAGCGATTa                                                                                                                    >  1:10572529/1‑150 (MQ=255)
                                      aCGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAAc                                                                                                                 >  2:3206984/1‑150 (MQ=255)
                                      aCGACCAGATCGCCGCCGTCGTAGCTTTGTGGATCGCTAAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACAGCGCCATCCACCTGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAAc                                                                                                                 >  2:10268169/1‑150 (MQ=255)
                                          ccAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAgcg                                                                                                             <  1:401942/150‑1 (MQ=255)
                                            aGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAgcggc                                                                                                           <  1:1859058/150‑1 (MQ=255)
                                            aGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACAAGGAAACCATAGGTTTCATTGTTATGATAGCGATTAAACAgaggc                                                                                                           >  1:6638383/1‑150 (MQ=255)
                                               tCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAAAAGCGTGGCAGAATGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCTTCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTg                                                                                                        <  1:475661/150‑1 (MQ=255)
                                                   cgccgTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGaaa                                                                                                    >  2:9159165/1‑150 (MQ=255)
                                                   cgccgTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGaaa                                                                                                    <  1:2446729/150‑1 (MQ=255)
                                                    gccgTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAg                                                                                                   <  1:1636815/150‑1 (MQ=255)
                                                    gccgTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAg                                                                                                   <  2:1081407/150‑1 (MQ=255)
                                                     ccgTCGTAGCTTTGTGGATAGCTTAAAAACAGCGTTGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGTAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAgg                                                                                                  <  1:6857335/150‑1 (MQ=255)
                                                        tcgtAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTTGATTGCTTCGTAAAGCCCCAACCACAATGAAACCAAAGGTTTAATTGTTATGATATCGATTAAAAAGCGGCGGGGAAAGGGTa                                                                                               >  1:5478779/1‑150 (MQ=255)
                                                        tcgtAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTTGATGGCTGCGTACGGCGCCAACCACAATGAAACCAAAGGTTTCATTGTTATGATATCGATTAAAAAGCGGCGTGGGAAGGGTa                                                                                               >  1:2785089/1‑150 (MQ=255)
                                                            aGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTAcgcg                                                                                           <  1:6960104/150‑1 (MQ=255)
                                                               tttGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGGTGGCTGCGTACTGCGCCCACCACAAGGAAACCAAAGGTTTAATTTTTATTATATCGATTTAAAAGCGGCGGGGGAAGGGTACGCGGcc                                                                                        >  1:4361142/1‑149 (MQ=255)
                                                               tttGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCa                                                                                        >  2:5180613/1‑150 (MQ=255)
                                                               tttGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCa                                                                                        >  2:9650420/1‑150 (MQ=255)
                                                                ttGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCaa                                                                                       <  1:7242184/150‑1 (MQ=255)
                                                                 tgtgGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAg                                                                                      <  1:4855084/150‑1 (MQ=255)
                                                                    ggATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGcc                                                                                   >  2:6170812/1‑150 (MQ=255)
                                                                       tCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCcgcg                                                                                <  2:6852867/150‑1 (MQ=255)
                                                                           ttAACAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCaa                                                                            <  2:8123318/150‑1 (MQ=255)
                                                                           ttAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCaa                                                                            >  2:378403/1‑150 (MQ=255)
                                                                           ttAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCaa                                                                            >  2:5664409/1‑150 (MQ=255)
                                                                               aaaCAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTTTGGATGGCTGCGTACTGCGCCAACCACATGGAAACCATAGGTTTAATTGTTCATATAGCGATTAAACAGCGGGGTCGAAAGGGTACGCCGCAAGGCCGCGGCAAAGaa                                                                        >  2:8401939/1‑150 (MQ=255)
                                                                                  cAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATaa                                                                     <  2:3646455/150‑1 (MQ=255)
                                                                                   aGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATaaa                                                                    <  1:3098008/150‑1 (MQ=255)
                                                                                   aGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATaaa                                                                    <  1:3687677/150‑1 (MQ=255)
                                                                                    gCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAg                                                                   <  2:7463404/150‑1 (MQ=255)
                                                                                    gCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAg                                                                   <  2:8697485/150‑1 (MQ=255)
                                                                                    gCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAg                                                                   <  2:4751352/150‑1 (MQ=255)
                                                                                     cGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGc                                                                  <  2:7129127/150‑1 (MQ=255)
                                                                                       tGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCAt                                                                <  2:8204675/150‑1 (MQ=255)
                                                                                             aaaGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGt                                                          >  2:3397478/1‑150 (MQ=255)
                                                                                             aaaGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGt                                                          >  2:11014135/1‑150 (MQ=255)
                                                                                             aaaGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACAGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGt                                                          >  2:5428602/1‑150 (MQ=255)
                                                                                              aaGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGtt                                                         >  2:8580528/1‑150 (MQ=255)
                                                                                                   cAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCCAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACcg                                                    >  2:4133183/1‑150 (MQ=255)
                                                                                                    aGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACcgc                                                   <  1:3237192/150‑1 (MQ=255)
                                                                                                     gTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACcgcg                                                  <  2:2751759/150‑1 (MQ=255)
                                                                                                      tACGCATCCAGCCGTTTTGTGGATTGCTGCGTACTTCGCCATCCACATGGAAACGATAGGTTTGATTGTTGTGTTTGGTATTAAACAGCGGCGTGGAAAGGGTTAGCGGCAAGGCCGCGGCAAAGAATAAATCATGTTGGTTATCCGCGt                                                 >  2:1889514/1‑150 (MQ=255)
                                                                                                      tACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGt                                                 <  2:7627279/150‑1 (MQ=255)
                                                                                                         gCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCa                                              <  2:6381646/150‑1 (MQ=255)
                                                                                                          cATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAg                                             <  2:4365866/150‑1 (MQ=255)
                                                                                                          cATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAg                                             <  1:7380256/150‑1 (MQ=255)
                                                                                                          cATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAg                                             <  1:5778718/150‑1 (MQ=255)
                                                                                                           aTCCAGCCGTTATGTGCATGGCTGCGTACTGCGCCAACCACAAGGACACCCAAGGGTTCATTGTTCTGATAGCGATAAAACACCGGCGTCGACAGGGAACGCGGCAACGCCGCGGCCAAGAAAAAAGCAGGGTGGATAACCCCGGTCAGc                                            >  2:3069759/1‑150 (MQ=255)
                                                                                                            tCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCa                                           >  2:1765082/1‑150 (MQ=255)
                                                                                                               aGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCt                                        >  2:6567202/1‑150 (MQ=255)
                                                                                                               aGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCt                                        >  2:5122717/1‑150 (MQ=255)
                                                                                                                      ttGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTg                                 <  2:2225642/150‑1 (MQ=255)
                                                                                                                           gATGGCTGCGTACTGCGCCATCCACATGGATACCATAGGTTTCATTGTTCTGATAGCGATTTAACAGCGGCGTGGAAAGGGTACGCGGCTAGGCCGCGGCAAAGAATACAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAggg                            <  2:5988476/150‑1 (MQ=255)
                                                                                                                             tGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGGGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCCCCTCATTTTGCAGggcg                          >  1:7730749/1‑150 (MQ=255)
                                                                                                                                cTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGgcggcg                       <  2:1181835/150‑1 (MQ=255)
                                                                                                                                 tGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCTGGGCGGCGt                      >  2:4934227/1‑150 (MQ=255)
                                                                                                                                  gCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGta                     >  1:2265686/1‑150 (MQ=255)
                                                                                                                                        tGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGt               >  1:5634089/1‑150 (MQ=255)
                                                                                                                                          cgcCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGtgc             >  1:3412426/1‑150 (MQ=255)
                                                                                                                                          cgcCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGtgc             <  1:9647192/150‑1 (MQ=255)
                                                                                                                                              aTCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCg         >  1:1317838/1‑150 (MQ=255)
                                                                                                                                              aTCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCg         >  1:9214526/1‑150 (MQ=255)
                                                                                                                                               tCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCgg        <  1:10017859/150‑1 (MQ=255)
                                                                                                                                                ccACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCggg       >  1:8275390/1‑150 (MQ=255)
                                                                                                                                                  acaTGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGgtg     >  1:9817754/1‑150 (MQ=255)
                                                                                                                                                  acaTGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGgtg     >  1:1042662/1‑150 (MQ=255)
                                                                                                                                                  acaTGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGgtg     >  2:5611300/1‑150 (MQ=255)
                                                                                                                                                  acaTGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGgtg     >  2:9407282/1‑150 (MQ=255)
                                                                                                                                                   caTGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGgtgt    <  1:2017613/150‑1 (MQ=255)
                                                                                                                                                    aTGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTc   >  1:8062880/1‑150 (MQ=255)
                                                                                                                                                     tGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCg  >  2:7809067/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCG  >  NC_000913/838772‑839070

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: