Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,107,876 A→C 5.7% I427S (ATT→AGT) ‡ speC ← ornithine decarboxylase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,107,8760AC 5.7% 277.2 / 3.0 92I427S (ATT→AGT) ‡speCornithine decarboxylase, constitutive
Reads supporting (aligned to +/- strand):  ref base A (34/52);  new base C (3/3);  total (37/55)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATATCCTTCAAAGCCGTGCCACTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTAGAAGCATGGAGCATAAAGGCGTT  >  NC_000913/3107727‑3108021
                                                                                                                                                      |                                                                                                                                                 
cATATCCTTCAAAGCCGTGCCACTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGaa                                                                                                                                                   <  2:519028/150‑1 (MQ=255)
 atatCCTTCAAAGCCGTGCCACTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAAt                                                                                                                                                  <  1:1014104/150‑1 (MQ=255)
 atatCCTTCAAAGCCGTGCCACTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAAt                                                                                                                                                  <  1:191607/150‑1 (MQ=255)
   atCCTTCAAAGCCGTGCCACTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAg                                                                                                                                                >  2:1154997/1‑150 (MQ=255)
    tCCTTCAAAGCCGTGCCACTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATACTCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACCGCTTACAGCGCGCAAGAATAGc                                                                                                                                               >  1:558866/1‑150 (MQ=255)
      cTTCAAAGCCGTGCCACTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACATTTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCTCGAATAGCCt                                                                                                                                             >  1:865585/1‑150 (MQ=255)
      cTTCAAAGCCGTGCCACTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGTGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCt                                                                                                                                             >  2:209233/1‑150 (MQ=255)
         cAAAGCCGTGCCACTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTgc                                                                                                                                          >  2:572602/1‑150 (MQ=255)
           aaGCCGTGCCACTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTgcgc                                                                                                                                        <  1:231927/150‑1 (MQ=255)
            aGCCGTGCCACTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTgcgcg                                                                                                                                       >  2:585478/1‑150 (MQ=255)
             gCCGTGCCACTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGGCTATTCTTgcgcgc                                                                                                                                      <  2:56523/150‑1 (MQ=255)
              ccGTGCCACTTCGCCCCCGGCTCAAAACTGAAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGGCTATTCTTgcgcgc                                                                                                                                      <  1:178577/150‑1 (MQ=255)
              ccGTGCCACTTCGCCCCCCGCTCAAAACTGAAAAAACGCCGGTCTCTGGCTATCACTGATGTCGGATAATCCTGCCACATTTTGCCATCATCAACTGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGGCTATTCTTGCGCGCt                                                                                                                                     <  1:33099/150‑1 (MQ=255)
               cGTGCCACTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCtt                                                                                                                                    <  1:1353630/150‑1 (MQ=255)
                gTGCCACTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTc                                                                                                                                   <  1:1078587/150‑1 (MQ=255)
                  gCCACTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCaa                                                                                                                                 >  2:1330660/1‑150 (MQ=255)
                   ccACTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAAt                                                                                                                                >  1:1643655/1‑150 (MQ=255)
                     aCTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATcc                                                                                                                              <  1:1310132/150‑1 (MQ=255)
                        tCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTa                                                                                                                           <  1:1459347/150‑1 (MQ=255)
                        tCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGACGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTa                                                                                                                           <  1:834552/150‑1 (MQ=255)
                        tCGACCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTa                                                                                                                           <  1:1662285/150‑1 (MQ=255)
                          gCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTAtc                                                                                                                         <  1:1215798/150‑1 (MQ=255)
                          gCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTAtc                                                                                                                         <  1:655252/150‑1 (MQ=255)
                             cccGGCTCAAAACTGATAAAACGTCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCa                                                                                                                      <  1:1682582/150‑1 (MQ=255)
                             cccGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCa                                                                                                                      <  1:42032/150‑1 (MQ=255)
                                ggCTCAAAAATGAAAAAACGGCGGTCGCTGGCAAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAaca                                                                                                                   <  2:629026/150‑1 (MQ=255)
                                   tCAAAACTGAAAAAACGGCGGTCGCTGCCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAAACGGCGGACCACATTACCGCGCGCATCAAAAGCCGAGCGCGCGTCACTCGCAAGCGCCacacac                                                                                                                >  1:221133/1‑150 (MQ=255)
                                    cAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACt                                                                                                               <  1:628620/150‑1 (MQ=255)
                                       aaCTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAg                                                                                                            >  1:1606536/1‑150 (MQ=255)
                                       aaCTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGACAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAg                                                                                                            <  2:1192278/150‑1 (MQ=255)
                                        aCTGAAAAAACGGCGGTCGCTGGCTAACTCTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGc                                                                                                           <  2:472960/150‑1 (MQ=255)
                                         cTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGcc                                                                                                          <  2:196094/150‑1 (MQ=255)
                                         cTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGcc                                                                                                          >  2:1696340/1‑150 (MQ=255)
                                         cTGAAAAAACGGCGGTCGCTGGCTAACACTGATCTCGGATAAACCTGCCACAATTAGCCATCCCCCACGGCCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCATGCGCTCTTCAAACCCAATCTCCCCACACACAGcc                                                                                                          >  1:230844/1‑150 (MQ=255)
                                             aaaaaaCGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCcac                                                                                                       <  2:354337/150‑1 (MQ=255)
                                               aaaaCGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAg                                                                                                     <  2:1471715/150‑1 (MQ=255)
                                                aaaCGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGc                                                                                                    >  2:661101/1‑150 (MQ=255)
                                                  aCGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCg                                                                                                  <  1:1719933/150‑1 (MQ=255)
                                                  aCGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCg                                                                                                  <  2:1321954/150‑1 (MQ=255)
                                                         tCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCa                                                                                           <  2:894018/150‑1 (MQ=255)
                                                          cGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCAc                                                                                          <  1:392470/150‑1 (MQ=255)
                                                              ggCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACtctc                                                                                      >  1:467347/1‑150 (MQ=255)
                                                                cTAACACTGATGTCGGATACTCCAGCCACAAATTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAAAAGAGTAGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGAAGCCCACTCTacc                                                                                    >  2:972729/1‑150 (MQ=255)
                                                                cTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTccc                                                                                    <  2:1263983/150‑1 (MQ=255)
                                                                 tAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTcccc                                                                                   <  2:543006/150‑1 (MQ=255)
                                                                  aaCACTGATGGCGGATAATCCTGCCACAATTTGCCCTCAACAACGTGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCATGCGCGCTTCAATCCCTATCTCAACCCACTGAGCCCACAGCCGACGCCCACTCTCCCCt                                                                                  <  1:82609/150‑1 (MQ=255)
                                                                    cacTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTc                                                                                <  2:420131/150‑1 (MQ=255)
                                                                     acTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCa                                                                               <  1:1489283/150‑1 (MQ=255)
                                                                      cTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCAt                                                                              <  2:930078/150‑1 (MQ=255)
                                                                      cTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCAt                                                                              <  2:1299326/150‑1 (MQ=255)
                                                                          tGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAAt                                                                          >  1:306062/1‑150 (MQ=255)
                                                                             cGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAAtttt                                                                       >  2:1386494/1‑150 (MQ=255)
                                                                              ggATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTg                                                                      <  2:298767/150‑1 (MQ=255)
                                                                              ggATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTg                                                                      >  1:1559265/1‑150 (MQ=255)
                                                                              ggATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTg                                                                      <  2:393563/150‑1 (MQ=255)
                                                                                  aaTCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGt                                                                  <  1:1685248/150‑1 (MQ=255)
                                                                                   aTCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGtt                                                                 <  2:441954/150‑1 (MQ=255)
                                                                                     ccTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTaa                                                               >  1:1340747/1‑150 (MQ=255)
                                                                                     ccTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTaa                                                               >  2:1272958/1‑150 (MQ=255)
                                                                                       tGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATATTGGCGTTAACa                                                             >  1:911959/1‑150 (MQ=255)
                                                                                        gCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACAt                                                            >  1:1728745/1‑150 (MQ=255)
                                                                                              aTTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCCGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTg                                                      >  1:1458795/1‑150 (MQ=255)
                                                                                                 tGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTgcag                                                   <  1:1306737/150‑1 (MQ=255)
                                                                                                   ccATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTgcagca                                                 >  1:587532/1‑150 (MQ=255)
                                                                                                      tCAACAACGGGCTGGATAAACGGGCGGAACAGCTTACTGCGCGCTAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAAGACACTCAGCCCACAGCCGACGCCCACTCTCACCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCATAc                                              <  1:384733/150‑1 (MQ=255)
                                                                                                        aacaacGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCCACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAg                                            >  1:1256933/1‑150 (MQ=255)
                                                                                                         acaacGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCGTTCATGAATTTTGGCGTGAACATCCAGTGCAGCACACAGc                                           >  2:819802/1‑150 (MQ=255)
                                                                                                          caacGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCg                                          <  1:1407608/150‑1 (MQ=255)
                                                                                                           aacGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCgg                                         <  1:1078411/150‑1 (MQ=255)
                                                                                                              gggcgggATAAACGGGCGTAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATa                                      <  1:267614/150‑1 (MQ=255)
                                                                                                               ggcgggATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAg                                     <  2:1384948/150‑1 (MQ=255)
                                                                                                               gcccgcaTAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCATCATGAATTATGGCGATAACATCCAGTGCAGCAAACAGCTCATAg                                     >  1:1106967/7‑150 (MQ=255)
                                                                                                                 cgggATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGAGAGaa                                   >  2:1208006/1‑150 (MQ=255)
                                                                                                                    gATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGATTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAgg                                <  2:98852/150‑1 (MQ=255)
                                                                                                                       aaaCGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCt                             >  2:705545/1‑150 (MQ=255)
                                                                                                                       aaaCGGGAGGAACAGCTTACAGCTCGCAAGAAAAGACTTGCGCGCTTAAATCCCTATCTCAACACTCTCAGCCGACAGCAGCCGCCCACTCTCACCTTAATGAATTTGGGCGTTAACATCCAGTGCAGCAAACAGCGGATTGTAAGGGCt                             >  2:102193/1‑150 (MQ=255)
                                                                                                                         aCGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTgg                           <  2:382520/150‑1 (MQ=255)
                                                                                                                           gggCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTa                         >  2:1250326/1‑150 (MQ=255)
                                                                                                                            ggaggAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGATTAAATCCCTATGTCAACACACTCAGACCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTAg                        <  2:1497147/147‑1 (MQ=255)
                                                                                                                               ggAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTAGAAg                     <  1:540738/150‑1 (MQ=255)
                                                                                                                               ggAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACAAACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGATTTTTGGCGTTAACTTCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTAGAAg                     <  1:1696620/150‑1 (MQ=255)
                                                                                                                                gAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTAGAAGc                    <  2:1247560/150‑1 (MQ=255)
                                                                                                                                    aGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACCCTCCGCCCACAGCCGCCGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCCGCAAACAGCGGCTAGAAAGGGCTGGTTGAAGCATgg                >  1:177290/1‑150 (MQ=255)
                                                                                                                                    aGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTAGAAGCATgg                <  2:1559265/150‑1 (MQ=255)
                                                                                                                                      cTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTAGAAGCATGGAg              <  2:306062/150‑1 (MQ=255)
                                                                                                                                        tACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGCCGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTAGAAGCATGGAGCa            >  2:262597/1‑150 (MQ=255)
                                                                                                                                         gaaGCGCGCAAGGCTATTCTTTCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACATCGGATAGAAAGGGCTGGTAGAAGAATGTAGCAt           >  2:277341/3‑150 (MQ=255)
                                                                                                                                         gaaGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGCAGCCCACTCTCCCCTTCATGAATTTAGGCGGTAACCTCCAGAGCAGCAAACAGCGGATAGAAAGGGCTGGTAGAAGCATGGAGCAt           >  1:711819/3‑150 (MQ=255)
                                                                                                                                         aCAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTAGAAGCATGGAGCAt           <  1:982870/150‑1 (MQ=255)
                                                                                                                                           aGCGCGCAAGGCTATTCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAATGGGCTGGTAGAAGCATGGTGCATaa         >  1:342069/1‑150 (MQ=255)
                                                                                                                                           aGCGCGCAAGGCTATTCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTAGAAGCATGGAGCATaa         >  2:1709464/1‑150 (MQ=255)
                                                                                                                                              gggcAAGAATTGGCTTGCGCTCTGCAATACCTATCTATAAACCCTCAGCCCCCAGCCGACGCCCACTCTCCCCTTCAGGAATTTTGGCGTTAACATCCAGTGCAGCAACCAGCGGATAGAAAGGGCTGGTAGAAGCATGGAGCATAAAgg      <  2:221133/148‑1 (MQ=255)
                                                                                                                                                  aaGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCGACATCCGACGCCCACTCTCCCCTACATGAATTTTGGCGTTAACATACAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTAGAAGCATGGAGCATAAAGGCGtt  <  1:1535467/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
CATATCCTTCAAAGCCGTGCCACTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTAGAAGCATGGAGCATAAAGGCGTT  >  NC_000913/3107727‑3108021

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: