Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,491,705 A→C 14.2% intergenic (‑85/+20) 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,7050AC14.2% 169.4 / 15.9 80intergenic (‑85/+20)yhfG/ppiAputative antitoxin for Fic/peptidyl‑prolyl cis‑trans isomerase A (rotamase A)
Reads supporting (aligned to +/- strand):  ref base A (20/48);  new base C (5/7);  total (25/55)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATTACGATTACGCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAA  >  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)
                atgCTCCCTGATACGGGACTTTTTCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGTAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCAt                                                                                                                                  <  2:1107169/148‑1 (MQ=255)
                aaGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCAt                                                                                                                                  <  1:642954/150‑1 (MQ=255)
                 cgCTCCCAGAGACGGGACTTTAGCGTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGAGTCCATATGCCACAACTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCAtc                                                                                                                                 <  1:283419/149‑2 (MQ=255)
                 aGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCCtt                                                                                                                                 >  1:27132/1‑150 (MQ=255)
                      ccAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGc                                                                                                                            <  1:1164178/150‑1 (MQ=255)
                        cgagaCGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAg                                                                                                                          <  2:103703/149‑1 (MQ=255)
                         gagaCGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAAGCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAgg                                                                                                                         <  2:579205/150‑1 (MQ=255)
                              gggACTTTTGCTTATCGGTTAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGTGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACttt                                                                                                                    <  2:1703249/150‑1 (MQ=255)
                                gACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAAATTGCCCGCGCGAGAAATAATTAAGGGAgggcgtttg                                                                                                                  >  1:782271/1‑143 (MQ=255)
                                   atttGCTTATCGGTGAGATTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTAAAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCgg                                                                                                               <  1:340659/149‑1 (MQ=255)
                                       gCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAg                                                                                                           <  1:1704882/150‑1 (MQ=255)
                                        cTTATCCGTTAGGTGCTTCACTTGAGGCCTCCCTGAAAGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGATACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAgg                                                                                                          <  2:1407319/150‑1 (MQ=255)
                                         ttATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGa                                                                                                         <  2:1289724/150‑1 (MQ=255)
                                            tCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATaa                                                                                                      <  1:1253130/150‑1 (MQ=255)
                                             cGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCCAGCTTTCCCCCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATaaa                                                                                                     >  1:1061424/1‑149 (MQ=255)
                                                       cttcACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGCTTGCCCGCGCCACAAATGAGTACGGCAGGGCTGTATCGGAAAGGATAACTCCCCGgttt                                                                                           >  2:163662/1‑150 (MQ=255)
                                                       cttcACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGtttt                                                                                           <  2:973386/150‑1 (MQ=255)
                                                        ttcACTTGATGCCTCCCTGACCGTGTCCATTTGCCACAAGTATACGCGGCAAATCCGGGGTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTg                                                                                          <  1:1941066/150‑1 (MQ=255)
                                                        ttcACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTg                                                                                          <  2:1254158/150‑1 (MQ=255)
                                                         tcACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTTCCCGCGCGAGAAATCATTACGGCAGGACTTTGGCGGAAAGGATAACTACCGGTTTTgc                                                                                         >  2:1383334/1‑149 (MQ=255)
                                                          cACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGAc                                                                                        <  1:603406/150‑1 (MQ=255)
                                                           aCTTGATGCCTCCCTGAACGTTTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACg                                                                                       <  1:1877786/150‑1 (MQ=255)
                                                            cTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACgg                                                                                      <  2:753595/150‑1 (MQ=255)
                                                              tGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTTCGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGca                                                                                    >  1:1161498/1‑150 (MQ=255)
                                                                  gCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAATATTTGCCCCGGGTGAAATACATAAGGGCGGACGTGTACGGAAAAGAAATCTACCCGTTTTGTCGGGACAtt                                                                                >  2:1225978/1‑150 (MQ=255)
                                                                      cccTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTgg                                                                            <  1:1748574/150‑1 (MQ=255)
                                                                        cTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTa                                                                          <  1:1856110/150‑1 (MQ=255)
                                                                        cTGAACGTGTCCATTACCAACAAGTATAAGCGGCAAATCCGGGGAACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTCGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTa                                                                          <  2:1438650/150‑1 (MQ=255)
                                                                         tGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTAc                                                                         <  1:1582328/150‑1 (MQ=255)
                                                                            ccGTGGCCATTTGCCACAAGAATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCACGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGa                                                                      <  2:1532258/149‑1 (MQ=255)
                                                                             cGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGAc                                                                     <  2:1065117/150‑1 (MQ=255)
                                                                              gtgtCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGAcc                                                                    <  2:458565/150‑1 (MQ=255)
                                                                                 tCCATTTGCCACAAGTATAAGCTGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATTACTACCGGTTTTGACGGCACATTCTGGTACGGACCaag                                                                 <  1:363504/150‑2 (MQ=255)
                                                                                  ccATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAACATTTCCCCCGCCGGAAAAAATTAAGCCACGCCATTTTCGGAAAGGTAAACGAACGGTTTTTACAGCACATTATGGGACGGCCCacca                                                                >  2:1233568/1‑147 (MQ=255)
                                                                                    aTTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTc                                                              <  2:1060194/150‑1 (MQ=255)
                                                                                     tttGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCa                                                             <  2:571446/150‑1 (MQ=255)
                                                                                       ttccaacaGTATAAGCGTCAAATCCGGGTTACGCCGGGAGAGCATAAGCGCGGGCATGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTTCCGGTTTTGACGGGACATTCTGGTACGGACCAACGTCATg                                                           <  1:1783855/143‑1 (MQ=255)
                                                                                       tGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATg                                                           <  1:1959108/150‑1 (MQ=255)
                                                                                        gCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAATGAAAACTACCGGTTTTGACGGCACATTCTTGTACGGACCACCGTCATGa                                                          >  2:1185441/1‑150 (MQ=255)
                                                                                          cacaAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTAAGGCAGGACTTTAGCGGAAAGGATAACTAACGGGTTTTACGGCCAATTCTGGTAAGGACCAACGTCATGAGt                                                        >  2:1138864/1‑150 (MQ=255)
                                                                                               gTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCa                                                   <  2:1896010/150‑1 (MQ=255)
                                                                                                 ataAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCAcc                                                 <  1:1810787/150‑1 (MQ=255)
                                                                                                      cGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGGGGGCAAGATTGCCCGCGCGAGAAATCATTCCGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGa                                            <  1:77179/150‑1 (MQ=255)
                                                                                                        gCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGaaa                                          <  2:1673701/150‑1 (MQ=255)
                                                                                                           aaTCCGGGTTACGCCGGGAGAGCAGATGCGCGGGCAAGATGGCCCGCGGGAGAAATCATTACGGCAGGACTTTAGCGGATAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCt                                       <  2:310287/150‑1 (MQ=255)
                                                                                                           aaTCCGGGGTCCGCCGGGAGAGCAGACGCGCGGGCTAGCTTGCCCGGGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATATCTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCt                                       <  2:272675/150‑1 (MQ=255)
                                                                                                            atcagttttACGCCGGGAGAGCAGAAGCGCTGTCAAGATTGCCCGCGCGAGAAATCATTACGGCATGTCTTTAGCGTAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCTACGTCATGAGTCGGCACCTGGGAATTCtt                                      <  2:810311/143‑1 (MQ=255)
                                                                                                            aTCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACCGCAGGACATTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGACGGCACCTGGGAAATCtt                                      >  1:277914/1‑150 (MQ=255)
                                                                                                             tCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAAATACCGGTTTTGACGGCACATTATGGTACGGACCAAAGTCATGAGTCGGCACCTGGGAAATCTTa                                     >  1:6382/1‑150 (MQ=255)
                                                                                                                 ggTTCCGCCGGGAGAGCAGCAGCGCGGGCAAGATGGCCCGCGCGAGAAATCATTCCCGCAGGACTCTAGCGGAAAGGATAACTACCCGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCgg                                 <  2:185306/150‑1 (MQ=255)
                                                                                                                   ttACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAACTCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCa                               <  2:23418/150‑1 (MQ=255)
                                                                                                                   ttACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCa                               <  2:1858814/150‑1 (MQ=255)
                                                                                                                      cGCCGGGAGAGCAGAAGCGCGGGCAAGAATGCCCGCGCGAGAAATCATTACGGCAGGACTTTTACGGAAAGGATAAATACCCGTTTTGACGGCACATTCTGGTACGGCCCAACGTGATGAGTCGGCACCGGGGAAAGCATATCGGGAAcc                            >  2:1169392/1‑149 (MQ=255)
                                                                                                                        ccGGGAGAGCAGAAGCGCGGGGAAGATTTCCCGCGCGAGCAATAAATAAGGCAGGCCTGTAGCGGAAAGGATAAAGACCGGGTTTGACGGCCCAGTCTGGTACGGCCCCACCGCATTAGTGGGCCCCAGGGGTATCTTTTAGGCCACCTc                          >  2:572949/1‑150 (MQ=255)
                                                                                                                             agagCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTTGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGc                     <  1:1429191/150‑1 (MQ=255)
                                                                                                                               atttctcgcgcGGGCAATCTTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTAGCGGTTTTTACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCCCCTGGGAAATCTTATCGGCAACGTCCATGCCt                   >  2:413653/8‑150 (MQ=255)
                                                                                                                                  aGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGTCGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTc                <  2:1476620/150‑1 (MQ=255)
                                                                                                                                  aGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTc                <  2:782271/150‑1 (MQ=255)
                                                                                                                                     cgcgcGGGCAATCTTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTcacc             >  1:1090753/2‑150 (MQ=255)
                                                                                                                                     aGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTcacc             >  1:1471491/1‑150 (MQ=255)
                                                                                                                                      gcgcGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCGTGcaccc            >  1:1745980/1‑149 (MQ=255)
                                                                                                                                         cgGGCAAGATTGCCCGCGCGAGAACTCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGCCCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCGTGCCTTTcaccacgt         >  2:291620/1‑148 (MQ=255)
                                                                                                                                            gCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTAc      >  1:385585/1‑150 (MQ=255)
                                                                                                                                            gCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTAc      <  1:517012/150‑1 (MQ=255)
                                                                                                                                            gCAAGATGGCCCGCGCGAGAAATCATTACGGCATGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACTGACCAACTTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTAc      <  2:275003/150‑1 (MQ=255)
                                                                                                                                                gATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCaaa  >  2:1380962/1‑150 (MQ=255)
                                                                                                                                                gATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCaaa  >  2:804394/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
AATTACGATTACGCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAA  >  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: