Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,085,902 C→T 5.3% D4N (GAC→AAC)  speA ← biosynthetic arginine decarboxylase, PLP‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,085,9020CT 5.3% 245.8 / 2.3 75D4N (GAC→AAC) speAbiosynthetic arginine decarboxylase, PLP‑binding
Reads supporting (aligned to +/- strand):  ref base C (36/35);  new base T (2/2);  total (38/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

CGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGCACATGGATTTTCCAGCAGTGAATGCTGACGCTCCAACTGCGCAACCGGTTTCTTTTT  >  NC_000913/3085761‑3086048
                                                                                                                                             |                                                                                                                                                  
cGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGATTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCTCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGa                                                                                                                                              <  5:349966/148‑1 (MQ=255)
 gTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGAc                                                                                                                                             <  4:79183/148‑1 (MQ=255)
      agtagtTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGc                                                                                                                                        >  3:274727/1‑148 (MQ=255)
        tagtTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTTGTCAGACATAGCg                                                                                                                                       >  8:213595/1‑147 (MQ=255)
        tagtTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTTGTCAGACATAGCg                                                                                                                                       >  8:208253/1‑147 (MQ=255)
          gtTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAAc                                                                                                                                    >  6:304376/1‑148 (MQ=255)
               tGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCa                                                                                                                                >  4:291490/1‑147 (MQ=255)
                    caccaGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAattat                                                                                                                          <  6:141222/148‑1 (MQ=255)
                      ccaGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAAttattt                                                                                                                        <  8:150467/148‑1 (MQ=255)
                       caGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAAttattt                                                                                                                        >  4:70851/1‑147 (MQ=255)
                        aGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAAttatttta                                                                                                                      >  1:325090/1‑148 (MQ=255)
                          gCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATgtcgtc                                                                                                                                                 <  6:134305/119‑1 (MQ=255)
                          gCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATgtcgtc                                                                                                                                                 >  5:134305/1‑119 (MQ=255)
                          gCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAAttattttatt                                                                                                                    <  3:70851/148‑1 (MQ=255)
                          gCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAAttattttat                                                                                                                     >  4:235498/1‑147 (MQ=255)
                           cAATATTGTAAGTACGCAGCATCTTGCTGGCTTGCTGGTAGCTCATTGCAAACTACTGCAGGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTttatta                                                                                                                   <  6:232050/148‑1 (MQ=255)
                           cAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAAttattttatt                                                                                                                    >  8:210263/1‑147 (MQ=255)
                           cAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTttatta                                                                                                                   <  6:181430/148‑1 (MQ=255)
                           cAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTttatta                                                                                                                   <  2:73258/148‑1 (MQ=255)
                           cAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTttatta                                                                                                                   <  2:78547/148‑1 (MQ=255)
                           cAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTttatta                                                                                                                   <  4:155118/148‑1 (MQ=255)
                                 tGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAgtgt                                                                                                              >  2:216401/1‑147 (MQ=255)
                                 tGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAgtgt                                                                                                              >  2:210653/1‑147 (MQ=255)
                                   tAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTaaa                                                                                                           <  1:356301/148‑1 (MQ=255)
                                        aCGCAGCATCTTGCTGGCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGt                                                                                                      >  5:95076/1‑148 (MQ=255)
                                         cGCAGCATCTTGCTGGCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGt                                                                                                      <  4:274727/147‑1 (MQ=255)
                                          gcagcaTCTTGCTGGCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTa                                                                                                    <  3:291490/148‑1 (MQ=255)
                                             gcaTCTTGCTGGCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAAc                                                                                                  >  3:328107/1‑147 (MQ=255)
                                                 cTTGCTGGCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTa                                                                                             >  1:177214/1‑148 (MQ=255)
                                                 cTTGCTGGCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTa                                                                                             >  1:78266/1‑148 (MQ=255)
                                                       ggCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTaa                                                                                                   >  1:125944/1‑136 (MQ=255)
                                                       ggCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTaa                                                                                                   <  2:125944/136‑1 (MQ=255)
                                                       ggCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACCAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAg                                                                                       >  4:56345/1‑148 (MQ=255)
                                                                   gCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCg                                                                           <  3:202551/148‑1 (MQ=255)
                                                                     tCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAAttatttt                                                                                                                       <  7:130161/102‑1 (MQ=255)
                                                                     tCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAAttatttt                                                                                                                       >  8:130161/1‑102 (MQ=255)
                                                                     tCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGa                                                                          >  7:78304/1‑147 (MQ=255)
                                                                             aCCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAAcc                                                                  >  6:161587/1‑147 (MQ=255)
                                                                             aCCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAAcc                                                                  >  1:277416/1‑147 (MQ=255)
                                                                                  ctGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATaa                                                             <  6:360644/147‑1 (MQ=255)
                                                                                    gCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGt                                                                                                      >  7:251829/1‑104 (MQ=255)
                                                                                    gCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGt                                                                                                      <  8:251829/104‑1 (MQ=255)
                                                                                     cATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGcac                                                         <  5:161587/148‑1 (MQ=255)
                                                                                     cATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGc                                                           <  2:173114/146‑1 (MQ=255)
                                                                                      aTGGAGCGTACTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCAGGAACTGCGATAACAACCCATAAGcaca                                                        >  4:96784/1‑148 (MQ=255)
                                                                                        ggAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGa                                                                                                                                      <  5:65536/68‑1 (MQ=255)
                                                                                        ggAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGa                                                                                                                                      >  6:65536/1‑68 (MQ=255)
                                                                                        ggAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATaacaa                                                                    <  8:173838/134‑1 (MQ=255)
                                                                                        ggAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATaacaa                                                                    >  7:173838/1‑134 (MQ=255)
                                                                                            cGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGCACATGGAtt                                                  <  7:34502/148‑1 (MQ=255)
                                                                                             gtagtaCACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGCACATGGAttt                                                 >  6:183088/1‑148 (MQ=255)
                                                                                              tagtaCACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGCACATGGAtttt                                                >  8:254624/1‑148 (MQ=255)
                                                                                               agtaCACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGCACATGGATTTTc                                               >  7:167448/1‑148 (MQ=255)
                                                                                                  acacCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGCACATGGATTTTCcag                                            >  6:260293/1‑148 (MQ=255)
                                                                                                  acacCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGCACATGGATTTTCcag                                            >  3:21114/1‑148 (MQ=255)
                                                                                                   cacCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGCACATGGATTTTCcagc                                           <  2:78266/148‑1 (MQ=255)
                                                                                                   cacCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGCACATGGATTTTCcagc                                           <  8:148574/148‑1 (MQ=255)
                                                                                                   cacCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGCACATGGATTTTCcagc                                           <  7:210263/148‑1 (MQ=255)
                                                                                                     ccGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGca                                                          <  1:301359/131‑1 (MQ=255)
                                                                                                     ccGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGca                                                          >  2:301359/1‑131 (MQ=255)
                                                                                                     ccGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGCACATGGATTTTCcagcag                                         >  8:47431/1‑148 (MQ=255)
                                                                                                        tgtTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGCACATGGATTTTCCAGCAGTGa                                      <  3:360763/148‑1 (MQ=255)
                                                                                                          ttCGCCCGCTGACGAAGGCAAACCCATAGACATGTTGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGCACATGGATTTTCCAGCAGTGaa                                     <  7:208253/147‑1 (MQ=255)
                                                                                                          ttCGCCCGCTGACGAAGGCAAACCCATAGACATGTTGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGCACATGGATTTTCCAGCAGTGaa                                     <  7:213595/147‑1 (MQ=255)
                                                                                                          ttCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGCACATGGATTTTCCAGCAGTGAAt                                    <  2:277416/148‑1 (MQ=255)
                                                                                                          ttCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGCACATGGATTTTCCAGCAGTGAAt                                    <  4:328107/148‑1 (MQ=255)
                                                                                                             gCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGCACATGGATTTTCCAGCAGTGAATGCt                                 <  8:78304/148‑1 (MQ=255)
                                                                                                              cccGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGCACATGGATTTTCCAGCAGTGAATGCTg                                <  7:47431/148‑1 (MQ=255)
                                                                                                                        aaGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTTGAACTGCGATAACAACCCATAAGCACATGGATTTTCCAGCAGTGAATGCTGACGCTCCAAc                      <  3:241495/148‑1 (MQ=255)
                                                                                                                                cccATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGCACATGGATTTTCCAGCAGTGAATGCTGACGCTCCAACTGCGCAAc              >  8:207217/1‑148 (MQ=255)
                                                                                                                                    tAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAAACCATAAGCACATGGATTTTCCAGCAGTGAATGCTGACGCTCCAACTGCGCAACCGGt          >  4:17956/1‑148 (MQ=255)
                                                                                                                                     aGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGCACATGGATTTTCCAGCAGTGAATGCTGACGCTCCa                        >  1:254519/1‑133 (MQ=255)
                                                                                                                                     aGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGCACATGGATTTTCCAGCAGTGAATGCTGACGCTCCa                        <  2:254519/133‑1 (MQ=255)
                                                                                                                                     aGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGCACATGGATTTTCCAGCAGTGAATGCTGACGCTCCAACTGCGCAACCGGtt         >  4:236728/1‑148 (MQ=255)
                                                                                                                                       aCATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGCACATGGATTTTCCAGCAGTGAATGCTGACGCTCCAACTGCGCAACCGGTTTc       >  6:152001/1‑148 (MQ=255)
                                                                                                                                         aTGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGCACATGGTTTTTGCAGCAGTGAATGCTGACGCTCCAACTGCGCAACCGGtttttt     <  3:96784/148‑1 (MQ=255)
                                                                                                                                          tGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGCACATGGATTTTCCAGCAGTGAATGCTGACGCTCCAACTGCGCAACCGGTTTCttt    >  6:256365/1‑148 (MQ=255)
                                                                                                                                           gtcgtcAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGCACATGGATTTTcc                                              >  8:12580/1‑105 (MQ=255)
                                                                                                                                           gtcgtcAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGCACATGGATTTTcc                                              <  7:12580/105‑1 (MQ=255)
                                                                                                                                           gtcgtcAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGCACATGGATTTTCCAGCAGTGAATGCTGACGCTCCAACTGCGCAACCGGTTTCtttt   >  2:295518/1‑148 (MQ=255)
                                                                                                                                           gtcgtcAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGCACATGGATTTTCCAGCAGTGAATGCTGACGCTCCAACTGCGCAACCGGTTTCtttt   >  7:37860/1‑148 (MQ=255)
                                                                                                                                            tcgtcAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGCACATGGATTTTCCAGCAGTGAATGCTGACGCTCCAACTGCGCAACCGGTTTCttttt  <  1:278615/148‑1 (MQ=255)
                                                                                                                                            tcgtcAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGCACATGGATTTTCCAGCAGTGAATGCTGACGCTCCAACTGCGCAACCGGTTTCttttt  <  2:177214/148‑1 (MQ=255)
                                                                                                                                            tcgtcAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGCACATGGATTTTCCAGCAGTGAATGCTGACGCTCCAACTGCGCAACCGGTTTCttttt  >  2:85611/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
CGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTATCGCAGCCTGGAACTGCGATAACAACCCATAAGCACATGGATTTTCCAGCAGTGAATGCTGACGCTCCAACTGCGCAACCGGTTTCTTTTT  >  NC_000913/3085761‑3086048

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