Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,681,614 G→A 5.5% V213V (GTG→GTA folM → dihydromonapterin reductase, NADPH‑dependent; dihydrofolate reductase isozyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,681,6140GA 5.5% 341.5 / 4.7 104V213V (GTG→GTAfolMdihydromonapterin reductase, NADPH‑dependent; dihydrofolate reductase isozyme
Reads supporting (aligned to +/- strand):  ref base G (45/53);  new base A (2/4);  total (47/57)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 104 total aligned reads displayed.

CTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAAC  >  NC_000913/1681467‑1681759
                                                                                                                                                   |                                                                                                                                                 
cTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGa                                                                                                                                                  <  1:141548/149‑1 (MQ=255)
cTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGa                                                                                                                                                  <  1:38985/149‑1 (MQ=255)
  cGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATc                                                                                                                                                <  3:558765/149‑1 (MQ=255)
  cGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATc                                                                                                                                                <  2:379102/149‑1 (MQ=255)
   gTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCg                                                                                                                                               >  1:393036/1‑149 (MQ=255)
    tttGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGa                                                                                                                                              <  4:743252/149‑1 (MQ=255)
       gCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCt                                                                                                                                           <  1:480640/149‑1 (MQ=255)
        cccGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTg                                                                                                                                          <  2:252547/149‑1 (MQ=255)
        cccGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTg                                                                                                                                          <  1:393707/149‑1 (MQ=255)
         ccGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGACGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTgg                                                                                                                                         >  4:251974/1‑149 (MQ=255)
          cGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGt                                                                                                                                        >  2:635420/1‑149 (MQ=255)
              aGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAt                                                                                                                                    <  2:604095/149‑1 (MQ=255)
              aGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAt                                                                                                                                    <  4:491967/149‑1 (MQ=255)
                 tGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttac                                                                                                                                 <  1:823465/149‑1 (MQ=255)
                 tGGCACCGGAAGTGAAAGTGAATTCTAATGCGCCATCGCTGATCGTGTTTAATGAACATGACGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttac                                                                                                                                 <  3:644819/149‑1 (MQ=255)
                  ggCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTAATCGACCAGGTCGAtcact                                                                                                                                <  3:785060/149‑1 (MQ=255)
                   gCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTAATCGACCAGGTCGAtcac                                                                                                                                 <  2:833613/147‑1 (MQ=255)
                   gCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTAATCGACCAGGTCGAtcac                                                                                                                                 >  1:833613/1‑144 (MQ=255)
                    cACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttactta                                                                                                                              <  3:257077/149‑1 (MQ=255)
                    cACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTAATCGACCAGGTCGAtcacttc                                                                                                                              <  4:330479/149‑2 (MQ=255)
                     aCCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttacttac                                                                                                                             >  4:796710/1‑149 (MQ=255)
                     aCCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttacttac                                                                                                                             <  4:609289/149‑1 (MQ=255)
                        ggAAGTGAAAGAGAATTCTATTGCGCCATCGCTGATCCTTTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATGAATCACTGATTAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttacttactta                                                                                                                          <  2:261072/149‑1 (MQ=255)
                          aagtgaaagtgaATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTAcc                                                                                                                        <  4:387074/149‑1 (MQ=255)
                           agtgaaagtgaaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCa                                                                                                                       <  3:473657/149‑1 (MQ=255)
                            gtgaaagtgaaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAg                                                                                                                      >  1:833242/1‑149 (MQ=255)
                                   tgaaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCttt                                                                                                               >  2:789677/1‑149 (MQ=255)
                                   tgaaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCttt                                                                                                               <  2:518961/149‑1 (MQ=255)
                                   tgaaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCttt                                                                                                               <  1:618321/149‑1 (MQ=255)
                                        tCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCAc                                                                                                          <  3:726472/149‑1 (MQ=255)
                                                 ccATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAg                                                                                                 <  1:408401/149‑1 (MQ=255)
                                                  cATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGt                                                                                                >  4:436385/1‑149 (MQ=255)
                                                   aTCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATTACTGATGAAAATCGCGCCTGGCGAGAAAGAAGAGATCGACCTAGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGtt                                                                                               >  4:423305/1‑149 (MQ=255)
                                                      gCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTcc                                                                                            >  3:598597/1‑149 (MQ=255)
                                                      gCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTcc                                                                                            >  3:388346/1‑149 (MQ=255)
                                                         gATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCttt                                                                                                               >  4:353589/1‑127 (MQ=255)
                                                         gATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCttt                                                                                                               <  3:353589/127‑1 (MQ=255)
                                                          aTCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATc                                                                                                                                                >  3:199485/1‑93 (MQ=255)
                                                              tGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGAt                                                                                    <  1:700182/149‑1 (MQ=255)
                                                              tGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGAt                                                                                    <  1:635420/149‑1 (MQ=255)
                                                               gTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGGGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATg                                                                                   >  1:116836/1‑149 (MQ=255)
                                                                  tAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTAGTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTCGATGGCg                                                                                >  3:760093/1‑149 (MQ=255)
                                                                   aaTGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGGTTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCgg                                                                               >  4:76316/1‑149 (MQ=255)
                                                                    aTGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGgt                                                                              >  4:344219/1‑149 (MQ=255)
                                                                     tGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTAATCGACCAGGTCGAtcactt                                                                                                                               <  1:836788/99‑1 (MQ=255)
                                                                     tGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCACCTGGCGAGAAAGAAGTAATCGACCAGGTCGAtcactt                                                                                                                               >  2:836788/1‑99 (MQ=39)
                                                                      gAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGt                                                                                                             >  1:430604/1‑116 (MQ=255)
                                                                      gAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGt                                                                                                             <  2:430604/116‑1 (MQ=255)
                                                                        aCATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGgtcgtc                                                                          >  3:242877/1‑149 (MQ=255)
                                                                         cATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCa                                                                         >  1:749862/1‑149 (MQ=255)
                                                                          atgatgatgCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCAt                                                                        <  2:699628/149‑1 (MQ=255)
                                                                           tgatgatgCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCa                                                                                                                       >  4:432105/1‑101 (MQ=255)
                                                                           tgatgatgCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCa                                                                                                                       <  3:432105/101‑1 (MQ=255)
                                                                            gatgatgCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCt                                                                      >  4:327455/1‑149 (MQ=255)
                                                                             atgatgCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTg                                                                     >  1:269920/1‑149 (MQ=255)
                                                                                atgCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGt                                                                  <  4:512675/149‑1 (MQ=255)
                                                                                   ccGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTaa                                                               <  1:106991/149‑1 (MQ=255)
                                                                                   ccGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTaa                                                               <  2:831133/149‑1 (MQ=255)
                                                                                    cGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCttt                                                                                                               <  3:741850/100‑1 (MQ=255)
                                                                                    cGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCttt                                                                                                               >  4:741850/1‑100 (MQ=255)
                                                                                    cGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAAt                                                              <  2:179710/149‑1 (MQ=255)
                                                                                    cGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAAt                                                              <  4:242877/149‑1 (MQ=255)
                                                                                    cGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAAt                                                              <  2:484049/149‑1 (MQ=255)
                                                                                     gAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATg                                                             >  3:498276/1‑149 (MQ=255)
                                                                                     gAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATg                                                             >  2:468181/1‑149 (MQ=255)
                                                                                      aaTATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGCCGTCATCTGCGTTAATGc                                                            >  4:695580/1‑149 (MQ=255)
                                                                                            gACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCttt                                                      >  3:301407/1‑149 (MQ=255)
                                                                                            gACAACAGGCGCTGAATAAATCACAGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACCTACCAGTTGCTTTGTCACCGGACGCCGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGTAGCttt                                                      >  2:733484/1‑149 (MQ=255)
                                                                                                acaGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATcc                                                  <  1:468181/149‑1 (MQ=255)
                                                                                                 caGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCa                                                 <  2:576059/149‑1 (MQ=255)
                                                                                                 caGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCa                                                 >  3:856237/1‑149 (MQ=255)
                                                                                                      gcTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCata                                            <  2:285769/149‑1 (MQ=255)
                                                                                                         gAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATaa                                         <  1:387013/149‑1 (MQ=255)
                                                                                                         gAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATCACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCAGCGTTAATGCAGCTTTATCCAGCATATaa                                         <  3:48981/149‑1 (MQ=255)
                                                                                                             aaaTCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACGAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTCTATCCAGCATATAATCaa                                     <  1:259659/149‑1 (MQ=255)
                                                                                                               aTCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAAt                                                              <  1:20241/122‑1 (MQ=255)
                                                                                                               aTCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAAt                                                              >  2:20241/1‑122 (MQ=255)
                                                                                                                tCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAacac                                  >  4:579244/1‑149 (MQ=255)
                                                                                                                 cACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACAcc                                 <  3:847008/149‑1 (MQ=255)
                                                                                                                   cTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCa                               >  3:273185/1‑149 (MQ=255)
                                                                                                                    tggtgaCAACCGCGCCTGTCGAGAACGTAGTGATCGACCTGGGCGCCTACTTACTTACCAGTTGCTTTGTCACCGTAGGCAGTTTCCCACTTGATGGCGGTCGTCATCTGAGTTAATGCAGCTTTATCCAGCTTATAATCAACACCGAt                              >  2:812849/4‑149 (MQ=255)
                                                                                                                    tgatgaAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttact                                                                                                                                >  1:509050/1‑51 (MQ=255)
                                                                                                                    tgatgaAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttact                                                                                                                                <  2:509050/51‑1 (MQ=255)
                                                                                                                         aaaaCCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCt                         >  2:597385/1‑149 (MQ=255)
                                                                                                                         aaaaCCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCt                         >  2:537994/1‑149 (MQ=255)
                                                                                                                         aaaaCCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTCATCCAGCATATAATCAACACCCATGCGCt                         >  3:394175/1‑149 (MQ=255)
                                                                                                                           aaCCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTa                       >  4:829649/1‑149 (MQ=255)
                                                                                                                            aCCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTAt                      <  4:161572/149‑1 (MQ=255)
                                                                                                                                 gcCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCCTGCGCTTATccccc                 <  1:350639/149‑1 (MQ=255)
                                                                                                                                 gcCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATccccc                 <  3:796710/149‑1 (MQ=255)
                                                                                                                                 gcCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATccccc                 >  1:517149/1‑149 (MQ=255)
                                                                                                                                    tGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTccc                                                                                           <  4:722584/72‑1 (MQ=255)
                                                                                                                                    tGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTccc                                                                                           >  3:722584/1‑72 (MQ=255)
                                                                                                                                        gacaaagaaGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCACCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGc          >  4:591720/4‑149 (MQ=255)
                                                                                                                                          caaagaaGTGATCGACCTGGTCGATTACTTACTTACCAGTCGCTTTGTCACCGCACGCAGTTTCCCACTTGACGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCAAACACGCaa        <  4:760093/148‑1 (MQ=255)
                                                                                                                                           aaagaaGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAAc       <  3:544552/149‑1 (MQ=255)
                                                                                                                                              gaagTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGa    >  4:229410/1‑149 (MQ=255)
                                                                                                                                                agTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAAc  <  4:388346/149‑1 (MQ=255)
                                                                                                                                                  tGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAg                                                          <  1:163040/91‑1 (MQ=255)
                                                                                                                                                  tGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAg                                                          >  2:163040/1‑91 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAAC  >  NC_000913/1681467‑1681759

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