Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,681,679 G→A G235D (GGC→GAC)  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,6790GA98.6% 286.6 / ‑4.3 73G235D (GGC→GAC) folMdihydromonapterin reductase, NADPH‑dependent; dihydrofolate reductase isozyme
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (31/41);  total (31/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

TTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATA  >  NC_000913/1681531‑1681826
                                                                                                                                                    |                                                                                                                                                   
tttAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACg                                                                                                                                                   >  1:954437/1‑151 (MQ=255)
   aaTGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGgtc                                                                                                                                                >  1:622401/1‑151 (MQ=255)
    aTGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGgtcg                                                                                                                                               >  1:1016581/1‑151 (MQ=255)
       aaCATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGgtcgtc                                                                                                                                             >  2:1078710/1‑150 (MQ=255)
       aaCATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCa                                                                                                                                            <  1:599009/151‑1 (MQ=255)
         cATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATc                                                                                                                                          <  2:1197876/151‑1 (MQ=255)
            gatgatgCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGc                                                                                                                                       <  1:151556/151‑1 (MQ=255)
               gatgCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGtt                                                                                                                                    <  2:519263/151‑1 (MQ=255)
                  gCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAAt                                                                                                                                 <  2:1043434/151‑1 (MQ=255)
                  gCCGAATATCGACAACAGGCGCAGACTAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAAt                                                                                                                                 >  2:180019/1‑151 (MQ=255)
                   ccGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATg                                                                                                                                <  2:779639/151‑1 (MQ=255)
                     gAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGc                                                                                                                               >  2:187369/1‑150 (MQ=255)
                          tCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCttt                                                                                                                         >  1:600812/1‑151 (MQ=255)
                               aacaGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATcc                                                                                                                     >  1:976459/1‑150 (MQ=255)
                                acaGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATcc                                                                                                                     >  1:47890/1‑149 (MQ=255)
                                 caGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGc                                                                                                                  <  2:1188249/151‑1 (MQ=255)
                                      gcTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCAtata                                                                                                             <  1:942409/151‑1 (MQ=255)
                                       cTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCAtata                                                                                                             >  2:302190/1‑150 (MQ=255)
                                           aTAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCa                                                                                                         >  1:1136503/1‑150 (MQ=255)
                                             aaaTCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAaca                                                                                                      <  1:830641/151‑1 (MQ=255)
                                             aaaTCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAac                                                                                                       >  1:708951/1‑150 (MQ=255)
                                               aTCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACAcc                                                                                                    <  1:455653/151‑1 (MQ=255)
                                                 cACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCa                                                                                                  <  1:510748/151‑1 (MQ=255)
                                                    tgatgaAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATgc                                                                                               >  2:413483/1‑151 (MQ=255)
                                                     gatgaAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATgc                                                                                               <  2:256373/150‑1 (MQ=255)
                                                     gatgaAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATg                                                                                                <  1:562931/149‑1 (MQ=255)
                                                       tgaAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCt                                                                                            >  2:374542/1‑151 (MQ=255)
                                                          aaaCCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTAt                                                                                         <  1:496511/151‑1 (MQ=255)
                                                           aaCCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATc                                                                                        >  1:965215/1‑151 (MQ=255)
                                                             ccGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGCTTACTTACTTACCCGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATccc                                                                                      >  1:277422/1‑151 (MQ=255)
                                                             ccGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATccc                                                                                      <  2:344739/151‑1 (MQ=255)
                                                               gcgcCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATccccc                                                                                    >  2:393615/1‑151 (MQ=255)
                                                                 gcCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCaa                                                                                  <  1:684056/151‑1 (MQ=255)
                                                                  ccTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAac                                                                                 <  2:272277/151‑1 (MQ=255)
                                                                     ggCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACg                                                                              <  2:798581/151‑1 (MQ=255)
                                                                      gCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGc                                                                             <  2:281317/151‑1 (MQ=255)
                                                                       cGAGAAAGAAGTGATCGACCTGGTCGATCACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGc                                                                             >  1:629993/1‑150 (MQ=255)
                                                                        gagaAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCaa                                                                           <  2:1118108/151‑1 (MQ=255)
                                                                        gagaAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCaa                                                                           <  2:466374/151‑1 (MQ=255)
                                                                         agaaagaaGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAAc                                                                          >  2:234230/1‑151 (MQ=255)
                                                                              gaagTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAAc                                                                     >  1:1150217/1‑151 (MQ=255)
                                                                                agTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAAcc                                                                    >  1:543511/1‑150 (MQ=255)
                                                                                  tGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACcaac                                                                 >  1:598490/1‑151 (MQ=255)
                                                                                   gATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACcaaca                                                                <  2:1019443/151‑1 (MQ=255)
                                                                                        aCCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAgc                                                           >  2:484152/1‑151 (MQ=255)
                                                                                         ccTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAgcg                                                          >  2:1093562/1‑151 (MQ=255)
                                                                                         ccTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAgcg                                                          <  1:627078/151‑1 (MQ=255)
                                                                                          cTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAgcgg                                                         >  1:405109/1‑151 (MQ=255)
                                                                                            ggTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAgcggc                                                        <  1:236781/150‑1 (MQ=255)
                                                                                               cGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCggg                                                     >  1:96167/1‑150 (MQ=255)
                                                                                                   tacttacttacCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAgc                                                >  1:833576/1‑151 (MQ=255)
                                                                                                   tacttacttacCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAg                                                 <  2:447621/150‑1 (MQ=255)
                                                                                                    acttacttacCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAgcg                                               <  1:437284/151‑1 (MQ=255)
                                                                                                    acttacttacCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAgc                                                <  2:105476/150‑1 (MQ=255)
                                                                                                    acttacttacCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAgc                                                <  1:437654/150‑1 (MQ=255)
                                                                                                       tacttacCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAgcgc                                              <  1:298642/149‑1 (MQ=255)
                                                                                                              cAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCg                                     <  1:1078683/151‑1 (MQ=255)
                                                                                                                   gCTTTGTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGa                                 <  2:817431/150‑1 (MQ=255)
                                                                                                                        gTCACCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATAcc                           >  2:212051/1‑151 (MQ=255)
                                                                                                                           aCCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACcaca                        <  1:302428/151‑1 (MQ=255)
                                                                                                                           aCCGGACGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACcaca                        <  2:509203/151‑1 (MQ=255)
                                                                                                                                aCGCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGa                    >  2:1052201/1‑150 (MQ=255)
                                                                                                                                  gCAGTTTCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGacac                 <  2:869011/151‑1 (MQ=255)
                                                                                                                                      tttCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAgg             >  2:281916/1‑151 (MQ=255)
                                                                                                                                       ttCCCACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAgg             >  1:155564/1‑150 (MQ=255)
                                                                                                                                         cccACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACTCCAGGTa           <  2:474369/150‑1 (MQ=255)
                                                                                                                                         cccACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTa           <  2:224394/150‑1 (MQ=255)
                                                                                                                                         cccACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTa           <  1:19687/150‑1 (MQ=255)
                                                                                                                                         cccACTTGATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTa           <  2:702105/150‑1 (MQ=255)
                                                                                                                                                gATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACaaaa    <  1:652696/150‑1 (MQ=255)
                                                                                                                                                gATGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAAt   <  2:1016605/151‑1 (MQ=255)
                                                                                                                                                 aTGACGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATa  <  1:797492/151‑1 (MQ=255)
                                                                                                                                                 aTGACGGTCGTCATCTACGTTAATGCAGCTTTAGCCAGCATATAACCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATACCGGTGAAAGACCACAGCGACACCAGGTAAGCaaaa    <  1:1084361/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATA  >  NC_000913/1681531‑1681826

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