Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,987,174 C→G 5.4% D237H (GAT→CAT)  hemY ← protein HemY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,987,1740CG 5.4% 505.8 / 14.0 148D237H (GAT→CAT) hemYprotein HemY
Reads supporting (aligned to +/- strand):  ref base C (66/73);  new base G (4/4);  total (71/77)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

Only 100 of 154 total aligned reads displayed.

GATAATTTGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTAATACTTCCG  >  NC_000913/3987025‑3987323
                                                                                                                                                     |                                                                                                                                                     
gatAATTTGCCGGGCAGTATCATGATCGTCACATTCAATAAGATGGTCCGCCATTGCCACCTGCAACGCTACCTGATCACGCGGTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATc                                                                                                                                                       <  2:725695/150‑1 (MQ=255)
 atAATTTGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATcc                                                                                                                                                      <  2:1189207/150‑1 (MQ=255)
 atAATTTGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACATGCAACGCTACCTGATGACGCGTTTTCCTGGTTTGGTTTTTCGACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGGGCCTGATcc                                                                                                                                                      <  2:1106914/150‑1 (MQ=255)
  tAATTTGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCa                                                                                                                                                     >  1:954615/1‑150 (MQ=255)
    aTTTGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATc                                                                                                                                                   >  1:154070/1‑150 (MQ=255)
      ttGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAg                                                                                                                                                 <  2:843368/150‑1 (MQ=255)
      ttGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAg                                                                                                                                                 <  1:839974/150‑1 (MQ=255)
        gCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGc                                                                                                                                               <  1:433118/150‑1 (MQ=255)
        gCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGc                                                                                                                                               >  2:2332607/1‑150 (MQ=255)
         cTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGcc                                                                                                                                              <  1:747903/150‑1 (MQ=255)
         cTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGcc                                                                                                                                              >  2:1825536/1‑150 (MQ=255)
         cTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATACATCAGGcc                                                                                                                                              >  2:1502924/1‑150 (MQ=255)
          tGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCa                                                                                                                                             >  2:1801642/1‑150 (MQ=255)
           gggCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATGGATCAGGcgcg                                                                                                                                            <  2:65160/150‑4 (MQ=255)
            ggCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAAt                                                                                                                                           >  1:832479/1‑150 (MQ=255)
            ggCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAAt                                                                                                                                           <  2:448700/150‑1 (MQ=255)
             gCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTTGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATc                                                                                                                                          <  1:2118375/150‑1 (MQ=255)
              cAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCCCCTGCAACGCTACCTGATGACGCGTTTTCCGGCTGTGCGTTGTCCCCCGGATACGCAAACCTTCGCTACCGTGATCGGCACGCGCCAGATCCATCAGGCCAATcc                                                                                                                                         >  2:1363759/1‑150 (MQ=255)
              cAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATGGATCAGGcgcgtgc                                                                                                                                         <  1:558201/150‑7 (MQ=255)
                   tCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGcc                                                                                                                                    <  1:1933352/150‑1 (MQ=255)
                   tCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGcc                                                                                                                                    <  2:1457544/150‑1 (MQ=255)
                   tCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGcc                                                                                                                                    >  1:531390/1‑150 (MQ=255)
                        aTCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCtgttg                                                                                                                               <  2:1695743/150‑1 (MQ=255)
                         tCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCtgttgt                                                                                                                              >  1:676476/1‑150 (MQ=255)
                            tCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTcc                                                                                                                           <  1:674823/150‑1 (MQ=255)
                              acaTTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAg                                                                                                                         <  2:663809/150‑1 (MQ=255)
                              acaTTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAg                                                                                                                         <  2:1089651/150‑1 (MQ=255)
                               caTTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGc                                                                                                                        >  1:253539/1‑150 (MQ=255)
                                  tCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCAtt                                                                                                                     >  1:1821866/1‑150 (MQ=255)
                                  tCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCAtt                                                                                                                     <  1:1314317/150‑1 (MQ=255)
                                   cAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTg                                                                                                                    >  1:1471178/1‑150 (MQ=255)
                                     aTAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCa                                                                                                                  >  1:984728/1‑150 (MQ=255)
                                           tGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGt                                                                                                            <  2:615002/150‑1 (MQ=255)
                                            gTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGtt                                                                                                           >  1:1822230/1‑150 (MQ=255)
                                             ttCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGttc                                                                                                          <  1:2477423/150‑1 (MQ=255)
                                              tCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGttct                                                                                                         <  1:1210763/150‑1 (MQ=255)
                                               ccGCCATCGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGttctt                                                                                                        >  2:406961/1‑150 (MQ=255)
                                                 gCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTtca                                                                                                      <  2:1422824/150‑1 (MQ=255)
                                                   cATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTtcatc                                                                                                    <  1:934803/150‑1 (MQ=255)
                                                       gCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCa                                                                                                >  1:1660638/1‑150 (MQ=255)
                                                          aCCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACa                                                                                             >  1:1842184/1‑150 (MQ=255)
                                                             tGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATgg                                                                                          >  1:2236638/1‑150 (MQ=255)
                                                              gCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATggg                                                                                         <  2:1034339/150‑1 (MQ=255)
                                                              gCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATggg                                                                                         <  2:1822230/150‑1 (MQ=255)
                                                              gCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATggg                                                                                         <  2:272080/150‑1 (MQ=255)
                                                              gCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATggg                                                                                         >  2:778631/1‑150 (MQ=255)
                                                               cAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGGCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGc                                                                                        <  2:134236/150‑1 (MQ=255)
                                                                     tACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGc                                                                                  >  1:368035/1‑150 (MQ=255)
                                                                         tgatgaCGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCAtt                                                                              <  2:253539/150‑1 (MQ=255)
                                                                            tgaCGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGAt                                                                           >  2:577410/1‑150 (MQ=255)
                                                                             gaCGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATg                                                                          >  1:2150442/1‑150 (MQ=255)
                                                                                 cgTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGat                                                                      >  2:936376/1‑150 (MQ=255)
                                                                                 cgTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGat                                                                      <  2:1272181/150‑1 (MQ=255)
                                                                                 cgTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGat                                                                      <  1:60449/150‑1 (MQ=255)
                                                                                   ttttCCGGCTTTTGTTTTTCCTCCAGTTAGGCAAACCTTCGATACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATGCATGCCTGTTGTTCCAGCATTGCACGATGTTATTCATCACCAACATGGGCTTTCGCCATTGATGGTataa                                                                    <  2:1185964/150‑1 (MQ=255)
                                                                                   ttttCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGataa                                                                    <  1:103099/150‑1 (MQ=255)
                                                                                   ttttCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGataa                                                                    <  1:121640/150‑1 (MQ=255)
                                                                                    tttCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGataat                                                                   <  2:954615/150‑1 (MQ=255)
                                                                                       ccGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATc                                                                >  2:1711496/1‑150 (MQ=255)
                                                                                         ggCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATcaa                                                              >  2:2295500/1‑148 (MQ=255)
                                                                                         ggCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCca                                                              >  2:1281621/1‑150 (MQ=255)
                                                                                         ggCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCca                                                              >  1:1716320/1‑150 (MQ=255)
                                                                                          gCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCcag                                                             >  1:314520/1‑150 (MQ=255)
                                                                                           cTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGACACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCcagc                                                            >  2:1586939/1‑150 (MQ=255)
                                                                                            tttGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCcagca                                                           <  2:2262330/150‑1 (MQ=255)
                                                                                            tttGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCcagca                                                           >  2:1893668/1‑150 (MQ=255)
                                                                                                gTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGa                                                       <  1:2423015/150‑1 (MQ=255)
                                                                                                gTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGa                                                       >  1:1372330/1‑150 (MQ=255)
                                                                                                 tttttCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGaa                                                      >  2:155093/1‑150 (MQ=255)
                                                                                                    ttCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTc                                                   <  2:1842184/150‑1 (MQ=255)
                                                                                                     tCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTcc                                                  <  2:902243/150‑1 (MQ=255)
                                                                                                      ccaccaGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCa                                                 >  2:1469049/1‑150 (MQ=255)
                                                                                                      ccaccaGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCa                                                 <  1:811588/150‑1 (MQ=255)
                                                                                                        accaGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATg                                               <  2:531390/150‑1 (MQ=255)
                                                                                                         ccaGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGc                                              >  2:1511182/1‑150 (MQ=255)
                                                                                                           aGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCAc                                            >  1:885140/1‑150 (MQ=255)
                                                                                                            gTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGAGATAATATCCAGCAGCGAACTCCATGCAcc                                           <  2:980226/150‑1 (MQ=255)
                                                                                                               aCGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCtgt                                        >  1:1424176/1‑150 (MQ=255)
                                                                                                                  cAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCg                                     <  1:1032291/150‑1 (MQ=255)
                                                                                                                  cAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCg                                     <  2:1759225/150‑1 (MQ=255)
                                                                                                                  cAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCg                                     >  1:791776/1‑150 (MQ=255)
                                                                                                                    aaCCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGa                                   <  1:1801642/150‑1 (MQ=255)
                                                                                                                     aCCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGat                                  >  1:888576/1‑150 (MQ=255)
                                                                                                                          cGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACg                             <  1:1132087/150‑1 (MQ=255)
                                                                                                                           gCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGc                            >  1:2340163/1‑150 (MQ=255)
                                                                                                                            cTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGcc                           <  2:2236638/150‑1 (MQ=255)
                                                                                                                             tACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCt                          <  1:150846/150‑1 (MQ=255)
                                                                                                                                 gTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTc                      <  2:2150442/150‑1 (MQ=255)
                                                                                                                                   tATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCg                    <  2:1936970/150‑1 (MQ=255)
                                                                                                                                       ggCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAg                >  1:2044136/1‑150 (MQ=255)
                                                                                                                                        ggattggcCTGATGGATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGa               >  2:82260/7‑150 (MQ=255)
                                                                                                                                        gCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGa               <  1:1316261/150‑1 (MQ=255)
                                                                                                                                         cACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGAc              >  1:1195092/1‑150 (MQ=255)
                                                                                                                                         cACGCGCCTGATCCATCAGGCAAATACATGCCTGTTGTTCCAGCATTGCACGATGTTCTTAGTCACAAACATGGGCTTTCGCCATTGATGGGATAATATACAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCAAGAc              >  2:274800/1‑150 (MQ=255)
                                                                                                                                          attggcCTGATGGATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTATTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATACAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGaag             >  2:833741/5‑148 (MQ=255)
                                                                                                                                            tggcCTGATGGATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTa           >  2:1853604/3‑150 (MQ=255)
                                                                                                                                             cgcCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTaa          >  1:1849630/1‑150 (MQ=255)
                                                                                                                                              gcCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTAAt         <  1:988740/150‑1 (MQ=255)
                                                                                                                                                    tCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTAATACTTcc   <  1:18076/150‑1 (MQ=255)
                                                                                                                                                     ccATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTAATACTTCCg  <  2:2041284/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATAATTTGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTAATACTTCCG  >  NC_000913/3987025‑3987323

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