Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 299,377 A→T 7.2% Y520N (TAC→AAC)  paoC ← PaoABC aldehyde oxidoreductase, Moco‑containing subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913299,3770AT 7.2% 281.6 / 7.0 82Y520N (TAC→AAC) paoCPaoABC aldehyde oxidoreductase, Moco‑containing subunit
Reads supporting (aligned to +/- strand):  ref base A (34/42);  new base T (2/4);  total (36/46)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGC  >  NC_000913/299235‑299523
                                                                                                                                              |                                                                                                                                                  
gTCAGTCTGCCGCCTGCGGTGGCTTCAGGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGcc                                                                                                                                              <  4:763190/149‑1 (MQ=255)
   tgTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGACGCGTCAGGTACCGTTGGTAATCTTGCCGTCGGCAAACGGCGACAGCGCAGGATCGAACCCGACTGCCGAGGCAATCGATTCGCGAAGCTTCATACAGGCGGCGTTTACGCCgcc                                                                                                                                           <  1:246603/148‑3 (MQ=255)
    gTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTTTACGCCgcct                                                                                                                                          <  1:317381/149‑4 (MQ=255)
        gccgccTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCggaggtgg                                                                                                                                      <  3:217680/149‑1 (MQ=255)
          cgccTGCGGTGGCTTGATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCTACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCggaggtggag                                                                                                                                    <  2:756743/149‑1 (MQ=255)
          cgccTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCCTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGCTCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCgttggtggag                                                                                                                                    <  2:697738/149‑1 (MQ=255)
                    ggCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCgcgc                                                                                                                          <  3:616038/149‑1 (MQ=255)
                     gCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGcc                                                                                                                         >  4:442512/1‑149 (MQ=255)
                       ttCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTGGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAGTCATTTCGTGAGGCTTCATACAGGCGGCGTAAACGGCGGAGGTGCAAGTACTCGCGccgt                                                                                                                       >  2:840813/1‑147 (MQ=255)
                       ttCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTTTACGCCgcct                                                                                                                                          <  4:699485/130‑4 (MQ=255)
                       ttCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTTTACGCCgcct                                                                                                                                          >  3:699485/1‑127 (MQ=255)
                       ttCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGcccc                                                                                                                       >  1:221627/1‑149 (MQ=255)
                          aTGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCAtt                                                                                                                    <  4:487180/149‑1 (MQ=255)
                          aTGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCAtt                                                                                                                    <  4:783731/149‑1 (MQ=255)
                           tGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTg                                                                                                                   <  4:425552/149‑1 (MQ=255)
                                cGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTCTTCGCGCCCCATTGTccac                                                                                                              <  2:178142/149‑1 (MQ=255)
                                  tGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTccacca                                                                                                            <  1:438792/149‑1 (MQ=255)
                                   ggCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCTATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAg                                                                                                           >  2:233205/1‑149 (MQ=255)
                                     cgcTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGaa                                                                                                         >  4:793717/1‑149 (MQ=255)
                                          cGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGc                                                                                                    <  2:450100/149‑1 (MQ=255)
                                                 cGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCggaggtggag                                                                                                                                    <  2:581437/110‑1 (MQ=255)
                                                 cGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCggaggtggag                                                                                                                                    >  1:581437/1‑110 (MQ=255)
                                                  gTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCg                                                                                            >  4:622895/1‑149 (MQ=255)
                                                  gTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCg                                                                                            >  1:740364/1‑149 (MQ=255)
                                                   ttGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCgg                                                                                           <  1:616632/149‑1 (MQ=255)
                                                      gTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGaa                                                                                        >  3:133704/1‑149 (MQ=255)
                                                        aaTCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAAc                                                                                      <  2:113406/149‑1 (MQ=255)
                                                        aaTCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAAc                                                                                      <  3:442512/149‑1 (MQ=255)
                                                          tCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTg                                                                                    >  2:143944/1‑149 (MQ=255)
                                                          tCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTg                                                                                    <  3:439364/149‑1 (MQ=255)
                                                              gCCGTCGGCAAACTGCGACTGCTTAGGATCAAAACCGACTGCCGAGGAAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAAt                                                                                <  2:131178/149‑1 (MQ=255)
                                                               ccGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATc                                                                               <  2:508500/149‑1 (MQ=255)
                                                                    ggCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGccca                                                                          <  2:182274/149‑3 (MQ=255)
                                                                    ggCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGa                                                                          >  2:147603/1‑149 (MQ=255)
                                                                     gCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAAc                                                                                      <  1:210624/136‑1 (MQ=255)
                                                                     gCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAAc                                                                                      >  2:210624/1‑136 (MQ=255)
                                                                     gCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAg                                                                         <  1:101780/149‑1 (MQ=255)
                                                                       aaaCTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGt                                                                       >  1:47705/1‑149 (MQ=255)
                                                                           tGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAAc                                                                   <  4:580443/149‑1 (MQ=255)
                                                                            gCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAAcc                                                                  <  3:622895/149‑1 (MQ=255)
                                                                            gCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAAcc                                                                  <  4:35470/149‑1 (MQ=255)
                                                                                cTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCaa                                                              <  1:444737/149‑1 (MQ=255)
                                                                                 tGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAAc                                                             >  2:151516/1‑149 (MQ=255)
                                                                                    tCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTg                                                          >  3:250249/1‑149 (MQ=255)
                                                                                    tCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTg                                                          >  4:474841/1‑149 (MQ=255)
                                                                                        gATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTcc                                                      >  1:19608/1‑149 (MQ=255)
                                                                                         aTCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAAGCGCAACCTGCTCCa                                                     >  1:643973/1‑149 (MQ=255)
                                                                                         aTCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCa                                                     >  4:223654/1‑149 (MQ=255)
                                                                                          tCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGATGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGCGGTGAACCGCAACCTGCTCCAg                                                    >  3:882154/1‑149 (MQ=255)
                                                                                            aaaCCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCgg                                                                                           <  4:354464/108‑1 (MQ=255)
                                                                                            aaaCCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCgg                                                                                           >  3:354464/1‑108 (MQ=255)
                                                                                              aCCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGt                                                >  1:255061/1‑149 (MQ=255)
                                                                                              aCCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGt                                                >  3:88723/1‑149 (MQ=255)
                                                                                                ccGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTAc                                              <  2:189241/149‑1 (MQ=255)
                                                                                                     tGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCaa                                         <  2:182410/149‑1 (MQ=255)
                                                                                                         gAGGCAATCATTTCGAGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCAt                                     >  4:461472/1‑149 (MQ=255)
                                                                                                            gCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTc                                  >  1:596785/1‑149 (MQ=255)
                                                                                                              aaTCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCg                                >  2:429944/1‑149 (MQ=255)
                                                                                                               aTCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGc                               <  1:143944/149‑1 (MQ=255)
                                                                                                               aTCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGc                               <  3:879115/149‑1 (MQ=255)
                                                                                                                  aTTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGc                            <  1:233205/149‑1 (MQ=255)
                                                                                                                   tttCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGcc                           >  1:628403/1‑149 (MQ=255)
                                                                                                                   tttCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGcc                           >  4:639198/1‑149 (MQ=255)
                                                                                                                   tttCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGcc                           >  4:483765/1‑149 (MQ=255)
                                                                                                                   tttCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGcc                           >  2:258379/1‑149 (MQ=255)
                                                                                                                      cgcgAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTc                        <  3:3454/149‑1 (MQ=255)
                                                                                                                      cgcgAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTc                        <  2:19608/149‑1 (MQ=255)
                                                                                                                         gAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTgg                     <  2:547969/149‑1 (MQ=255)
                                                                                                                          aaGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTggg                    <  1:821108/149‑1 (MQ=255)
                                                                                                                               tCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGCCAGa               <  4:88723/149‑1 (MQ=255)
                                                                                                                                  tACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGCCAGAATg            <  3:786360/149‑1 (MQ=255)
                                                                                                                                     ctccggcgTTTACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGCCAGAATGgtg         >  3:276152/4‑149 (MQ=255)
                                                                                                                                     aGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGCCCGAATGgtg         >  2:845202/1‑149 (MQ=255)
                                                                                                                                     aGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGCCAGAATGgtg         <  1:808653/149‑1 (MQ=255)
                                                                                                                                      tccggcgTTTACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGaaa                                                                                       <  1:261854/68‑1 (MQ=255)
                                                                                                                                       gcggcgTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGCCAGAATGGTGTa       <  3:155917/149‑1 (MQ=255)
                                                                                                                                       gcggcgTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGCCAGAATGGTGTa       <  3:223654/149‑1 (MQ=255)
                                                                                                                                          gcgTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTc                                  >  2:604304/1‑119 (MQ=255)
                                                                                                                                          gcgTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTc                                  <  1:604304/119‑1 (MQ=255)
                                                                                                                                            gTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGCCAGAATgg           <  4:383917/140‑1 (MQ=255)
                                                                                                                                            gTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGCCAGAATgg           >  3:383917/1‑140 (MQ=255)
                                                                                                                                            gTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGc  >  3:825005/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
GTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGC  >  NC_000913/299235‑299523

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