Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,681,622 T→A 5.5% L216Q (CTG→CAG)  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,6220TA 5.5% 319.7 / 4.3 105L216Q (CTG→CAG) folMdihydromonapterin reductase, NADPH‑dependent; dihydrofolate reductase isozyme
Reads supporting (aligned to +/- strand):  ref base T (42/56);  new base A (2/4);  total (44/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

CCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCG  >  NC_000913/1681475‑1681770
                                                                                                                                                   |                                                                                                                                                    
cccGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTgg                                                                                                                                                    >  1:704238/1‑150 (MQ=255)
cccGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTgg                                                                                                                                                    >  1:692528/1‑150 (MQ=255)
 ccGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCAGGt                                                                                                                                                   <  2:1790859/150‑1 (MQ=255)
  cGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTc                                                                                                                                                  <  2:1301865/150‑1 (MQ=255)
    cAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGa                                                                                                                                                <  2:171063/150‑1 (MQ=255)
       gCTGGCACCGGAAGTGAAAGGGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAtta                                                                                                                                             <  1:800646/150‑1 (MQ=255)
        cTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttac                                                                                                                                            <  1:1072300/150‑1 (MQ=255)
         tGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATATATCACTGATTAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttact                                                                                                                                           <  2:1331864/150‑1 (MQ=255)
          ggCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttactt                                                                                                                                          >  1:1431025/1‑150 (MQ=255)
          ggCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTAATCGACCAGGTCGAtcactt                                                                                                                                          <  1:1206645/150‑1 (MQ=255)
          ggCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATACTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttactt                                                                                                                                          <  1:746537/150‑1 (MQ=255)
          ggCACCGGAAGTGAAAGGGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttactt                                                                                                                                          <  1:1598778/150‑1 (MQ=255)
           gCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttactta                                                                                                                                         <  2:827490/150‑1 (MQ=255)
           gCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttactta                                                                                                                                         <  2:834705/150‑1 (MQ=255)
            cACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttacttac                                                                                                                                        >  1:494790/1‑150 (MQ=255)
            cACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTAATCGACCAGGTCGAtcacttct                                                                                                                                        <  2:430262/150‑3 (MQ=255)
             aCCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGGGAGAAAGAAGTGATCGACCTGGTCGAttacttact                                                                                                                                       <  2:1132029/150‑1 (MQ=255)
             aCCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttacttact                                                                                                                                       <  2:677274/150‑1 (MQ=255)
             aCCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttacttact                                                                                                                                       <  1:1080941/150‑1 (MQ=255)
              ccGGACGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATTCAGATGATGATGCCGAATATCCACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTAATCGACCAGGTCGAtcacttcttt                                                                                                                                      <  2:286955/150‑5 (MQ=255)
                ggAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttacttacttac                                                                                                                                    <  1:799221/150‑1 (MQ=255)
                  aactgaaagtgaaTTCTAATCCGCCATCGCTGATCCTGTATAATGAACCATACGCTGCCGAATATCGACAACAGGCGCTGAATAAATCACTCACGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCGGGGCGATTACTTACTTACCa                                                                                                                                  <  1:104905/147‑1 (MQ=255)
                   agtgaaagtgaaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAg                                                                                                                                 >  1:1132472/1‑150 (MQ=255)
                       aaagtgaaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGc                                                                                                                             <  2:880895/150‑1 (MQ=255)
                         agtgaaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCtt                                                                                                                           <  1:64394/150‑1 (MQ=255)
                           tgaaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTg                                                                                                                         <  2:328013/150‑1 (MQ=255)
                           tgaaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTg                                                                                                                         <  1:824722/150‑1 (MQ=255)
                           tgaaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATAATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTg                                                                                                                         <  2:1017624/150‑1 (MQ=255)
                             aaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTc                                                                                                                       <  2:1706497/150‑1 (MQ=255)
                                   aTTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGa                                                                                                                 >  2:854920/1‑150 (MQ=255)
                                    ttGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGAc                                                                                                                <  1:623158/150‑1 (MQ=255)
                                       cgcCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCa                                                                                                             >  1:628678/1‑150 (MQ=255)
                                         ccATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGt                                                                                                           >  2:1000058/1‑150 (MQ=255)
                                             cGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTcc                                                                                                       <  2:1439293/150‑1 (MQ=255)
                                                tGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCAc                                                                                                    <  2:1444552/150‑1 (MQ=255)
                                                 gATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACt                                                                                                   <  2:222809/150‑1 (MQ=255)
                                                 gATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACt                                                                                                   <  2:1261351/150‑1 (MQ=255)
                                                     cTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGAt                                                                                               <  1:1064863/150‑1 (MQ=255)
                                                     cTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGAt                                                                                               >  1:1810995/1‑150 (MQ=255)
                                                        tttAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAACTCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGc                                                                                            <  2:850801/150‑1 (MQ=255)
                                                          tAATGAACATGATGATGCCGTATATCGACAACAGGCGCTGAATAAATCTCTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGGCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCgg                                                                                          <  2:162648/150‑1 (MQ=255)
                                                          tAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGTAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCgg                                                                                          <  1:31328/150‑1 (MQ=255)
                                                          tAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCgg                                                                                          >  2:1247551/1‑150 (MQ=255)
                                                           aaTGAACATGATGATGCCGAATATCGACAAGAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGgt                                                                                         <  1:822288/150‑1 (MQ=255)
                                                           aaTGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGgt                                                                                         <  2:418961/150‑1 (MQ=255)
                                                              ggaCATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGATACCCGCGCCTGGAGAGAAAGAAGTGATCGCCCTGGTCGATTACTTACTTACCAGTTGCTTGGTCACCGGACGCAGTTTCCCACGAGATGGCGgtcga                                                                                      >  1:1285632/3‑149 (MQ=255)
                                                              gAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCCCCGGACGCAGTTTCCCACGTGATGGCGgtcgt                                                                                      >  2:972510/1‑150 (MQ=255)
                                                              gAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACGTACCAGTTGCTTTGTCACCGGCCGCAGTTTCCCACTTGATGGCGgtcgt                                                                                      >  2:45422/1‑150 (MQ=255)
                                                               aaCATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGgtcgtc                                                                                     >  2:1444873/1‑150 (MQ=255)
                                                                   tgatgatgCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCt                                                                                 >  2:1032033/1‑150 (MQ=255)
                                                                         tgCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGGAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTa                                                                           <  1:1536311/150‑1 (MQ=255)
                                                                           ccGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTAACAGATGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAAt                                                                         >  2:1126979/1‑150 (MQ=255)
                                                                            cGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATg                                                                        <  1:933065/150‑1 (MQ=255)
                                                                            cGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATg                                                                        <  2:494790/150‑1 (MQ=255)
                                                                              aaTATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCa                                                                      <  1:1154448/150‑1 (MQ=255)
                                                                                 atCGACAACAGGCGCTGAATAAATCACTGATGAAACCCGCGCCTGGCGAGAAAGAAGTGAGCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCt                                                                   <  2:133876/150‑1 (MQ=255)
                                                                                       aacaGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATcc                                                             <  2:1061121/150‑1 (MQ=255)
                                                                                         caGGCGCTTAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATAGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAg                                                           <  2:167768/150‑1 (MQ=255)
                                                                                          cggCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGc                                                          <  2:204692/149‑1 (MQ=255)
                                                                                          aGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGc                                                          <  2:704238/150‑1 (MQ=255)
                                                                                          aGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGc                                                          >  2:9585/1‑150 (MQ=255)
                                                                                           ggCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCa                                                         >  2:245579/1‑150 (MQ=255)
                                                                                            gcgcTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCat                                                        <  2:701911/150‑1 (MQ=255)
                                                                                               cTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCAtata                                                     <  1:1393790/150‑1 (MQ=255)
                                                                                                tGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATaa                                                    >  2:827134/1‑150 (MQ=255)
                                                                                                 gAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTAGCCAGCATATAAt                                                   >  2:964428/1‑150 (MQ=255)
                                                                                                  aaTAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATc                                                  <  1:137778/150‑1 (MQ=255)
                                                                                                   aTAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCa                                                 <  1:1016871/150‑1 (MQ=255)
                                                                                                     aaaTCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATAAac                                               >  1:107510/1‑150 (MQ=255)
                                                                                                           cTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCAt                                         >  2:849790/1‑150 (MQ=255)
                                                                                                              atgaAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATgcg                                      >  2:1435522/1‑150 (MQ=255)
                                                                                                                  aaaCCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGCCGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCACCACCCATGCGCTTa                                  >  2:1206625/1‑150 (MQ=255)
                                                                                                                  aaaCCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTa                                  >  2:1485966/1‑150 (MQ=255)
                                                                                                                   aaCCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTAt                                 >  1:1522467/1‑150 (MQ=255)
                                                                                                                    aCCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATc                                <  2:1431025/150‑1 (MQ=255)
                                                                                                                     ccGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTGCTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATcc                               <  2:1154080/150‑1 (MQ=255)
                                                                                                                     ccGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATcc                               >  2:931049/1‑150 (MQ=255)
                                                                                                                     ccGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATcc                               >  1:1814483/1‑150 (MQ=255)
                                                                                                                      cgcgCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATccc                              <  1:1000058/150‑1 (MQ=255)
                                                                                                                          ccTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCaa                          <  2:1765338/150‑1 (MQ=255)
                                                                                                                           cTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCTGACTCAGTTTACCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAAGACCCATGCGCTTATCCCCCAac                         <  1:1261867/150‑1 (MQ=255)
                                                                                                                           cTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAac                         <  1:608647/150‑1 (MQ=255)
                                                                                                                            tGGCGAGAAAGAAGTGATCGAACTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTACCACTTGATGGCGGTCGTCATCTGCGTTAATGAAGCTTTATCAAGCATATAATCAACACACATGCGATTATCCCCCAaca                        >  1:1022215/1‑150 (MQ=255)
                                                                                                                             ggCGAGATAGAAGTTATCGACCTGGTCGATTACTTAATTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAacac                       <  2:170166/150‑1 (MQ=255)
                                                                                                                               cGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCCCTTGATGGCGGTCGTCATCTGCGTTAATGCCGCTTTATCCAGCATATAATCAACACCCATGCGCGTAGCCCCCAACCCGc                     >  2:319248/1‑150 (MQ=255)
                                                                                                                                gagaAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCa                    <  1:9585/150‑1 (MQ=255)
                                                                                                                                  gaatgcagTGATCGACCTGGTCGATTACGTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAAc                  <  2:935149/144‑1 (MQ=255)
                                                                                                                                    aagaagTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAg                >  2:493946/1‑150 (MQ=255)
                                                                                                                                     gtaagTAATCGACCAGGTCGAGGACTTACTTACCAGTTGCTTTGTCACCGGAGGCAGTTTCCCACTTGATGGCGGGCGTCATCTGCGTTAATGAAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATGCCCCTAGACGCATCAGa               >  2:700856/3‑150 (MQ=255)
                                                                                                                                      gaagTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATTGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGaa              >  2:787384/1‑150 (MQ=255)
                                                                                                                                       aagTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATTATCACCACCCCTGCGCTTATCCCCCAACACGCAACCGAAc             >  2:893574/1‑150 (MQ=255)
                                                                                                                                        agTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAAcc            >  1:1445298/1‑150 (MQ=255)
                                                                                                                                        agTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCCGTCGTCATCTGCGTTATTGCAGCTTTATCCAGCATATAATCAACACCCATGCGCGTATCCCCCAACACGCAAAAGAAcc            >  2:639709/1‑150 (MQ=255)
                                                                                                                                        agTAATCGACCAGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCACCAGAAcc            >  2:750841/1‑150 (MQ=255)
                                                                                                                                         ttGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACca           <  1:1126979/149‑1 (MQ=255)
                                                                                                                                         gTGATCGACCTGGTCGATTACTTACTTTCCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGGAACAGAGccc           >  2:1023709/1‑149 (MQ=255)
                                                                                                                                         gTGATCGACATGGTCGATTACTTACTTACCAGTAGCTTTTTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTATTGAAGCTTATTCCAGCTTATAATCAACACCCATGCGCTTAGACCCCAACCCCCAACAGAACca           >  2:1150622/1‑150 (MQ=255)
                                                                                                                                           gATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCCTGCGCTTATCCCCCAACACGCAACAGAACcaaa         >  1:1744495/1‑149 (MQ=255)
                                                                                                                                             tCGACCTTGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGAATATAATCAAAACCCAAGCGCTTATCCCCCAACACGCAACAGAACcaaaaa       >  2:791300/1‑150 (MQ=255)
                                                                                                                                                  cTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCGTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAgcg  >  2:968976/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCG  >  NC_000913/1681475‑1681770

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