Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 615,829 C→A 6.2% T558N (ACC→AAC)  entF → apo‑serine activating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913615,8290CA 6.2% 271.8 / 5.6 80T558N (ACC→AAC) entFapo‑serine activating enzyme
Reads supporting (aligned to +/- strand):  ref base C (41/34);  new base A (3/2);  total (44/36)
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.22e-01

GCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCG  >  NC_000913/615681‑615978
                                                                                                                                                    |                                                                                                                                                     
gCGGTGGCACTACCGCGCTCGGTCTTTTTGTCCCTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTacc                                                                                                                                                      >  2:448785/1‑150 (MQ=255)
           aCCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCa                                                                                                                                           >  1:292267/1‑150 (MQ=255)
            ccGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCaa                                                                                                                                          >  1:1510399/1‑150 (MQ=255)
                gcTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGc                                                                                                                                      <  2:278100/150‑1 (MQ=255)
                     gTCTTTTTGACCCTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCgcgc                                                                                                                                 >  2:1420771/1‑150 (MQ=255)
                       cTTTTTGACCCTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCtt                                                                                                                               >  2:1965294/1‑150 (MQ=255)
                                ccTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGt                                                                                                                      <  2:1217145/150‑1 (MQ=255)
                                   ggCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTcc                                                                                                                   >  1:713787/1‑150 (MQ=255)
                                      aCTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCaa                                                                                                                >  2:778956/1‑150 (MQ=255)
                                      aCTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCaa                                                                                                                <  2:49424/150‑1 (MQ=255)
                                       cTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAAt                                                                                                               >  2:136466/1‑150 (MQ=255)
                                        tCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCTACGATCAACTGCCTCGCTTTAGCGATGTTCCCAAtt                                                                                                              <  2:19433/150‑1 (MQ=255)
                                            tGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAAc                                                                                                          <  1:1007800/150‑1 (MQ=255)
                                            tGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAAc                                                                                                          >  2:531247/1‑150 (MQ=255)
                                              cGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACaa                                                                                                        >  2:1318976/1‑150 (MQ=255)
                                                aTAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGc                                                                                                      >  2:1836171/1‑150 (MQ=255)
                                                   gTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCtt                                                                                                   <  1:1318976/150‑1 (MQ=255)
                                                    ttGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCttt                                                                                                  >  1:1287507/1‑150 (MQ=255)
                                                      gAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGc                                                                                                <  1:1809522/150‑1 (MQ=255)
                                                      gAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGc                                                                                                <  1:899770/150‑1 (MQ=255)
                                                      gAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATAACAACCCACGAGCAAATGCCGCGATTAAGCGATGTTCCCAATTTAACAAGCCTTTga                                                                                                >  2:567159/1‑149 (MQ=255)
                                                       aaGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCt                                                                                               <  2:1463507/150‑1 (MQ=255)
                                                       aaGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCt                                                                                               >  1:405268/1‑150 (MQ=255)
                                                        aGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCta                                                                                              <  2:154037/150‑1 (MQ=255)
                                                               gCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGcc                                                                                       <  1:1980994/150‑1 (MQ=255)
                                                               gCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGcc                                                                                       >  2:1980994/1‑150 (MQ=255)
                                                                  gCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCCCCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCg                                                                                    >  1:1366643/1‑150 (MQ=255)
                                                                  gCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCg                                                                                    >  1:592232/1‑150 (MQ=255)
                                                                    cTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCt                                                                                  <  1:1589160/150‑1 (MQ=255)
                                                                    cTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCt                                                                                  >  1:1137823/1‑150 (MQ=255)
                                                                     tGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCtt                                                                                 <  1:1738656/150‑1 (MQ=255)
                                                                       gCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTac                                                                               <  2:552756/150‑1 (MQ=255)
                                                                        cTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTaca                                                                              >  2:1437334/1‑150 (MQ=255)
                                                                             gCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCa                                                                         >  1:1127999/1‑150 (MQ=255)
                                                                               tGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAAc                                                                                                          >  2:889371/1‑115 (MQ=255)
                                                                               tGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAAc                                                                                                          <  1:889371/115‑1 (MQ=255)
                                                                                 gATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGc                                                                     <  2:292267/150‑1 (MQ=255)
                                                                                 gATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGc                                                                     <  1:338791/150‑1 (MQ=255)
                                                                                  aTACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCa                                                                    >  2:1984626/1‑150 (MQ=255)
                                                                                   tACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAg                                                                   <  1:134846/150‑1 (MQ=255)
                                                                                   tACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAg                                                                   >  2:640979/1‑150 (MQ=255)
                                                                                     ccGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTg                                                                 <  1:709832/150‑1 (MQ=255)
                                                                                      cGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTgc                                                                >  2:989122/1‑150 (MQ=255)
                                                                                             ccGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTg                                                         <  2:713787/150‑1 (MQ=255)
                                                                                             ccGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTg                                                         >  2:1555545/1‑150 (MQ=255)
                                                                                              cGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGc                                                        <  2:121459/150‑1 (MQ=255)
                                                                                                  cGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCtt                                                                                 <  2:1257586/121‑1 (MQ=255)
                                                                                                  cGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCtt                                                                                 >  1:1257586/1‑121 (MQ=255)
                                                                                                   gATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCCACtt                                                   >  2:72605/1‑150 (MQ=255)
                                                                                                    aTCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACttt                                                  >  2:1487879/1‑150 (MQ=255)
                                                                                                     tCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTc                                                 <  2:446230/150‑1 (MQ=255)
                                                                                                             aaaTGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGc                                         <  1:640979/150‑1 (MQ=255)
                                                                                                                tgatgCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATc                                      >  2:1658482/1‑150 (MQ=255)
                                                                                                                   tgCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATcaca                                   >  2:1951979/1‑150 (MQ=255)
                                                                                                                     cTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGAAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACg                                 >  2:1304222/1‑150 (MQ=255)
                                                                                                                       ggAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCt                                                                                  >  2:716503/1‑99 (MQ=255)
                                                                                                                       ggAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCt                                                                                  <  1:716503/99‑1 (MQ=255)
                                                                                                                       ggAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGc                               >  1:1615584/1‑150 (MQ=255)
                                                                                                                       ggAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGc                               >  1:1196743/1‑150 (MQ=255)
                                                                                                                          aGATGTGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTta                            <  1:1474594/150‑1 (MQ=255)
                                                                                                                             tGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTatat                         >  1:668560/1‑150 (MQ=255)
                                                                                                                              gcgcGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTgc                                                                <  2:13315/110‑1 (MQ=255)
                                                                                                                              gcgcGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTgc                                                                >  1:13315/1‑110 (MQ=255)
                                                                                                                               cgcgTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATatca                       <  1:1965294/150‑1 (MQ=255)
                                                                                                                                 cgtccgtcGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATatcatc                     <  2:1137823/150‑1 (MQ=255)
                                                                                                                                 cgtccgtcGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATatcatc                     <  1:110472/150‑1 (MQ=255)
                                                                                                                                  gtccgtcGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCt                    <  1:1420771/150‑1 (MQ=255)
                                                                                                                                   tccgtcGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCtt                   >  2:96138/1‑150 (MQ=255)
                                                                                                                                   atcgtcGGTGGTAATTAACAGCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCtt                   >  2:262136/3‑150 (MQ=255)
                                                                                                                                   atcgtcGGTGGTAATTAACAGCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCtt                   >  1:1863040/3‑150 (MQ=255)
                                                                                                                                      gtcGGTGGTAATTAACAGCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTAc                >  2:424931/1‑150 (MQ=255)
                                                                                                                                           tggtAATTAACAGCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTg           <  1:262136/147‑1 (MQ=255)
                                                                                                                                            gTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTgg          >  2:251820/1‑150 (MQ=255)
                                                                                                                                             ttAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGc         >  2:1228051/1‑150 (MQ=255)
                                                                                                                                             gtAATTAACAGCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGc         <  1:424931/149‑1 (MQ=255)
                                                                                                                                              tAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCt        <  1:124563/150‑1 (MQ=255)
                                                                                                                                              tAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCt        >  2:1689435/1‑150 (MQ=255)
                                                                                                                                               aaTTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTc       >  1:106152/1‑150 (MQ=255)
                                                                                                                                                aTTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTcc      >  1:605159/1‑150 (MQ=255)
                                                                                                                                                aTTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTcc      <  2:1207142/150‑1 (MQ=255)
                                                                                                                                                aTTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTcc      <  2:1214595/150‑1 (MQ=255)
                                                                                                                                                   accaccGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCAcc   <  2:1127999/150‑1 (MQ=255)
                                                                                                                                                    ccaccGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCg  <  1:567159/150‑1 (MQ=255)
                                                                                                                                                    ccaccGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCg  >  1:1169580/1‑150 (MQ=255)
                                                                                                                                                    ccaccGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCg  >  2:1227374/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCG  >  NC_000913/615681‑615978

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