Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 838,930 T→A 7.0% Y93F (TAT→TTT)  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,9300TA 7.0% 912.6 / 15.6 290Y93F (TAT→TTT) ybiXFe(II)‑dependent oxygenase superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (140/117);  new base A (24/6);  total (166/124)
Fisher's exact test for biased strand distribution p-value = 1.03e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 294 total aligned reads displayed.

TACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGA  >  NC_000913/838782‑839079
                                                                                                                                                    |                                                                                                                                                     
tACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATa                                                                                                                                                      >  1:1122842/1‑150 (MQ=255)
tACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATa                                                                                                                                                      <  2:872466/150‑1 (MQ=255)
tACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATa                                                                                                                                                      <  1:2119371/150‑1 (MQ=255)
tACACGATGTTGTCAGAATGTTTCATTATCGAACAGTTCGCCGGAGTGGTAGGTTTGTGGATCGCTTTTTAACAGGGTGGCAGATAGGTCAGTACGCATCCAGCCGTTTTGTGTATTGCTGAGTACTGCGCCATCCCCATGGAAACCTTa                                                                                                                                                      <  2:2159209/150‑1 (MQ=255)
   ccGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGt                                                                                                                                                   >  2:1182911/1‑150 (MQ=255)
    cGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGtt                                                                                                                                                  >  1:615442/1‑150 (MQ=255)
       tgttgtCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCa                                                                                                                                               <  1:3654731/150‑1 (MQ=255)
        gttgtCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCAt                                                                                                                                              <  1:760900/150‑1 (MQ=255)
         gtgtCCGACGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCTATGGTTTccat                                                                                                                                             <  2:535602/149‑4 (MQ=255)
            tCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGtt                                                                                                                                          <  1:1505131/150‑1 (MQ=255)
            tCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAAACATAGGTTTCATtttt                                                                                                                                          >  1:1155666/1‑150 (MQ=255)
             ccGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTc                                                                                                                                         <  2:3884059/150‑1 (MQ=255)
              cGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCt                                                                                                                                        >  1:3745435/1‑150 (MQ=255)
              cGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCt                                                                                                                                        >  2:1715419/1‑150 (MQ=255)
              cGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCt                                                                                                                                        >  1:2430798/1‑150 (MQ=255)
                      tCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGa                                                                                                                                >  2:2568110/1‑150 (MQ=255)
                      tCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGa                                                                                                                                <  1:460648/150‑1 (MQ=255)
                         ttAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTa                                                                                                                             >  2:4282073/1‑150 (MQ=255)
                          tAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAAAAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTTGATGGCTTCGTACTTAGCCATCCACAAGGAAACCATAGGTAAAAAAGTTCTGATATCAATTaa                                                                                                                            >  1:1873346/1‑150 (MQ=255)
                           aaCGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGACAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTaaa                                                                                                                           >  1:2607740/1‑150 (MQ=255)
                             cGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACa                                                                                                                         >  1:2496706/1‑150 (MQ=255)
                             cGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACa                                                                                                                         <  2:310177/150‑1 (MQ=255)
                             cGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAAACAAAGGTTTCATTGTTCTGATAGCGATTaaaaa                                                                                                                         >  1:1332091/1‑148 (MQ=255)
                              gACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCAACCACATGGAAACCATAGGTTTCATTGTTCTTATAGCGATTAAACAg                                                                                                                        >  1:254308/1‑150 (MQ=255)
                              gACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTGTTGTGGTAGGCTGCGTACTGAGCCAACCACAAGGAAACCAAAGGTTTCAATGTTCTGATAGCGATTAAACAg                                                                                                                        >  1:4154630/1‑150 (MQ=255)
                               aCCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAgc                                                                                                                       <  1:529872/150‑1 (MQ=255)
                                 cAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTTTTGATTGCTGCGTACTGCGCCATCCAAAATGAAAACAAAGGTTTAAATGTTATGATAGCGATTAAAAAgcgg                                                                                                                     >  1:679566/1‑150 (MQ=255)
                                 cAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTTTGGATGGCTGCGTACTGCGCCATACACAATGAAAACAAAGGGTTCATTGTTTTGATAGCGATTAAAAAgagg                                                                                                                     >  1:386435/1‑150 (MQ=255)
                                 cAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTTGATTGCTGCGTACTGCGCCAACCACAAGGAAACCAAAGGTTTAATTTTTATGTTATCGATTAAACAgcgg                                                                                                                     >  1:3956080/1‑150 (MQ=255)
                                 cAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAACCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACAAGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAgcgg                                                                                                                     >  1:410429/1‑150 (MQ=255)
                                  aGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAgcggc                                                                                                                    <  2:3471932/150‑1 (MQ=255)
                                  aGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAgcggc                                                                                                                    <  2:3707980/150‑1 (MQ=255)
                                      cgccgcCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTTGATGGCTGGGTACTTCGCCCAACACAAGGAAAACAAATGGTTAAATGTTATTTTATCGATTAAACAACGGCgggg                                                                                                                >  1:1935002/1‑150 (MQ=255)
                                      cgccgcCGTCGGAGCTTTGTGGATCGCTTAACAACAGCGTGGCAGATAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAACCCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTgg                                                                                                                <  2:2972771/150‑1 (MQ=255)
                                         cgcctTCTTAGCTTTGTTGATCGATTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGaaa                                                                                                             <  2:3420087/150‑1 (MQ=255)
                                         cgccgTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGCTGGCTGCGTACCGCGCCAACCCCCAGGAAACCCAAGGTTTCAATGTTCTGATAGCGAGTAAACAGCGGCGTGGaca                                                                                                             >  1:2725831/1‑148 (MQ=255)
                                          gccgTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAg                                                                                                            <  2:2866642/150‑1 (MQ=255)
                                                gtAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTAcg                                                                                                      >  2:1133837/1‑150 (MQ=255)
                                                  aGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTAcgcg                                                                                                    <  1:3314586/150‑1 (MQ=255)
                                                      ttgttgATCTCTTAAAAACAGCGTGGCAGATAGGTCAGTACGCATCCAGCCTTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCaa                                                                                                <  2:864041/150‑1 (MQ=255)
                                                      ttGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCaa                                                                                                >  2:2411955/1‑150 (MQ=255)
                                                       tgtgGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGGGGACAGGGTACGCGGCAAg                                                                                               >  2:4076146/1‑150 (MQ=255)
                                                         tgGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGc                                                                                             <  2:2469126/150‑1 (MQ=255)
                                                            aTCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCcgc                                                                                          <  1:3504638/150‑1 (MQ=255)
                                                               gCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCCTTGCcgcgta                                                                                       <  2:2230534/150‑3 (MQ=255)
                                                               gCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGc                                                                                       >  2:4241719/1‑150 (MQ=255)
                                                                 ttAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGGTGGCTGCGTGCTGCGCCAACCACAAGGGAACCCAAGGGTTGATTGTTTTTTTAGAGAGTTAAAAGCGGCGGGGGAAGGGTACGCGGCCAGGCCGCGGCaa                                                                                     >  1:4055127/1‑150 (MQ=255)
                                                                 ttAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTGTTGTTGATTGCTGCGTACTTCGCCATCCAAAAGGAAAACAAAGGTTTAATTTTTGTGATAGCGTTTTAACAGGGGCGGGGGTAGGGGTAGCGGCCAGGCCGCGGCaa                                                                                     >  1:1309085/1‑150 (MQ=255)
                                                                   aaaaaCAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTTTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAg                                                                                   >  2:809025/1‑150 (MQ=255)
                                                                       aCAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATa                                                                               <  2:72521/150‑1 (MQ=255)
                                                                         aGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAATAATaaa                                                                             <  2:60905/150‑1 (MQ=255)
                                                                          gCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAg                                                                            <  1:1741357/150‑1 (MQ=255)
                                                                            gTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCa                                                                          <  1:337093/150‑1 (MQ=255)
                                                                                 aGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGTTGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGTATAAAGCATGTTg                                                                     <  2:3263918/150‑1 (MQ=255)
                                                                                  gAAAGGTCAGTACGCATCCAGTCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTgg                                                                    <  2:3980438/150‑1 (MQ=255)
                                                                                      ggTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTaa                                                                <  1:1674977/150‑1 (MQ=255)
                                                                                       gTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAAc                                                               <  2:3745435/150‑1 (MQ=255)
                                                                                         cAGTACGCATCCAGCCATTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACcg                                                             >  1:3679639/1‑150 (MQ=255)
                                                                                            tACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGt                                                          <  1:2220575/150‑1 (MQ=255)
                                                                                             aCGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGtt                                                         <  1:2843824/150‑1 (MQ=255)
                                                                                                cATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTAAGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAg                                                      >  2:923810/1‑150 (MQ=255)
                                                                                                  tCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGACAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCa                                                    <  2:4156749/150‑1 (MQ=255)
                                                                                                     aGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCt                                                 >  2:1162095/1‑150 (MQ=255)
                                                                                                     aGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCt                                                 >  2:1743201/1‑150 (MQ=255)
                                                                                                     aGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGGTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCt                                                 >  2:1635047/1‑150 (MQ=255)
                                                                                                          ttttGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCAGAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCAttt                                            >  2:2696771/1‑150 (MQ=255)
                                                                                                              gtgGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCa                                        <  1:992980/150‑1 (MQ=255)
                                                                                                              gtgGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCa                                        <  1:3204885/150‑1 (MQ=255)
                                                                                                               tgGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAg                                       >  2:3340033/1‑150 (MQ=255)
                                                                                                               tgGATGGCTGAGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATATCGATTAAACAGCTGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAATAATAAAGCATGTTGGTTAACCGCGTTAAGCACATCATTTTGCAg                                       >  2:1902534/1‑150 (MQ=255)
                                                                                                                  aTGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGggc                                    >  2:3472119/1‑150 (MQ=255)
                                                                                                                   tGGCTGCGGACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGggcg                                   <  2:3595369/150‑1 (MQ=255)
                                                                                                                    ggCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGggcgg                                  >  1:93914/1‑150 (MQ=255)
                                                                                                                    ggCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGggcgg                                  >  2:1783819/1‑150 (MQ=255)
                                                                                                                     gCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGggcggc                                 >  1:2373701/1‑150 (MQ=255)
                                                                                                                       tGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGt                               >  1:1553175/1‑150 (MQ=255)
                                                                                                                       tGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGt                               >  1:1005462/1‑150 (MQ=255)
                                                                                                                           tACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATaa                           >  2:1098874/1‑150 (MQ=255)
                                                                                                                              tGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGt                        <  1:322028/150‑1 (MQ=255)
                                                                                                                               gcgcCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGtg                       <  2:1427910/150‑1 (MQ=255)
                                                                                                                                  ccATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGtgctg                    >  1:291678/1‑150 (MQ=255)
                                                                                                                                   cATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTACCCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGtgctgc                   >  2:2740937/1‑150 (MQ=255)
                                                                                                                                   cATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGtgctgc                   >  1:3457215/1‑150 (MQ=255)
                                                                                                                                    aTCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCg                  >  1:2531701/1‑150 (MQ=255)
                                                                                                                                    aTCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCg                  >  1:771270/1‑150 (MQ=255)
                                                                                                                                     tCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCgg                 <  2:2257735/150‑1 (MQ=255)
                                                                                                                                     tCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCgg                 <  1:1404374/150‑1 (MQ=255)
                                                                                                                                       cacaTGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGgt               >  2:208973/1‑150 (MQ=255)
                                                                                                                                        acaTGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGgtg              >  2:351598/1‑150 (MQ=255)
                                                                                                                                        acaTGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGACGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGTCGTATAACGTGCTGCGGgtg              >  2:46075/1‑150 (MQ=255)
                                                                                                                                         caTGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGgtgt             <  1:4199808/150‑1 (MQ=255)
                                                                                                                                          aTGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTc            <  2:200632/150‑1 (MQ=255)
                                                                                                                                             gAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGAc         >  1:1636260/1‑150 (MQ=255)
                                                                                                                                             gAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGAc         <  1:226461/150‑1 (MQ=255)
                                                                                                                                               aaCCATAGTTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGAGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCt       <  1:297471/150‑1 (MQ=255)
                                                                                                                                               aaCCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCt       >  2:1831492/1‑150 (MQ=255)
                                                                                                                                                 ccATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGt     >  2:379878/1‑150 (MQ=255)
                                                                                                                                                 ccATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGAACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGt     >  2:1566759/1‑150 (MQ=255)
                                                                                                                                                    tAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGa  >  2:385342/1‑150 (MQ=255)
                                                                                                                                                    tAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGa  <  2:1347037/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGA  >  NC_000913/838782‑839079

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: