Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,299,162 G→A Y139Y (TAC→TAT napH ← ferredoxin‑type protein NapH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,299,1620GA98.6% 249.4 / ‑6.8 73Y139Y (TAC→TATnapHferredoxin‑type protein NapH
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (43/29);  total (44/29)
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.55e-01

CGGTGGCCGCAACGGTAATCACACCTTTGCTACCCAGCACGCCATACAGCGCGCCAACCGGACAAATGTGCCCGCACCAGCCGTGTTCAACGACCAGTAGATCAAATAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACGTAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACGGTGTCGAATA  >  NC_000913/2298915‑2299403
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 
cGGTGGCCGCAACGGTAATCACACCTTTGCTACCCAGCACGCCATACAGCGCGCCAACCGGACAAATGTGCCCGCACCAGCCGTGTTCAACGACCAGTAGATCAAATAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAg                                                                                                                                                                                                                                                <  1:107239/251‑1 (MQ=255)
 ggTGGCCGCAACGTTAAACACACCTTTGCTACCCAGCACGCCATACAGCGCGCCAACCGGACAAATGTGCCCGCACCAGCCGTGTTCAACGACCAGTAGATCAAATAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGc                                                                                                                                                                                                                                               <  2:465366/251‑1 (MQ=255)
          aaCGGTAATCACACCTTTGCTACCCAGCACGCCATACAGCGCGCCAACCGGACAAATGTGCCCGCACCAGCCGTGTTCAACGACCAGTAGATCAAATAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGc                                                                                                                                                                                                                                      <  1:115956/251‑1 (MQ=255)
          aaCGGTAATCACACCTTTGCGACCCAGCACGCCATACAGCGCGCCAACCGGACAAATGTGCCCGCACCAGCCGTGTTCAACGACCAGTAGATCAAATAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATAAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTAAAGGCTGATCCCACCAGAATGACCACCAACAGAACATAGCGAAGGTGGc                                                                                                                                                                                                                                      >  2:317133/1‑251 (MQ=255)
              gTAATCACACCTTTACTAACCAGCACGCCATACAGCGCGCCAACAGGACACATGTGCCCGCACCAGCCGTGTTCAACGACCAGTAGATCAAATAAACACAGCGCGAGACTAAGCAGCGCGCCGCTGCCGACGCCCATTACCAGGCTACGCCCCCTCAGAGAAACGGGGATGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATATGGCgcgg                                                                                                                                                                                                                                  >  2:282023/1‑248 (MQ=255)
                         tttGCTACCCAGCACGCCATACAGCGCGCCAACCGGACAAATGTGCCCGCACCAGCCGTGTTCAACGACCAGTAGATCCAATAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGc                                                                                                                                                                                                                         >  1:86542/1‑249 (MQ=255)
                         tttGCTACCCAGCACGCCATACAGCGCGCCAACCGGACAAATGTGCCCGCACCAGCCGTGTTCAACGACCAGTAGATCAAATAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCACATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCCAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGACTGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGc                                                                                                                                                                                                                         >  1:382128/1‑249 (MQ=255)
                              tACCCAGCACGCCATACATCGCGCCAACCGGACAAATGTGCCCGCACCAGCCGTGTTCAACGACCAGTAGATCAAATAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTAACAGGCTACGCCCCATCAGAGCAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCGGATCCCACCAGAAAGACCACCAACAGCACATAGCGAATGTGGCGCGGGATCGTCGCAGACTgg                                                                                                                                                                                                                  >  2:156899/1‑250 (MQ=255)
                                              cAGCGCGCCAACCGGACAAATGTGCCCGCACCAGCCGTGTTCAACGACCAGTAGATCAAATAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGcc                                                                                                                                                                                                  <  1:300461/251‑1 (MQ=255)
                                                       aaCCGGACAAATGTGCCCGCACCAGCCGTGTTCAACGACCAGTAGATCAAATAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAAGGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAAcc                                                                                                                                                                                         <  1:232909/251‑1 (MQ=255)
                                                             aCAAATGTGCCCGCACCAGCCGTGTTCAACGACCAGTAGATCAAATAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGATCACCAGCAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCg                                                                                                                                                                                   <  1:354338/251‑1 (MQ=255)
                                                                      cccGCACCAGCCGTGTTCAACGACCAGTAGATCAAATAAAAACAGCGCGAGAATACGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGGCTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCgg                                                                                                                                                                          >  2:161382/1‑251 (MQ=255)
                                                                       ccGCACCAGCCGTGTTCAACGACCAGTAGATCAAATAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGt                                                                                                                                                                         >  1:553410/1‑251 (MQ=255)
                                                                         gCACCAGCCGTGTTCAACGACCAGTAGATCAAATAAAAACAGCGCGCGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTa                                                                                                                                                                        >  1:227256/1‑250 (MQ=255)
                                                                                   tgtTCAACGACCAGTAGATCAAATAAAAACAGCGCGAGAATAAGCAGCGCGCCGCAGCCGAAGCCCATTACCAGGACACGCCCCATCAGAGAAACGGGGCTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCAACCAGAATGACCACAACCAGCACATAGCGAATGTGGCGAGAGATCGTCGAAGACTGATTGAGGACAAAGCGCCTGCGTAACCAGCTCGCTAAGTCGGTAATAGGGTTCa                                                                                                                                                              >  2:390009/1‑250 (MQ=255)
                                                                                    gtttaaCGACCAGTAGATCAAATAAAAACAGCGCGAGAAAAAGAAGCGCGACGCTGCCGAAGTTCATTACCAGGCTACGCCCCAGCAGAAAAATGTGGTTGATACATTCCCAAAGGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCAACAACAGCACATATCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGTTAATCGGGTTCat                                                                                                                                                             <  2:457054/246‑2 (MQ=255)
                                                                                          cGACCAGTAGATCAAATAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGc                                                                                                                                                        >  1:414906/1‑249 (MQ=255)
                                                                                                gTAGATCAAATAAAAACAGCGCGAGAATAAGAAGCGCGCCGCTGCCGAAGCCCCTTACCAGGCTACGCACCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGGGAAACCAGTTCGCCAAGACGGTCATCGGGTTCAGCGGGCAGAccc                                                                                                                                                  >  2:580754/1‑249 (MQ=255)
                                                                                                 tAGATCAAATAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGACGACGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCACATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGGGGCGAGGAATCGTCGCAGACTGGTTGAGGTCAAAGCGCCTGCCTAACCAGTTCGCTAAGTCGATAATCGGGTTCAGCGGGCAGACCCAGc                                                                                                                                               >  2:395654/1‑251 (MQ=255)
                                                                                                     tCAAATAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGc                                                                                                                                            >  1:65104/1‑250 (MQ=255)
                                                                                                       aaaTAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGATGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGcca                                                                                                                                          >  1:97825/1‑248 (MQ=255)
                                                                                                            aaaaCAGAGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCAATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGATAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGCCGCAGACTGATTGAGGTCAAAACGCCTGCGTAACCAGTTAGCTAAGTCGGTAATCGGGTTAAGCGGGCAGACCCCGCTGCAAactaat                                                                                                                                    >  2:524383/1‑251 (MQ=255)
                                                                                                               aCAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGc                                                                                                                                 >  1:322853/1‑251 (MQ=255)
                                                                                                               aCAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGc                                                                                                                                 <  2:535510/251‑1 (MQ=255)
                                                                                                                      gagaATAAGCAGCGGGCCGCTGCCGAAGCCCATTACCAGGATACGCCCCATCAGAGAAACGGGGTTGGTCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCCGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGACGGTAATCGGGTTCAGCGGGCAGCCCCAGCTGCAAAATAATCGCTTTccct                                                                                                                          >  2:358164/1‑250 (MQ=255)
                                                                                                                       agaATAAGCAGCGCGCCGCTGCCGATGCCCAGTACCAGGCTACGCCCCATCAGAGAAACGGGGTTAATCCATTCCCAAATGAGCGGGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTACGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGc                                                                                                                         <  2:39032/251‑1 (MQ=255)
                                                                                                                         aaTAAGCAGCGCTCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGGCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGcc                                                                                                                        >  2:438694/1‑250 (MQ=255)
                                                                                                                         aaTAAGCAGCGCGCCGCTGCCGAAGCCCATTCCCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCa                                                                                                                       <  1:329647/251‑1 (MQ=255)
                                                                                                                          aTAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTACGTCGGGAATCGGGTTCAGCGGGCAGACCCAGCTGCCAAATAATCGCTTTCCCGcc                                                                                                                        >  2:265280/1‑249 (MQ=255)
                                                                                                                           taatcagcGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGAGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGAGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCACAATAATCGCTTTCCCACCAgg                                                                                                                     >  2:116362/5‑251 (MQ=255)
                                                                                                                               cagcGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGt                                                                                                                 >  1:206185/1‑251 (MQ=255)
                                                                                                                                agcGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGt                                                                                                                 >  1:481424/1‑250 (MQ=255)
                                                                                                                                 gggcgcCGCTGCCGAAGCCCATTCCCAGGCTACGCCCCAACAGAGAAACGGGGTTGATCCATTACCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTag                                                                                                               <  2:156009/249‑1 (MQ=255)
                                                                                                                                   gcgcCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTag                                                                                                               >  1:20669/1‑249 (MQ=255)
                                                                                                                                   gcgcCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTag                                                                                                               >  1:354012/1‑249 (MQ=255)
                                                                                                                                          tGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAg                                                                                                                                                                                                                                                >  1:587500/1‑113 (MQ=255)
                                                                                                                                           gCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTg                                                                                                      <  2:60526/250‑1 (MQ=255)
                                                                                                                                                aGCCCATTACAAGGCTACGCACCATCAGAGAAACGGAGTTGATGAATTCCGAAATGAGCGTGCCAGTTAACACTGATCCCACCAGAATGACCACCAACAGCACATATCGAATGTGGCGAGGGATCGGAGCAGGCTGATTTATGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGGCGGTAATCGGGTTCAGCATGCAGACCCCTGTGCGAAATCCACTGTTTCCCGCCAGGGCGTGGAGCACGGCGATAATc                                                                                                >  2:1222/1‑251 (MQ=255)
                                                                                                                                                   ccATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATcacc                                                                                             <  1:127550/251‑1 (MQ=255)
                                                                                                                                                       tACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATcaccgca                                                                                          <  1:317133/250‑1 (MQ=255)
                                                                                                                                                        aCCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAAGGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATcaccgcac                                                                                         <  1:317890/250‑1 (MQ=255)
                                                                                                                                                          cAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCgg                                                                                      <  2:65104/251‑1 (MQ=255)
                                                                                                                                                            ggCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCAACCAGAATGCCCACCAACAGCACATACCGAATGTGGCGCGGGCTCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACACAGCTGCAAAATAATCTCTTTCACGCCAGGACGCAGGACACGGTGGTAATCACCGCACCGGTc                                                                                    >  2:241088/1‑251 (MQ=255)
                                                                                                                                                                   ccccATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCAACGcc                                                                              >  1:256135/1‑250 (MQ=255)
                                                                                                                                                                     ccATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAAGGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCGAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCAACGCCAcc                                                                           <  1:15833/251‑1 (MQ=255)
                                                                                                                                                                                  ggggTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCAACGCCACCGTggcaggc                                                                  >  1:230435/1‑247 (MQ=255)
                                                                                                                                                                                    ggtggacccATTCCCAAACGAGCGTGCCATTGAAGGCGGATCCCCCCAGAATGACCACCAACAGCACATATCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTTCCGCTAGGGCGTAGAGCACGGTGATAATCACCGCACCGGGCAACGCCACCGTGGCAGGCAGATGa                                                            <  2:49874/244‑1 (MQ=255)
                                                                                                                                                                                    ggTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGATCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGATGCAAAATAATCGCTTTCCCGCCAGGGAGTAGAGCACGGTGATAATCACCGCAACGGTCAACGCCACCGTGGCAGGCAGATg                                                             >  1:8122/1‑250 (MQ=255)
                                                                                                                                                                                       tGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCAACGCCACCGTGGCAGGCAGATGAcc                                                          >  1:587510/1‑250 (MQ=255)
                                                                                                                                                                                       tGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCATGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCg                                                                                                                                  >  1:568533/1‑178 (MQ=255)
                                                                                                                                                                                       tGATCCAGTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAATCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCg                                                                                                                                  <  2:568533/178‑1 (MQ=255)
                                                                                                                                                                                        gATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGACGCAGACCGATTGAGGTCAAAGCGCCTGCGTCACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCAACGCCACCGTGGCAGGCAGATGACCAc                                                        >  1:399204/1‑251 (MQ=255)
                                                                                                                                                                                           ccATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCTTCCATAAGGACCACCAACAGCACATAGCGAATGTGGCGAGGGAGCGTCGCAGACTTATTGAGGTCAAAGCGCCTGCGTAACCAGTGCGCTAGGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAAGAATCGCTTGCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCAACGCCACCGTGGCAGGCAGATGACCACTgg                                                     <  1:573369/251‑1 (MQ=255)
                                                                                                                                                                                            cATTCCCAAATGAGCGTTCCAGTTAAGGCTGATCCCTCCATAAGGCCCACCATCAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGGAACCAGTTCGCTAAGTCGGTAAGCGGGTTCAGCTGGCAGACCCAGCTGCAAAAGAATCGCTTGCCCGCCAGGGCGTAGAGCACGGTGATAAGCACCGCACCGTTCAACGCCACCGGGGCAGTCAGATGAccgctagc                                                    <  2:97825/251‑7 (MQ=255)
                                                                                                                                                                                                  cAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCt                                                                                                                                              >  1:577532/1‑155 (MQ=255)
                                                                                                                                                                                                  cAAATGAGCGTGCCAGTGAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGGCTGAGTGAGGGCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACACAGCt                                                                                                                                              <  2:577532/155‑1 (MQ=255)
                                                                                                                                                                                                       gAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGGCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGc                                          <  2:108154/250‑1 (MQ=255)
                                                                                                                                                                                                          cGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCCGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCCACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCg                                      >  2:208275/1‑251 (MQ=255)
                                                                                                                                                                                                           gTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGt                                     >  1:548061/1‑251 (MQ=255)
                                                                                                                                                                                                            gtccAGTTAAGGCCGATCCCACCAGAAAGACCACCAACAGCCAATAGCGAAAGTGGCGAGGGATCTTCGCAGACTGAGTTAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGGCGGTAATCGGGTTCAGCGGGCAGCCCCAGCTGCAAAATAACCGTTTGCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCAACGCCACCGTGGCAGGCAGATGACCCCTGGCAAGACTTGGCAGCGtc                                    <  2:92277/249‑1 (MQ=255)
                                                                                                                                                                                                              ccAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGtcat                                  >  1:159805/1‑251 (MQ=255)
                                                                                                                                                                                                               cAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCACGACTTTGCAGCGtcatc                                 >  1:125687/1‑251 (MQ=255)
                                                                                                                                                                                                                 gTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCACGCCTTTGCAGCGtcatc                                 >  1:382950/1‑249 (MQ=255)
                                                                                                                                                                                                                    aaGGCTGATCCCACCAGAATGACCACCAACAGCACATGGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGGTCAGCGGGCAGACGCATCTGCAAAATAATCCCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCTCCGCACCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGGCTTTGCAGCGTCATCAGCg                             <  2:12521/250‑1 (MQ=255)
                                                                                                                                                                                                                    aaGGCTGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCgg                            >  1:542560/1‑251 (MQ=255)
                                                                                                                                                                                                                         tGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCgg                       >  1:54012/1‑251 (MQ=255)
                                                                                                                                                                                                                         tGATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGGGATAATCACCGCACCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCgg                       >  1:570455/1‑251 (MQ=255)
                                                                                                                                                                                                                          gATCCCACCAGAATGACCACCAATAGCACATAGCGAATGTGGCGAAGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGt                      >  1:287583/1‑251 (MQ=255)
                                                                                                                                                                                                                          gATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCAACGCCACCGTGTCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCgg                       <  1:224510/250‑1 (MQ=255)
                                                                                                                                                                                                                          gATCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCg                        <  1:438694/249‑1 (MQ=255)
                                                                                                                                                                                                                           aTCCCACCAGAATGACCACCAACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTc                     <  1:263135/251‑1 (MQ=255)
                                                                                                                                                                                                                                     aaTGACCACCACCAGCACGTAGCGAATGTGGCGAGGGATCGTCGCAGACTGCTTGAGGTCAAAGCGCCTGCGTAACCAGATCGCTACGTCGGTAATCGGGTTAAGCGGGCAGACCCAGCTGCAAACTAATCGCTTTCCCGCAAGGGCGTAGAGCACGGGGATAATCACCGCACCGGTCAACGCCACCGTGGCAGGCAGACGCCCACTGACCAGCCTTTGCAGCATCATCAGCGGATCGGTCAGCGGAACgg           >  2:77274/1‑251 (MQ=255)
                                                                                                                                                                                                                                              caACAGCACATAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCATGGTGATAATCACCGCACCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACGGTGTCGAATa  <  1:145245/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 
CGGTGGCCGCAACGGTAATCACACCTTTGCTACCCAGCACGCCATACAGCGCGCCAACCGGACAAATGTGCCCGCACCAGCCGTGTTCAACGACCAGTAGATCAAATAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCATTACCAGGCTACGCCCCATCAGAGAAACGGGGTTGATCCATTCCCAAATGAGCGTGCCAGTTAAGGCTGATCCCACCAGAATGACCACCAACAGCACGTAGCGAATGTGGCGAGGGATCGTCGCAGACTGATTGAGGTCAAAGCGCCTGCGTAACCAGTTCGCTAAGTCGGTAATCGGGTTCAGCGGGCAGACCCAGCTGCAAAATAATCGCTTTCCCGCCAGGGCGTAGAGCACGGTGATAATCACCGCACCGGTCAACGCCACCGTGGCAGGCAGATGACCACTGGCAAGACTTTGCAGCGTCATCAGCGGATCGGTCAGCGGAACGGTGTCGAATA  >  NC_000913/2298915‑2299403

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