Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,107,879 G→T 6.7% A426D (GCT→GAT) ‡ speC ← ornithine decarboxylase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,107,8790GT 6.7% 264.7 / 7.2 90A426D (GCT→GAT) ‡speCornithine decarboxylase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (33/51);  new base T (3/3);  total (36/54)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

ATCCTTCAAAGCCGTGCCACTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTAGAAGCATGGAGCATAAAGGCGTTATTCAAC  >  NC_000913/3107730‑3108028
                                                                                                                                                      |                                                                                                                                                     
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:655252/150‑1 (MQ=255)
                       gCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTAtc                                                                                                                                <  1:1215798/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:1696340/1‑150 (MQ=255)
                                      cTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGcc                                                                                                                 <  2:196094/150‑1 (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:1299326/150‑1 (MQ=255)
                                                                   cTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCAt                                                                                     <  2:930078/150‑1 (MQ=255)
                                                                       tGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAAt                                                                                 >  1:306062/1‑150 (MQ=255)
                                                                          cGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAAtttt                                                                              >  2:1386494/1‑150 (MQ=255)
                                                                           ggATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTg                                                                             >  1:1559265/1‑150 (MQ=255)
                                                                           ggATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTg                                                                             <  2:393563/150‑1 (MQ=255)
                                                                           ggATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTg                                                                             <  2:298767/150‑1 (MQ=255)
                                                                               aaTCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGt                                                                         <  1:1685248/150‑1 (MQ=255)
                                                                                aTCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGtt                                                                        <  2:441954/150‑1 (MQ=255)
                                                                                  ccTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTaa                                                                      >  2:1272958/1‑150 (MQ=255)
                                                                                  ccTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTaa                                                                      >  1:1340747/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)
                                                                                                                                                     aGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTAGAAGCATGGAGCATAAAGGCGTTATTCaa   <  2:484922/150‑1 (MQ=255)
                                                                                                                                                      gCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTAGAAGCATGGAGCATAAAGGCGTTATTCAAc  <  1:1076524/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
ATCCTTCAAAGCCGTGCCACTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTAGAAGCATGGAGCATAAAGGCGTTATTCAAC  >  NC_000913/3107730‑3108028

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: