Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,491,704 G→T 12.3% intergenic (‑84/+21) yhfG ← / ← ppiA putative antitoxin for Fic/peptidyl‑prolyl cis‑trans isomerase A (rotamase A)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,491,7040GT12.3% 162.9 / 10.5 78intergenic (‑84/+21)yhfG/ppiAputative antitoxin for Fic/peptidyl‑prolyl cis‑trans isomerase A (rotamase A)
Reads supporting (aligned to +/- strand):  ref base G (16/49);  new base T (3/6);  total (23/55)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.22e-01

AATTACGATTACGCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAA  >  NC_000913/3491560‑3491853
                                                                                                                                                |                                                                                                                                                      
aaTTACGATTACGCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGc                                                                                                                                                   <  1:854389/150‑1 (MQ=255)
 attacgattacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTGCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGcc                                                                                                                                                  <  2:839071/150‑1 (MQ=255)
     cgattacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTTAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAATCTTGCCcgcg                                                                                                                                              <  2:109985/150‑1 (MQ=255)
     cgattacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGGGGGCGAAAtttgacccc                                                                                                                                              >  2:1810787/1‑143 (MQ=255)
     cgattacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAATCTTGCCcgcg                                                                                                                                              <  1:455441/150‑1 (MQ=255)
     cgattacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAATCTTGCCcgcg                                                                                                                                              <  2:1520579/150‑1 (MQ=255)
      gattacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAATCTTGCCcgcgc                                                                                                                                             <  1:182267/150‑1 (MQ=255)
       attacgCTGAAGCTCCCAGAGGCGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGCACGTGTCCATTAGCCACAAGTCTAAGCGGCAAATCCGGGGTACTCCGGGAGAGCAGAAGCGCGGGCAATCTTGCCcgcgct                                                                                                                                            <  1:893767/150‑2 (MQ=255)
       attacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATATGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGCGCAACAATGCCcccggg                                                                                                                                            >  1:810311/1‑148 (MQ=255)
          acgCTGAAGCGCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCgaga                                                                                                                                         <  2:289291/150‑1 (MQ=255)
           cgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAATCTTGCCcgcgcttctg                                                                                                                                        <  2:45564/150‑6 (MQ=255)
           cgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGCTTGCCCGCGCGAGaa                                                                                                                                        >  2:1055294/1‑150 (MQ=255)
           cgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGaa                                                                                                                                        <  1:1571559/150‑1 (MQ=255)
                atgCTCCCTGATACGGGACTTTTTCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGTAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CAt                                                                                                                                  <  2:1107169/148‑1 (MQ=255)
                aaGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CAt                                                                                                                                  <  1:642954/150‑1 (MQ=255)
                 cgCTCCCAGAGACGGGACTTTAGCGTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGAGTCCATATGCCACAACTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CAtc                                                                                                                                 <  1:283419/149‑2 (MQ=255)
                 aGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CCtt                                                                                                                                 >  1:27132/1‑150 (MQ=255)
                      ccAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGc                                                                                                                            <  1:1164178/150‑1 (MQ=255)
                        cgagaCGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAg                                                                                                                          <  2:103703/149‑1 (MQ=255)
                         gagaCGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAAGCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAgg                                                                                                                         <  2:579205/150‑1 (MQ=255)
                              gggACTTTTGCTTATCGGTTAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGTGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACttt                                                                                                                    <  2:1703249/150‑1 (MQ=255)
                                gACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAAATTGCCCGCGCGAGAAAT‑AATTAAGGGAgggcgtttg                                                                                                                  >  1:782271/1‑143 (MQ=255)
                                   atttGCTTATCGGTGAGATTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTAAAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCgg                                                                                                               <  1:340659/149‑1 (MQ=255)
                                       gCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAg                                                                                                           <  1:1704882/150‑1 (MQ=255)
                                        cTTATCCGTTAGGTGCTTCACTTGAGGCCTCCCTGAAAGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGATACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAgg                                                                                                          <  2:1407319/150‑1 (MQ=255)
                                         ttATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGa                                                                                                         <  2:1289724/150‑1 (MQ=255)
                                            tCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATaa                                                                                                      <  1:1253130/150‑1 (MQ=255)
                                             cGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCCAGCTTTCCCCCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATaaa                                                                                                     >  1:1061424/1‑149 (MQ=255)
                                                       cttcACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGCTTGCCCGCGCCACAAAT‑GAGTACGGCAGGGCTGTATCGGAAAGGATAACTCCCCGgttt                                                                                           >  2:163662/1‑150 (MQ=255)
                                                       cttcACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGtttt                                                                                           <  2:973386/150‑1 (MQ=255)
                                                        ttcACTTGATGCCTCCCTGACCGTGTCCATTTGCCACAAGTATACGCGGCAAATCCGGGGTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTg                                                                                          <  1:1941066/150‑1 (MQ=255)
                                                        ttcACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTg                                                                                          <  2:1254158/150‑1 (MQ=255)
                                                         tcACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTTCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTGGCGGAAAGGATAACTACCGGTTTTgc                                                                                         >  2:1383334/1‑149 (MQ=255)
                                                          cACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGAc                                                                                        <  1:603406/150‑1 (MQ=255)
                                                           aCTTGATGCCTCCCTGAACGTTTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACg                                                                                       <  1:1877786/150‑1 (MQ=255)
                                                            cTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACgg                                                                                      <  2:753595/150‑1 (MQ=255)
                                                              tGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTTCGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGca                                                                                    >  1:1161498/1‑150 (MQ=255)
                                                                  gCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAATATTTGCCCCGGGTGAAATACATAAGGGC‑GGACGTGTACGGAAAAGAAATCTACCCGTTTTGTCGGGACAtt                                                                                >  2:1225978/1‑150 (MQ=255)
                                                                      cccTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTgg                                                                            <  1:1748574/150‑1 (MQ=255)
                                                                        cTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTa                                                                          <  1:1856110/150‑1 (MQ=255)
                                                                        cTGAACGTGTCCATTACCAACAAGTATAAGCGGCAAATCCGGGGAACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTCGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTa                                                                          <  2:1438650/150‑1 (MQ=255)
                                                                         tGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTAc                                                                         <  1:1582328/150‑1 (MQ=255)
                                                                            ccGTGGCCATTTGCCACAAGAATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCACGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGa                                                                      <  2:1532258/149‑1 (MQ=255)
                                                                             cGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGAc                                                                     <  2:1065117/150‑1 (MQ=255)
                                                                              gtgtCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGAcc                                                                    <  2:458565/150‑1 (MQ=255)
                                                                                 tCCATTTGCCACAAGTATAAGCTGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATTACTACCGGTTTTGACGGCACATTCTGGTACGGACCaag                                                                 <  1:363504/150‑2 (MQ=255)
                                                                                  ccATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAACATTTCCCCCGCCGGAAAA‑AATTAAGCCACGCCATTTTCGGAAAGGTAAACGAACGGTTTTTACAGCACATTATGGGACGGCCCacca                                                                >  2:1233568/1‑147 (MQ=255)
                                                                                    aTTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTc                                                              <  2:1060194/150‑1 (MQ=255)
                                                                                     tttGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCa                                                             <  2:571446/150‑1 (MQ=255)
                                                                                       ttccaacaGTATAAGCGTCAAATCCGGGTTACGCCGGGAGAGCATAAGCGCGGGCATGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATAACTTCCGGTTTTGACGGGACATTCTGGTACGGACCAACGTCATg                                                           <  1:1783855/143‑1 (MQ=255)
                                                                                       tGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATg                                                           <  1:1959108/150‑1 (MQ=255)
                                                                                        gCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAATGAAAACTACCGGTTTTGACGGCACATTCTTGTACGGACCACCGTCATGa                                                          >  2:1185441/1‑150 (MQ=255)
                                                                                          cacaAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTAAGGCAGGACTTTAGCGGAAAGGATAACTAACGGGTTTTACGGCCAATTCTGGTAAGGACCAACGTCATGAGt                                                        >  2:1138864/1‑150 (MQ=255)
                                                                                               gTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCa                                                   <  2:1896010/150‑1 (MQ=255)
                                                                                                 ataAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCAcc                                                 <  1:1810787/150‑1 (MQ=255)
                                                                                                      cGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGGGGGCAAGATTGCCCGCGCGAGAAAT‑CATTCCGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGa                                            <  1:77179/150‑1 (MQ=255)
                                                                                                        gCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGaaa                                          <  2:1673701/150‑1 (MQ=255)
                                                                                                           aaTCCGGGTTACGCCGGGAGAGCAGATGCGCGGGCAAGATGGCCCGCGGGAGAAAT‑CATTACGGCAGGACTTTAGCGGATAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCt                                       <  2:310287/150‑1 (MQ=255)
                                                                                                           aaTCCGGGGTCCGCCGGGAGAGCAGACGCGCGGGCTAGCTTGCCCGGGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATATCTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCt                                       <  2:272675/150‑1 (MQ=255)
                                                                                                            atcagttttACGCCGGGAGAGCAGAAGCGCTGTCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCATGTCTTTAGCGTAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCTACGTCATGAGTCGGCACCTGGGAATTCtt                                      <  2:810311/143‑1 (MQ=255)
                                                                                                            aTCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACCGCAGGACATTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGACGGCACCTGGGAAATCtt                                      >  1:277914/1‑150 (MQ=255)
                                                                                                             tCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATAAATACCGGTTTTGACGGCACATTATGGTACGGACCAAAGTCATGAGTCGGCACCTGGGAAATCTTa                                     >  1:6382/1‑150 (MQ=255)
                                                                                                                 ggTTCCGCCGGGAGAGCAGCAGCGCGGGCAAGATGGCCCGCGCGAGAAAT‑CATTCCCGCAGGACTCTAGCGGAAAGGATAACTACCCGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCgg                                 <  2:185306/150‑1 (MQ=255)
                                                                                                                   ttACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAACT‑CATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCa                               <  2:23418/150‑1 (MQ=255)
                                                                                                                   ttACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCa                               <  2:1858814/150‑1 (MQ=255)
                                                                                                                      cGCCGGGAGAGCAGAAGCGCGGGCAAGAATGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTTACGGAAAGGATAAATACCCGTTTTGACGGCACATTCTGGTACGGCCCAACGTGATGAGTCGGCACCGGGGAAAGCATATCGGGAAcc                            >  2:1169392/1‑149 (MQ=255)
                                                                                                                        ccGGGAGAGCAGAAGCGCGGGGAAGATTTCCCGCGCGAGCAAT‑AAATAAGGCAGGCCTGTAGCGGAAAGGATAAAGACCGGGTTTGACGGCCCAGTCTGGTACGGCCCCACCGCATTAGTGGGCCCCAGGGGTATCTTTTAGGCCACCTc                          >  2:572949/1‑150 (MQ=255)
                                                                                                                             agagCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTTGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGc                     <  1:1429191/150‑1 (MQ=255)
                                                                                                                               atttctcgcgcGGGCAATCTTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATAACTAGCGGTTTTTACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCCCCTGGGAAATCTTATCGGCAACGTCCATGCCt                   >  2:413653/8‑150 (MQ=255)
                                                                                                                                  aGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGTCGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTc                <  2:1476620/150‑1 (MQ=255)
                                                                                                                                  aGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTc                <  2:782271/150‑1 (MQ=255)
                                                                                                                                     cgcgcGGGCAATCTTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTcacc             >  1:1090753/2‑150 (MQ=255)
                                                                                                                                     aGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTcacc             >  1:1471491/1‑150 (MQ=255)
                                                                                                                                      gcgcGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCGTGcaccc            >  1:1745980/1‑149 (MQ=255)
                                                                                                                                         cgGGCAAGATTGCCCGCGCGAGAACT‑CATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGCCCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCGTGCCTTTcaccacgt         >  2:291620/1‑148 (MQ=255)
                                                                                                                                            gCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTAc      >  1:385585/1‑150 (MQ=255)
                                                                                                                                            gCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTAc      <  1:517012/150‑1 (MQ=255)
                                                                                                                                            gCAAGATGGCCCGCGCGAGAAAT‑CATTACGGCATGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACTGACCAACTTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTAc      <  2:275003/150‑1 (MQ=255)
                                                                                                                                                gATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCaaa  >  2:1380962/1‑150 (MQ=255)
                                                                                                                                                gATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCaaa  >  2:804394/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                      
AATTACGATTACGCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAAT‑CATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAA  >  NC_000913/3491560‑3491853

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑