Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 547,891 C→T L404L (CTC→CTT fdrA → putative oxamate carbamoyltransferase FdrA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913547,8910CT94.9% 114.7 / ‑5.3 40L404L (CTC→CTTfdrAputative oxamate carbamoyltransferase FdrA
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (23/14);  minor base G (0/2);  total (24/16)
Fisher's exact test for biased strand distribution p-value = 1.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATAGCGCCTGTCAGCAGCGGATTTATTTGCGGTTTGTATACCGGCGGTACGCTGGCTGCCGAAGCGGCGGGATTACTTGCCGGACACCTTGGCGTGGAAGCCGACGATACCCATCAACATGGCATGATGCTGGACGCCGATAGCCACCAGATTATTGACCTCGGCGATGATTTCTACACCGTCGGGCGTCCCCATCCGATGATCGACCCAACCTTACGCAACCAGTTAATTGCCGATCTCGGCGCTAAACCGCAAGTGCGCGTGTTGCTGCTTGATGTCGTGATTGGCTTCGGTGCGACCGCCGATCCTGCCGCCTCGCTGGTGAGCGCCTGGCAAAAAGCCTGTGCCGCGCGTTTAGATAATCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGC  >  NC_000913/547651‑548133
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  
aaTAGCGCCTGTCAGCAGCGGATTTATTTGCGGTTTGTATACCGGCGGTACGCTGGCTGCCGAAGCGGCGGGATTACTTGCCGGACACCTTGGCGTGGAAGCCGACGATACCCATCAACATGGCATGATGCTGGACGCCGATAGCCACCAGATTATTGACCTCGGCGATGATTTCTACACCGTCGGGCGTCCCCATCCGATGATCGACCCAACCTTACGCAACCAGTTAATTGCCGATCTTGGCGCTAAAc                                                                                                                                                                                                                                          <  1:222694/251‑1 (MQ=255)
aaTAGCGCCTGTCAGCAGCGGATTTATTTGCGGTTTGTATACCGGCGGTACGCTGGCTGCCGAAGCGGCGGGATTACTTGCCGGACACCTTGGCGTGGAAGCCGACGATACCCATCAACATGGCATGATGCTGGACGCCGATAGCCACCAGATTATTGACCTCGGCGATGATTTCTACACCGTCGGGCGTCCCCATCCGATGATCGACCCAACCTTACGCAACCAGTTAATTGCCGATCTTGGCGCTAAAc                                                                                                                                                                                                                                          <  1:250732/251‑1 (MQ=255)
       ccTGTCAGCAGCGGATTTATTTGCGGTTTGTATACCGGCGGTACGCTGGCTGCCGAAGCGGCGGGATTACTTGCCGGACACCTTGGCGTGGAAGCCGACGATACCCATCAACATGGCATGATGCTGGACGCCGATAGCCACCAGATTATTGACCTCGGCGATGATTTCTACACCGTCGGGCGTCCCCATCCGATGATCGACCCAACCTTACGCAACCAGTTAATTGCCGATCTTGGCGCTAAACCGCAAGt                                                                                                                                                                                                                                   <  1:9959/251‑1 (MQ=255)
          gTCAGCAGCGGATTTATTTGCGGTTTGTATACCGGCGGTACGCTGGCTGCCGAAGCGGCGGGGTTACTTGCCGGGCACCTTGGCGGGGGAGCCGACGATACCCATCAAAATGGCATGATGCTGGACGCCGATAGCCACCAAATTATTGACCCCGGCGATGATTTTTACACCGTCGGGCGTCCCCATCCGATGATCCACCCAACCTTACGCAACCCGTTAATTTCCGATCTTTGCGCTAAACCCCaaaagcg                                                                                                                                                                                                                                >  1:272047/1‑246 (MQ=255)
                   ggATTTATTTGCGGTTTGTATACCGGCGGTACGCTGGCTGCCGAAGCGGCGGGATTACTTGCCGGACACCTTGGCGTGGAAGCCGACGATACCCATCAACATGGCATGATGCTGGACGCCGATAGCCACCAGATTATTGACCTCGGCGATGATTTCTACACCGGCGGGCGTCCCCATCCGATGATCGACCCAACCTTACGCAACCAGTTAATTGCCGATGTTGGCGCTAAACCGCAAGTGCGCGTGTtgct                                                                                                                                                                                                                       <  2:191528/251‑1 (MQ=255)
                             gCGGTTTGTATACCGGCGGTACGCTGGCTGCCGAAGCGGCGGGATTACTTGCCGGACACCTTGGCGTGGAAGCCGACGATACCCATCAACATGGCATGATGCTGGACGCCGATAGCCACCAGATTATTGACCTCGGCGATGATTTCTACACCGTCGGGCGTCCCCATCCGATGATCGACCCAACCTTACGCAACCAGTTAATTGCCGATCTTGGCGCTAAACCGCAAGTGCGCGTGTTGCTGCTTGATGt                                                                                                                                                                                                              <  1:214763/250‑1 (MQ=255)
                               ggTTTGTATACCGGCGGTACGCTGGCTGCCGAAGCGGCGGGATTACTTGCCGGGCACCCTGGGGGGGGAACCGACGATAACCCACAACATGGGATGATGCTGGGCGCCGATAGCCACCCGAATATTTACCCCGGCGATGATTTCTACACCGCCGGGCGTCCCCATCCGGTGATTGACCCCACCTTACGCCACCAGTTAAATTCCGCTCTTTGGGCTCAACCGCAAGGGCGCGCGTTTCTGCTTGCTGTCgc                                                                                                                                                                                                           >  1:180503/1‑250 (MQ=255)
                                gTTTTTATACCGGCGGTACGCTGGCTGCCGAAGCGGCGGGATTACTTGCCGGACACCTTGGCGTGGAAGCCGACGATACCCATCAACATGGCATGATGCTGGACGCCGATAGCCACCAGATTATTGACCGCGGCGATGATTTCTACACCGTCGGGCGTCCCCATCCGATGATCGACCCAACCTTACGCAACCAGTTAATTGCCGATCTTGGCGCTAAACCGCAAGTGCGCGTGTTGCTGCTTGATGTCGTg                                                                                                                                                                                                          <  2:117186/251‑1 (MQ=255)
                                  ttGTATACCGGTGGTACGCTGGCTGCCGAAGCGGCGGGATTACTTGCCGGACACCTTGGCGTGGAAGCCGACGATACCCATCAACATGGCATGATGCTGGACGCCGATAGCCACCAGATTATTGACCTCGGCGATGATTTCTACACCGTCGGGCGTCCCCATCCGATGATCGACCCAACCTTACGCAACCAGTTAATTGCCGATCTTGGCGCTAAACCGCAAGTGCGCGTGTTGCTGCTTGATGTCGtta                                                                                                                                                                                                         <  2:222480/250‑3 (MQ=255)
                                         ccGGCGGTACGCTGGCTGCCGAAGCGGCGGGATTACTTGCCGGACACCTTGGCGTGGAAGCCGACGATACCCATCAACATGGCATGATGCTGGACGCCGATAGCCACCAGATTATTGACCTCGGCGATGATTTCTACACCGTCGGGCGTCCCCATCCGATGATCGACCCAACCTTACGCAACCAGTTAATTGCCGATCTTGGCGCTAAACCGCAAGTGCGCGTGTTGCTGCTTGATGTCGTGATTGGCTTc                                                                                                                                                                                                 <  1:217031/251‑1 (MQ=255)
                                                                    cggGATTACTTGCCGGACACCTTGGCGTGGAAGCCGCCGATACCTATCAACATGGCATGAGGCTGGCCGCCGATAGCCCCCAGATTATTGACCTCGCCGATGATTTCTACACCGTCGGGCGTCCCCACCCGATGATCGCCCCAACTTTCCGCAACCAGTTAATTGCCGATCTTGGCGCTAAACCGAAAGTGCGCGTGTTGCGGCTTGATGTCGTGATGGGCTTCGGTGCGACCGCCGATCCTGCCGCCTCg                                                                                                                                                                      <  1:224094/251‑1 (MQ=255)
                                                                      ggATTACTTGCCGGACACCTTGGCGTGGAAGCCGACGATACCCATCAACATGGCATGATGCTGGACGCCGATAGCCACCAGCTTATTGACCTCGGCGATGATTTCTACACCGTCGGGCGTCCCCATCCGCTGATCGACCCAACCTTACGCAACCAGTTAATTGCCGATCTTGGCGCTAAACCGCAAGTGCGCGTGTTGCTGCTTGATGTCGTGATTGGCTTCGGTGCGACCGCCGATCCTGCCGCCTCGcc                                                                                                                                                                    >  2:58956/1‑250 (MQ=255)
                                                                              ggccggcctccTTGGCGTGGAACCCGACGATACCCTCCAACATGCCATGAGGCGGGACGCGGATAGCCACCAGATTATTGACCTCGGCGATGTTTTCTACACCGTCGGGCTTCCCCATCCGATGTCCGACCCAACCTTACGCAACCAGTTAATGGCCGATCTTGGCGCTAAACCGCAAGTGCGCTTGTTGCTGCTTGATGGCGTGATTGGCTTCGGTGCGACCGCCGATCCTGCCGCCTCGCTGGTGAgcg                                                                                                                                                            <  1:156664/242‑1 (MQ=255)
                                                                                   gACACCTTGGCGTGGAAGCCGACGATACCCATCAACATGGCATGATGCTGGACGCCGATAGCCACCAGATAATTGACCTCGGCGATGCTTTCTACACCGTCGGGCGTCCCCATCCGATGATCGACCCAACCTTACGCAACCAGTTAATTGCCGATCTTGGCGCTAAACCGCAAGTGCGCGTGTTGCTGCTTGATGTCGTGATTGGCTTCGGAGCGACAGCCGATCCTGCCGCCTCGCTGGTGAGCGCCTgg                                                                                                                                                       >  1:211071/1‑251 (MQ=255)
                                                                                                       gacgaTACCCATCAACATGGCATGATGATGGACGCCGATAGCCACCAGATTATTGACCTCGGCGATGATTTCTACACCGTCGGGCGTCCCCATCCGATGATCGACCCAACCTTACGCAACCAGTTAATTGCCGATCTTGGCGCTAAACCGCAAGTGCGCGTGTTGCTGCTTGATGTCGTGATTGGCTTCGGTGCGACCGCCGATCCTGCCGCCTCGCTGGTGAGCGCCTGGCAAAAAGCCTGTGCcgcgcg                                                                                                                                   >  2:215921/1‑251 (MQ=255)
                                                                                                                  aTCAACATGGCATGATGCTGGACGCCGATAGCCACCAGATTATTGACCTCGGCGATGATTTCTACACCGTCGGGCGTCCCCATCCGATGATCGACCCAACCTTACGCAACCAGTTAATTGCCGATCTTGGCGCTAAACCGCAAGTGCGCGTGTTGCTGCTTGATGTCGTGATTGGCTTCGGTGCGACCGCCGATCCTGCCGCCTCGCTGGTGAGCGCCTGGCAAAAAGCCTGTGCCGCGCGTTTAGATAAt                                                                                                                         >  1:299897/1‑251 (MQ=255)
                                                                                                                    cAACATGGCATGATGCTGGACGCCGATAGCCACCAGATTATTGACCTCGGCGATGATTTCTACACCGTCGGGCGTCCCCATCCGATGATCGACCCAACCTTACGCAACCAGTTAATTGCCGATCTTGGCGCTAAACCGCAAGTGCGCGTGTTGCTGCTTGATGTCGTGATTGGCTTCGGTGCGACCGCCGATCCTGCCGCCTCGCTGGTGAGCGCCTGGCAAAAAGCCTGTGCCGCGCGTTTAGATAATc                                                                                                                        >  2:208365/1‑250 (MQ=255)
                                                                                                                     aaCATGGCATGATGCTGGACGCCGATAGCCACCAGATTATTGACCTCGGCGATGATTTCTACACCGTCGGGCGTCCCCATCCGATGATCGACCCAACCTTACGCAACCAGTTAATTGCCGATCTTGGCGCTAAACCGCAAGTGCGCGTGTTGCTGCTTGATGTCGTGATTGGCTTCGGTGCGACCGCCGATCCTGCCGCCTCGCTGGTGAGCGCCTGGCAAAAAGCCTGTGCCGCGCGTTTAGATAATCa                                                                                                                       >  2:127375/1‑250 (MQ=255)
                                                                                                                                   cTGGACGCCGATAGCCACCAGATTATTGACCTCGGCGATGATTTCTACACCGTCGGGCGTCCCCATCCGATGATCGACCCAACCTTACGCAACCAGTTAATTGCCGATCTTGGCGCTAAACCGCAAGTGCGCGTGTTGCTGCTTGATGTCGTGATTGGCTTCGGTGCGACCGCCGATCCTGCCGCCTCGCTGGTGAGCGCCTGGCAAAAAGCCTGTGCCGCGCGTTTAGATAATCAACCACTGTATGCCAt                                                                                                        >  1:93809/1‑251 (MQ=255)
                                                                                                                                                 ccaccaGATTATTGACCTCGGCGATGATTTCTACACCGTCGGGCGTCCCCATCCGATGATCGACCCAACCTTACGCAACCAGTTAATTGCAGATCTTGGCGCTAAACCGCAAGTGCGCGTGTTGCTGCTTGATGTCGTGATTGGCTTCGGTGCGACCGCCGATCCTGCCGCCTCGCTGGTGAGCGCCTGGCAAAAAGCCTGTGCCGCGCGTTTAGATAATCAACCACTGTATGCCATTGCCACGGTGAc                                                                                            >  2:312779/1‑249 (MQ=255)
                                                                                                                                                 ccaccaGATTATTGACCTCGGCGATGATTTCTACACCGTAGGGCGTCCCCATACGCTAATCGACCCAACCTTACGCAACAAGTTAATTGCCGATCTTGGCGCTAAACCGCAAGTGCGCGTGTTGATGCTTGATGTCGTGATTGGCTTCGGTGCGACCGCCGATCCTGCCGCATCGCTGGTGCGCGCCTGGCAAAAAGCCTGTGCCGCGCGGTTAGCTAATCAACCACTGTATGCCATTGCCACGGTGAc                                                                                            >  2:161426/1‑249 (MQ=255)
                                                                                                                                                     caGATTAGTGCCCTGGGCGATGATTTCTACACCGCCGGGCGCCCCCATCCGATGATCGCCCCAACCTTACGCAACCAGTTAATTGCCGATCTGGGCGCAAACCCGCAAGTGCGCGTGTTGCTGCTTGATGTCGTGATTGGCTTCGGTGCGACCGCCGATCCTGCCGCCTCGCTGGTGAGCGCCTGGCAAAAAGCCTGTGCCGCGCGTTTAGATAATCAACCACTGTATGCCATTGCCACGGTGACAGGCAc                                                                                      <  1:284079/251‑1 (MQ=255)
                                                                                                                                                            ttgcccTCGGCGATGATTTCTCCCCCTTCGGGCGTCCCCATCCGATGATCGACCCAACCTTACGCACCCAGTTAATTGCCGATCTTGGCGCTAAACCGCAAGTGCGCGTTTTGCTGCTTGATGTCGTGATTGGCTTCGGTGCGACCGCCGATCCTGCCGCCTCGCTGGTGAGCGCCTGGCAAAAAGCCTGTGCCGCGCGTTTAGATAATCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGAACg                                                                                <  1:168537/246‑1 (MQ=255)
                                                                                                                                                              gACCTCGGCGATGATTTCTACACCGTCGGGCGTCCCCATCCGATGATCGACCCAACCTTACGCAACCAGTTAATTGCCGATCTTGGCGCTAAACCGCAAGTGCGCGTGTTGCTGCTTGATGTCGTGATTGGCTTCGGTGCGACCGCCGATCCTGCCGCCTCGCTGGTGAGCGCCTGGCAAAAAGCCTGTGCCGCGCGTTTAGATAATCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGCACGTg                                                                              >  2:175247/1‑250 (MQ=255)
                                                                                                                                                                             ttctcccccGTCGGGCGCCCCCACCCGATGATCGACCAAACCTTACGCAACCAGTTAATTGCCGTTCTTGGCGTTAAACCGCAAGGGCGCGTGTTGCTGCTTGATGGCGTGTTTGGCTCCGGTGCGACCGCCGATCCTGCCGCCTCGCTGGTGAGCGCCTGGCAAAAAGCCTGTGCCGCGCGTTTAGATAATCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCt                                                               <  2:166607/243‑1 (MQ=255)
                                                                                                                                                                             ttCTACACCGTCGGGCGTCCCCATCCGATGATCGACCCAACCTTACGCAACCAGTTAATTGCCGATCTTGGCGCTAAACCGCAAGTGCGCGTGTTGCTGCTTGATGTCGTGATTGGCTTCGGTGCGACCGCCGATCCTGCCGCCTCGCTGGTGAGCGCCTGGCAAAAAGCCTGTGCCGCGCGTTTAGATAATCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTc                                                              >  2:73971/1‑251 (MQ=255)
                                                                                                                                                                              tCTACACCGTCGGGCGTCCCCATCCGATGATCGACCCAACCTTACGCAACCAGTTAATTGCCGATCTTGGCGCTAAACCGCAAGTGCGCGTGTTGCTGCTTGATGTCGTGATTGGCTTCGGTGCGACCGCCGATCCTGCCGCCTCGCTGGTGAGCGCCTGGCAAAAAGCCTGTGCCGCGCGTTTAGATAATCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCg                                                             >  2:82846/1‑251 (MQ=255)
                                                                                                                                                                               cTACACCGTCGGGCGTCCCCATCCGATGATCGACCCAACCTTACGCCACCAGTTAATTGCCGATCTTGGCGCTAAACCGCAAGTGCGCGTGTTGCTGCTTGATGTCGTGATTGGCTTCGGTGCGACCGCCGATCCTGCCGCCTCGCTGGTGAGCGCCTGGCAAAAAGCCTGTGCCGCGCGTTTAGATAATCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCgc                                                            >  1:34763/1‑251 (MQ=255)
                                                                                                                                                                               cTACACCGTCGGGCGTCCCCATCCGATGATCGACCCAACCTTACGCAACCAGTTAATTGCCGATCTTGGCGCTAAACCGCAAGTGCGCGTGTTGCTGCTTGATGTCGTGATTGGCTTCGGTGCGACCGCCGATCCTGCCGCCTCGCTGGTGAGCGCCTGGCAAAAAGCCTGTGCCGCGCGTTTAGATAATCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCgc                                                            >  1:286523/1‑251 (MQ=255)
                                                                                                                                                                                  aacCGTCGGGCGTCCCCATCCGATGATCGACCCAACCTTACGCAACCAGTTCATTGCCGATCTTGGCGCTACACCGCACGTGCGCGTGTTGCTGCTTGCTGTCGTGATTGGCTTCGGTGCGACCGCCGATCCTGCCGCCTCGCTGGTGAGCGCCTGGCAAAAAGCCTGTGCCGCGCGTTTAGATAATCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCgcagc                                                         >  2:52455/2‑251 (MQ=255)
                                                                                                                                                                                      gTCGGGCGTCCCCATCCGATGATCGACCCAACCTTACGCAACCAGTTAATTGCCGATCTTGGCGCTAAACCGCAAGTGCGCGTGTTGCTGCTTGATGTCGTGATTGGCTTCGGTGCGACCGCCGATCCTGCCGCCTCGCTGGTGAGCGCCTGGCAAAAAGCCTGTGCCGCGCGTTTAGATAATCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAAt                                                     >  2:101209/1‑251 (MQ=255)
                                                                                                                                                                                         gggCGTCCCCATCCGATGATCGACCCAACCTTACGCAACCAGTTAATTGCCGATCTTGGCGCTAAACCGCAAGTGCGCGTGTTGCTGCTTGATGTCGTGATTGGCTTCGGTGCGACCGCCGATCCTGCCGCCTCGCTGGTGAGCGCCTGGCAAAAAGCCTGTGCCGCGCGTTTAGATAATCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGc                                                  >  1:47863/1‑251 (MQ=255)
                                                                                                                                                                                              tCCCCATCCGATGATCGCCCCAACCTTACGCAACCAGTTAATTGCCGATCTTGGCGCTAAACCGCAAGTGCGCGTGTTGCTGCTTGATGTCGTAATTGGCTTCGGTGCGACCGCCGATCCTGCCGCCTCGCTGGTGAGCGCCTGGCAAAAAGCCTGTGCCGCGCGTTTAGATAATCAACCACTGTATGCCATTGCCCCGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACg                                              >  1:241114/1‑250 (MQ=255)
                                                                                                                                                                                                cccATCCGATGATCGACCCAACCTTACGCAACCAGTTAATTGCCGATCTTGGCGCTAAACCGCAAGTGCGCGTGTTGCTGCTTGATGTCGTGATTGGCTTCGGTGCGACCGCCGATCCTGCCGCCTCGCTGGTGAGCGCCTGGCAAAAAGCCTGTGCCGCGCGTTTAGATAATCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTg                                           >  1:159628/1‑251 (MQ=255)
                                                                                                                                                                                                     ccGATGATCGACCCAACCTTACGCAACCAGTTAATTGCCGATCTTGGCGCTAAACCGCAAGTGCGCGTGTTGCTGCTTGATGTCGTGATTGGCTTCGGTGCGACCGCCGATCCTGCCGCCTCGCTGGTGAGCGCCTGGCAAAAAGCCTGTGCCGCGCGGTTAGATAATCAACCACTGTATGCCATTGCCACGGTGCCAGGCCCGGAACGTGACCAGCAATGCCGCTCGCAGCCACACGCCACGCTGGa                                         >  2:175438/1‑248 (MQ=255)
                                                                                                                                                                                                               ccccAACCTTACGCACCCGTTTAATTGCCGATCTTGGCGCAAAACCGCAAGTGCGCGTGTTGCTGCTTGATGTCGTGATTGGCTGCGGTGCGACCGCCGATCCTGCCGCCTCGCTGGTGAGCGCCTGGCAAAAAGCCTGTGCCGCGCGTTTAGATAATCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGAtt                            <  1:243121/250‑1 (MQ=255)
                                                                                                                                                                                                                     ccttccgaaacaaGTTAATTGCCTATCTTGGCGCTAACCCGCAAGTGCGCGTGTTGCGGCTTGATGTCTGGATTGGCTTCGGTGCGACCGCCGATCCTGCCGCCTCGCTGGTGAGCGCCTGGCAAAAAGCCTGTGCCGCGCGTTTAGATAATCAACCCCTGTATGCCATTGCCTCGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTc                      <  1:53354/239‑1 (MQ=255)
                                                                                                                                                                                                                       ttACGCAACCAGTTAATTGCCGATCTTGGCGCTAAACCGCAAGTGCGCGTGTTGCTGCTTGATGTCGTGATTGGCTTCGGTGCGACCGCCGATCCTGCCGCCTCGCTGGTGAGCGCCTGGCAAAAAGCCTGTGCCGCGCGTTTAGATAATCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGt                    >  1:66734/1‑251 (MQ=255)
                                                                                                                                                                                                                            cacccAGTTAATTGCCGAGCTGGGCGCAAACCCGCAGGTGCGCGTGTTTCTGCTTGAGGTCGTGCTTGGCTTCGGTGCGACCGCCGATCCTGCCGCCTCGCTGGTGAGCGCCTGGCAAAAAGCCTGTGCCGCGCGTTTAGATAATCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGtt               <  1:222405/248‑1 (MQ=255)
                                                                                                                                                                                                                                          ccGATCTCGGCGATAATCCGCAAGTGCGCGGGTTGCTGCTTGATGTCGTGATTGGCTTCGATGCGACCGCAGATCCTGCCGCCTCGCTGGTGAGCGCCTGACAAAAAGCCTGTGCCGCGCGTTTCGATAATCAACCACTGTATGCCATTGCCACGGTGACAGGCCCTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGGGGGGATTGCGGTCGTGCGCTCGATACCGGAAGc  >  2:127547/1‑250 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  
AATAGCGCCTGTCAGCAGCGGATTTATTTGCGGTTTGTATACCGGCGGTACGCTGGCTGCCGAAGCGGCGGGATTACTTGCCGGACACCTTGGCGTGGAAGCCGACGATACCCATCAACATGGCATGATGCTGGACGCCGATAGCCACCAGATTATTGACCTCGGCGATGATTTCTACACCGTCGGGCGTCCCCATCCGATGATCGACCCAACCTTACGCAACCAGTTAATTGCCGATCTCGGCGCTAAACCGCAAGTGCGCGTGTTGCTGCTTGATGTCGTGATTGGCTTCGGTGCGACCGCCGATCCTGCCGCCTCGCTGGTGAGCGCCTGGCAAAAAGCCTGTGCCGCGCGTTTAGATAATCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGC  >  NC_000913/547651‑548133

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