Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,107,880 C→T 6.9% A426T (GCT→ACT) ‡ speC ← ornithine decarboxylase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,107,8800CT 6.9% 265.7 / 5.8 90A426T (GCT→ACT) ‡speCornithine decarboxylase, constitutive
Reads supporting (aligned to +/- strand):  ref base C (30/51);  new base T (3/3);  total (35/55)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

TCCTTCAAAGCCGTGCCACTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTAGAAGCATGGAGCATAAAGGCGTTATTCAAC  >  NC_000913/3107731‑3108028
                                                                                                                                                      |                                                                                                                                                    
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: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:393563/150‑1 (MQ=255)
                                                                          ggATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTg                                                                             >  1:1559265/1‑150 (MQ=255)
                                                                          ggATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTg                                                                             <  2:298767/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)
                                                                                                                                                    aGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTAGAAGCATGGAGCATAAAGGCGTTATTCaa   <  2:484922/150‑1 (MQ=255)
                                                                                                                                                     gCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTAGAAGCATGGAGCATAAAGGCGTTATTCAAc  <  1:1076524/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
TCCTTCAAAGCCGTGCCACTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTAGAAGCATGGAGCATAAAGGCGTTATTCAAC  >  NC_000913/3107731‑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: