Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,305,086 C→G 11.5% intergenic (+693/+22) eco → / ← mqo ecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,305,0860CG11.5% 286.4 / 25.9 87intergenic (+693/+22)eco/mqoecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base C (39/38);  new base G (4/6);  total (43/44)
Fisher's exact test for biased strand distribution p-value = 7.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

GGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTACTGC  >  NC_000913/2304938‑2305228
                                                                                                                                                    |                                                                                                                                              
ggCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATAc                                                                                                                                                <  1:942842/149‑1 (MQ=255)
     gTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGc                                                                                                                                           >  1:855759/1‑149 (MQ=255)
        ttATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTATTATCAGCATACGCCACATCCGGCATACCATGCCgg                                                                                                                                        <  3:910427/149‑1 (MQ=255)
         tATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGa                                                                                                                                       <  3:199088/149‑1 (MQ=255)
         tATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGa                                                                                                                                       <  4:477302/149‑1 (MQ=255)
          aTCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCTGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGAt                                                                                                                                      <  4:692347/149‑1 (MQ=255)
          aTCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATGGTATGCCGGAt                                                                                                                                      <  3:1025319/149‑1 (MQ=255)
          aTCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATGGTATGCCGGAt                                                                                                                                      <  4:353930/149‑1 (MQ=255)
          aTCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGAt                                                                                                                                      >  2:798004/1‑149 (MQ=255)
          aTCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGAt                                                                                                                                      >  1:1068090/1‑149 (MQ=255)
          aTCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGAt                                                                                                                                      <  3:853007/149‑1 (MQ=255)
           tCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGAt                                                                                                                                      <  4:999251/148‑1 (MQ=255)
           tCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGAt                                                                                                                                      >  3:999251/1‑148 (MQ=255)
             aGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATGGTATGCCGGAtgtg                                                                                                                                   <  1:75303/149‑1 (MQ=255)
             aGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGAtgtg                                                                                                                                   <  2:292680/149‑1 (MQ=255)
             aGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGAtgtg                                                                                                                                   <  3:470234/149‑1 (MQ=255)
              ggCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATGGTATGCCGGATGTgg                                                                                                                                  <  1:148018/149‑1 (MQ=255)
              ggCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTgg                                                                                                                                  <  1:1097635/149‑1 (MQ=255)
               gCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGc                                                                                                                                 <  1:962940/149‑1 (MQ=255)
                 cTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGc                                                                                                                                           <  4:128072/137‑1 (MQ=255)
                 cTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGc                                                                                                                                           >  3:128072/1‑137 (MQ=255)
                   acacCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTAt                                                                                                                             <  2:229164/149‑1 (MQ=255)
                               aGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCa                                                                                                                 >  4:399562/1‑149 (MQ=255)
                                   cTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAatat                                                                                                             >  4:650405/1‑149 (MQ=255)
                                    tCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATGGTATGCCGGAt                                                                                                                                      >  2:471066/1‑123 (MQ=255)
                                    tCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATGGTATGCCGGAt                                                                                                                                      <  1:471066/123‑1 (MQ=255)
                                    tCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATc                                                                                                            >  1:150097/1‑149 (MQ=255)
                                    tCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATc                                                                                                            <  1:869793/149‑1 (MQ=255)
                                     ccACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATcc                                                                                                           <  4:118962/149‑1 (MQ=255)
                                      cACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCa                                                                                                                 <  4:810611/142‑1 (MQ=255)
                                      cACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCa                                                                                                                 >  3:810611/1‑142 (MQ=255)
                                      cACCCCGTAGGTTGGATAAGGCGGTTACACCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCg                                                                                                          >  1:26104/1‑149 (MQ=255)
                                          ccGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCAc                                                                                                      >  1:530353/1‑149 (MQ=255)
                                            gTAGGTTGGATAAGGCGTTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATAc                                                                                                                                                <  4:913648/105‑1 (MQ=255)
                                            gTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATAc                                                                                                                                                >  3:913648/1‑105 (MQ=255)
                                               ggTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCActtctt                                                                                                 <  1:1048516/149‑1 (MQ=255)
                                                gTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGCGCCGGTTGCCTGATGCTACGCTGGCGCGTCGTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCttt                                                                                                <  2:530353/149‑1 (MQ=255)
                                                  tGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCttttt                                                                                              >  4:91461/1‑149 (MQ=255)
                                                     aTAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATGGTATGCCGGAtgtg                                                                                                                                   <  3:106164/109‑1 (MQ=18)
                                                     aTAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATGGTATGCCGGAtgtg                                                                                                                                   >  4:106164/1‑109 (MQ=18)
                                                     aTAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGaa                                                                                           <  4:997549/149‑1 (MQ=255)
                                                                  gccgcATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGttt                                                                              <  3:291397/149‑1 (MQ=255)
                                                                          cAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAAc                                                                                                                    >  2:88179/1‑103 (MQ=255)
                                                                          cAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAAc                                                                                                                    <  1:88179/103‑1 (MQ=255)
                                                                           aaCATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCgg                                                                     >  4:302018/1‑149 (MQ=255)
                                                                              aTCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGttt                                                                  >  4:81785/1‑149 (MQ=255)
                                                                                 tAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCATATCCGACATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTGGAACGGGTTGCGGTTTAAACTGCGGTTTCgg                                                               >  3:668659/1‑149 (MQ=255)
                                                                                  aaCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGc                                                              >  4:895528/1‑149 (MQ=255)
                                                                                     gCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTa                                                           >  3:134026/1‑149 (MQ=255)
                                                                                       ccGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCg                                                             <  2:870281/145‑1 (MQ=255)
                                                                                       ccGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCg                                                             >  1:870281/1‑145 (MQ=255)
                                                                                         gAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTAt                                                       <  3:91461/149‑1 (MQ=255)
                                                                                           gCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATct                                                     <  3:183643/149‑1 (MQ=255)
                                                                                              ggTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGaa                                                                                           >  2:606006/1‑108 (MQ=255)
                                                                                              ggTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGaa                                                                                           <  1:606006/108‑1 (MQ=255)
                                                                                                    cTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCTACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGctgcctgc                                            >  1:787222/1‑149 (MQ=255)
                                                                                                    cTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGctgcccgc                                            <  1:798004/149‑1 (MQ=255)
                                                                                                         gCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCtt                                       >  1:634278/1‑149 (MQ=255)
                                                                                                           gACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGt                                     <  2:1068090/149‑1 (MQ=255)
                                                                                                           gACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGt                                     >  4:276395/1‑149 (MQ=255)
                                                                                                               cTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTa                                 <  4:1017347/149‑1 (MQ=255)
                                                                                                                tGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAg                                <  3:681198/149‑1 (MQ=255)
                                                                                                                       tCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAAcc                         >  2:903478/1‑149 (MQ=255)
                                                                                                                          tATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCa                      <  4:958090/149‑1 (MQ=255)
                                                                                                                          tATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCa                      >  4:504329/1‑149 (MQ=255)
                                                                                                                          tATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCa                      <  2:847164/149‑1 (MQ=255)
                                                                                                                           aTCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAg                     >  1:1067534/1‑149 (MQ=255)
                                                                                                                              aGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCAc                  >  3:543671/1‑149 (MQ=255)
                                                                                                                               gCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACt                 >  2:1087814/1‑149 (MQ=255)
                                                                                                                                cATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACtt                >  4:1060178/1‑149 (MQ=255)
                                                                                                                                 aTACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTc               >  2:429767/1‑149 (MQ=255)
                                                                                                                                  tACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCg              >  1:963818/1‑149 (MQ=255)
                                                                                                                                  tACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCg              <  1:518840/149‑1 (MQ=255)
                                                                                                                                   aCGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGcgaagtgc             <  2:150097/149‑8 (MQ=255)
                                                                                                                                     gCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTg           >  4:423028/1‑149 (MQ=255)
                                                                                                                                       cacaTCCGGCATACCATGCCGGATGTGGCGTATCTTTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGgt         <  1:903478/149‑1 (MQ=255)
                                                                                                                                       cacaTCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGgt         >  3:311332/1‑149 (MQ=255)
                                                                                                                                        acaTCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAg                                <  3:803877/125‑1 (MQ=255)
                                                                                                                                        acaTCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAg                                >  4:803877/1‑125 (MQ=255)
                                                                                                                                          aTCCGGCATGGTATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGTTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTa      >  2:1023842/1‑149 (MQ=255)
                                                                                                                                          aTCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTa      >  2:955273/1‑149 (MQ=255)
                                                                                                                                           tCCGGCATGGTATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTAc     >  1:883749/1‑149 (MQ=255)
                                                                                                                                           tCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTAc     >  4:4112/1‑149 (MQ=255)
                                                                                                                                              ggCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCg                                                                      <  1:194564/81‑1 (MQ=255)
                                                                                                                                              ggCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCg                                                                      >  2:194564/1‑81 (MQ=255)
                                                                                                                                              ggCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCgg                                                               <  2:400141/88‑1 (MQ=255)
                                                                                                                                              ggCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCgg                                                               >  1:400141/1‑88 (MQ=255)
                                                                                                                                              ggCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTACTGc  <  2:963818/149‑1 (MQ=255)
                                                                                                                                                   aCCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGtt                              >  2:1011782/1‑116 (MQ=255)
                                                                                                                                                   aCCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGtt                              <  1:1011782/116‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
GGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTACTGC  >  NC_000913/2304938‑2305228

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