Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 816,956 C→A 100% intergenic (‑309/‑88) ybhK ← / → moaA putative CofD superfamily transferase/molybdopterin biosynthesis protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913816,9560CA97.1% 10628.0 / 137.0 3645intergenic (‑309/‑88)ybhK/moaAputative CofD superfamily transferase/molybdopterin biosynthesis protein A
Reads supporting (aligned to +/- strand):  ref base C (58/55);  new base A (1773/1755);  total (1834/1811)
Fisher's exact test for biased strand distribution p-value = 8.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 3648 total aligned reads displayed.

TTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGCGCGTAAGTTTTACTACTTGCGCCTGTCGATTA  >  NC_000913/816808‑817104
                                                                                                                                                    |                                                                                                                                                    
    cccTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAACCACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTg                                                                                                                                                  >  8:142037/1‑149 (MQ=255)
    cccTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCATCTGCACCTGGGTCAACTGATACGGTg                                                                                                                                                  >  8:348712/1‑149 (MQ=255)
      cTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATATCCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCt                                                                                                                                                >  4:309461/1‑149 (MQ=255)
       ttATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCtt                                                                                                                                               >  7:307826/1‑149 (MQ=255)
         aTTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCtt                                                                                                                                               <  1:100235/147‑1 (MQ=255)
          ttCTGTATCGAAATAGCACGATCATGTCGCTATATACATGATTACATATCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTgg                                                                                                                                            <  6:48020/149‑1 (MQ=255)
               tAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGt                                                                                                                                        >  1:17838/1‑148 (MQ=255)
               tAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCAGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGcc                                                                                                                                          <  3:329796/146‑1 (MQ=255)
                 gCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGa                                                                                                                                      <  7:109695/148‑1 (MQ=255)
                     aaTAGCACGATCATGCCGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGACTCTGCACCTGGGTAAACTGATAAGGTGCTTTGGCCGTGACAATg                                                                                                                                 >  3:236322/1‑149 (MQ=255)
                     aaTAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATg                                                                                                                                 >  5:75592/1‑149 (MQ=255)
                      aTAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGc                                                                                                                                <  7:195499/149‑1 (MQ=255)
                         gCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCg                                                                                                                             <  3:85272/149‑1 (MQ=255)
                          cACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCg                                                                                                                             >  6:228445/1‑148 (MQ=255)
                           aCGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGt                                                                                                                            >  6:176879/1‑148 (MQ=255)
                             gATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTaaa                                                                                                                         <  8:283764/149‑1 (MQ=255)
                              aTCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAg                                                                                                                        <  8:194260/149‑1 (MQ=255)
                                 aTGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGAtt                                                                                                                     <  2:360007/149‑1 (MQ=255)
                                 aTGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGAtt                                                                                                                     <  1:119401/149‑1 (MQ=255)
                                 aTGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGACGTGAAAATGCTCTTAAAGAtt                                                                                                                     >  3:364388/1‑149 (MQ=255)
                                    aCGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGcc                                                                                                                  >  3:283737/1‑149 (MQ=255)
                                    aCGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGc                                                                                                                   <  5:164891/148‑1 (MQ=255)
                                    aCGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGc                                                                                                                   <  1:146900/148‑1 (MQ=255)
                                        tatataCATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGccac                                                                                                                >  2:24126/1‑147 (MQ=255)
                                             aCATGATTACATAGCGATAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCTCTATACACTCTAGCCTCTGCCGCTGGGTCAAGTGATAAGGTGCTTTGGCCGTGGCAATGCTCGTAAAGATTGCCACCAGGGC                                                                                                          >  5:279256‑M1/1‑148 (MQ=255)
                                              cATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCG                                                                                                         >  6:149673‑M1/1‑148 (MQ=255)
                                               aTGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGTGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAA                                                                                                       <  1:318070‑M1/149‑1 (MQ=255)
                                               aTGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAA                                                                                                       <  6:362279‑M1/149‑1 (MQ=255)
                                               aTGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAA                                                                                                       <  2:163266‑M1/149‑1 (MQ=255)
                                               aTGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAA                                                                                                       <  3:188192‑M1/149‑1 (MQ=255)
                                                      cATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGgaaga                                                                                                 >  7:138263/1‑148 (MQ=255)
                                                       aTAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATGGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAg                                                                                                             <  4:229275/135‑1 (MQ=255)
                                                       aTAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGc                                                                                                                   >  1:84903/1‑129 (MQ=255)
                                                       aTAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGaaa                                                                                               <  6:127135/149‑1 (MQ=255)
                                                        tAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCCCCAGGGCGAAGGAAGAAAt                                                                                              >  8:73555/1‑149 (MQ=255)
                                                        tAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTATAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAAt                                                                                              <  5:39896/149‑1 (MQ=255)
                                                         aGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCGTTGGCCGTGACAATGCTCGGAAAGCTTGACACCCGTGCGAAGGATGAAAt                                                                                              >  5:195709/1‑148 (MQ=255)
                                                         aGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACATGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAAt                                                                                              >  2:22141/1‑148 (MQ=255)
                                                          gCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGa                                                                                            >  4:259750/1‑149 (MQ=255)
                                                            gAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGc                                                                                                                   >  2:54832/1‑124 (MQ=255)
                                                            gAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGc                                                                                                                   <  1:54832/124‑1 (MQ=255)
                                                            gAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGgaagaa                                                                                                <  6:338858/143‑1 (MQ=255)
                                                                   tgGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGtcatttcttccttcgccct                                                                                   <  7:161370‑M1/149‑20 (MQ=255)
                                                                   tgGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTc                                                                                   <  5:298289/149‑1 (MQ=255)
                                                                     gATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTccc                                                                                 >  2:273316/1‑149 (MQ=255)
                                                                      aTGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGtcatttcttccttcgccctgg                                                                                 <  7:222353‑M1/148‑22 (MQ=255)
                                                                        ggTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTa                                                                              >  8:320129/1‑149 (MQ=255)
                                                                          tAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATc                                                                            >  6:206393/1‑149 (MQ=255)
                                                                             aaaTCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGG                                                                                                           <  3:16177‑M1/115‑1 (MQ=255)
                                                                                  tcATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTAGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAgg                                                                    >  5:135702/1‑149 (MQ=255)
                                                                                    aTTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGgtg                                                                  >  5:102209/1‑149 (MQ=255)
                                                                                    aTTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGgtg                                                                  >  1:377859/1‑149 (MQ=255)
                                                                                    aTTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGgtg                                                                  >  5:33304/1‑149 (MQ=255)
                                                                                    aTTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGgt                                                                   >  7:281043/1‑148 (MQ=255)
                                                                                        acacGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTa                                                                >  4:286445/1‑147 (MQ=255)
                                                                                        acacGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTAc                                                               <  6:81244/148‑1 (MQ=255)
                                                                                          acGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATg                                                            >  2:109743/1‑149 (MQ=255)
                                                                                            gCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGGCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGTAAGAAATGACTTCGCCTCCCGTATCTGGAAAg                                                                     <  3:259899/138‑1 (MQ=255)
                                                                                              tAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCtt                                                        >  8:111876/1‑149 (MQ=255)
                                                                                              tAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCtt                                                        >  4:200511/1‑149 (MQ=255)
                                                                                                gTATCGGCATAACCACTAAACACTATAGCCTGTGCACCTGTGTCATCTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAATGATTGCCACCAGGGCGAAGGTAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATgtattc                                                       <  2:193864/148‑6 (MQ=255)
                                                                                                    cGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAac                                                  >  2:326696/1‑149 (MQ=255)
                                                                                                     ggCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAac                                                  <  4:121377/148‑1 (MQ=255)
                                                                                                     ggCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAac                                                  <  3:113674/148‑1 (MQ=255)
                                                                                                     ggCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAac                                                  <  6:283900/148‑1 (MQ=255)
                                                                                                        aTAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAactgac                                              >  4:61636/1‑149 (MQ=255)
                                                                                                            gcACTAAAGACTCTAGCCTCTGCACCTGGGTCATCTGATAAGGTGCTTTGGCCGTTACAATGCTCGTATTGATTTGGTTCAGTTAGAAGGTAGAAATGAGTTCGCCTCCCGTATCTTGAATGGTGTACATTGCTTCACAACTGACTGAt                                          <  8:166243/148‑1 (MQ=255)
                                                                                                            ccACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGAt                                          <  2:296316/149‑1 (MQ=255)
                                                                                                            ccACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGAt                                          >  3:137964/1‑149 (MQ=255)
                                                                                                             cACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATg                                         >  8:46226/1‑149 (MQ=255)
                                                                                                             cACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATg                                         >  2:383290/1‑149 (MQ=255)
                                                                                                             cACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATg                                         >  4:174390/1‑149 (MQ=255)
                                                                                                              aCTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATg                                         >  1:209587/1‑148 (MQ=255)
                                                                                                               cTAAACACTCCAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCa                                       >  4:336690/1‑149 (MQ=255)
                                                                                                                tAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAActgactga                                           >  3:306892/1‑144 (MQ=255)
                                                                                                                tAAACACTCTAGCCTCTGCAAATGGGTCAACTGATAAGGTGCTTTGGCCGTGACATGCCACGTAAAGATTGCCACCCGGGCGAAGGAAGAAATGACTTCGCCTCCCGTTTATGGAAAGGTGTACATGGCTACACACCTGACAGATGCAt                                      >  8:252236/1‑149 (MQ=255)
                                                                                                                  aaCACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCAttt                                    >  1:232724/1‑149 (MQ=255)
                                                                                                                  aaCACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCAttt                                    <  1:353912/149‑1 (MQ=255)
                                                                                                                    cacTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTgc                                  <  8:378809/149‑1 (MQ=255)
                                                                                                                    cacTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCGACCAGGGCGAAGGAAGATATGACTTCGCCTGCCGTAGCTGTAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTgc                                  <  3:21170/149‑1 (MQ=255)
                                                                                                                    cacTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTgc                                  <  8:222587/149‑1 (MQ=255)
                                                                                                                       tctAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATgg                                                           >  5:353418/1‑121 (MQ=255)
                                                                                                                         tAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGCGCGt                              >  7:334733/1‑148 (MQ=255)
                                                                                                                         tAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGAGCGt                              >  5:219744/1‑148 (MQ=255)
                                                                                                                          aGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGActt                                                                                         <  4:150430‑M2/88‑4 (MQ=255)
                                                                                                                           gCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGAt                                          <  2:88235/134‑1 (MQ=255)
                                                                                                                             ctctGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGCGCGTAAGtt                         >  1:264628/1‑149 (MQ=255)
                                                                                                                                  cACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGCGCGTAAGTTTta                      <  6:101292/147‑1 (MQ=255)
                                                                                                                                    ccTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGCGCGTAAGTTTtactac                  <  7:53138/149‑1 (MQ=255)
                                                                                                                                          tCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGCGCGTAAGTTTTACTACTTgcgc            <  7:80527/149‑1 (MQ=255)
                                                                                                                                          tCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGCGCGTAAGTTTTACTACTTgcgc            <  5:250220/149‑1 (MQ=255)
                                                                                                                                          tCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGCGCGTAAGTTTTACTACTTgcgc            <  7:385985/149‑1 (MQ=255)
                                                                                                                                          gtaaCTGATAAGGTGCTTTGGACGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACTACTGACTGATGCATTTGCGCGTAAGTTTTACTACTTgcgc            <  4:363447/147‑1 (MQ=255)
                                                                                                                                           cAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGCGCGTAAGTTTTACTACTTGCGcc           >  3:55656/1‑149 (MQ=255)
                                                                                                                                            aaCTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGCGCGTAAGTTTTACTACTTGCGCCt          <  7:361238/149‑1 (MQ=255)
                                                                                                                                             aCTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGCGCGTAAGTTTTACTACTTGCGCCTg         >  1:289685/1‑149 (MQ=255)
                                                                                                                                              cTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCAt                                      <  6:309045/119‑1 (MQ=255)
                                                                                                                                              cTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGCGCGTAAGTTTTACTACTTGCGCCTGt        >  6:254489/1‑149 (MQ=255)
                                                                                                                                                 aTAAGGTGCTTTTGCCGTGAAAATGCTCGTAAAGAATGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGCTCGTAAGTTTTACTACTTTCGCCTGTc       <  2:134692/147‑1 (MQ=255)
                                                                                                                                                 aTAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGTCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGCGCGTAAGTTTTACTACTTGCGCCTGTCGa     >  7:247605/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGCGCGTAAGTTTTACTACTTGCGCCTGTCGATTA  >  NC_000913/816808‑817104

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