Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,681,460 T→G M162R (ATG→AGG)  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,4600TG97.8% 449.0 / 1.3 136M162R (ATG→AGG) folMdihydromonapterin reductase, NADPH‑dependent; dihydrofolate reductase isozyme
Reads supporting (aligned to +/- strand):  ref base T (3/0);  new base G (63/68);  total (68/68)
Fisher's exact test for biased strand distribution p-value = 1.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-01

Only 100 of 138 total aligned reads displayed.

GCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATA  >  NC_000913/1681363‑1681556
                                                                                                 |                                                                                                
gccgccAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGa                                                                                                <  1:2625428/100‑1 (MQ=255)
 ccgccAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTTTTTTTTTTTTTTTCAAAGCGGCACTGGATAATAGGAc                                                                                               >  2:190287/1‑100 (MQ=255)
  cgccAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGAcc                                                                                              <  1:987581/100‑1 (MQ=255)
   gccAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGAccc                                                                                             <  1:323738/100‑1 (MQ=255)
   gccAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGAccc                                                                                             <  1:3593289/100‑1 (MQ=255)
       gCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGCAAGCTGCAAGCCAAGCGGCACTAGATAATACGACCCGccc                                                                                         >  2:1170008/1‑98 (MQ=255)
         gATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGt                                                                                       >  2:2897034/1‑100 (MQ=255)
                   aCTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCaa                                                                             >  2:1992971/1‑100 (MQ=255)
                   aCTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCaa                                                                             >  1:3314689/1‑100 (MQ=255)
                    cTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAg                                                                            >  2:1185899/1‑100 (MQ=255)
                     tttACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGc                                                                           <  1:1774767/100‑1 (MQ=255)
                        aCCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTgg                                                                        >  1:3710415/1‑100 (MQ=255)
                         gcGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGc                                                                       >  2:245983/2‑100 (MQ=255)
                         ccGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGc                                                                       >  1:308810/1‑100 (MQ=255)
                          cGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCa                                                                      >  1:1469998/1‑100 (MQ=255)
                          cGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCa                                                                      <  1:875562/100‑1 (MQ=255)
                          cGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCa                                                                      >  1:2243091/1‑100 (MQ=255)
                          cGATTATGGGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACGCGCTCGTTTGCCCGCCAGCTGGCa                                                                      >  2:1370136/1‑100 (MQ=255)
                           gATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGacac                                                                     >  1:994977/1‑100 (MQ=255)
                           gATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCAc                                                                     <  2:736212/100‑1 (MQ=255)
                            aTTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGTGGCACAGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGGAcc                                                                    >  2:1426336/1‑100 (MQ=255)
                             ttATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCg                                                                   <  2:3607184/100‑1 (MQ=255)
                               aTGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGa                                                                 <  1:3139348/100‑1 (MQ=255)
                                tgtgGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGaa                                                                >  2:1164060/1‑100 (MQ=255)
                                 gtggtgGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGaag                                                               >  1:2324354/1‑100 (MQ=255)
                                 gtggtgGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGaag                                                               <  1:2504197/100‑1 (MQ=255)
                                 gtggtgGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGaag                                                               <  2:496670/100‑1 (MQ=255)
                                 gtggtgGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGaag                                                               >  1:1337185/1‑100 (MQ=255)
                                 gtggtgGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGaag                                                               <  2:712010/100‑1 (MQ=255)
                                   ggtggAGCGCGGTAGCGACAAACATATTCCGTATGCTGCAAGCCAAGCGGCACTGGATAATAAGACCCGCTCGTTTGCCCGCCAGGGGGCACCgggactg                                                             >  2:2456964/1‑95 (MQ=255)
                                    gtggAGCGCGATAGCGAAGAAACTATTGCGTATGCTGCACGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCCAGATGGCACCGGaagtga                                                            >  2:3992106/1‑100 (MQ=255)
                                     tggAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGaagtgat                                                           <  2:812540/100‑2 (MQ=255)
                                     tggAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGaagtgaa                                                           >  2:2586729/1‑100 (MQ=255)
                                      ggAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGaagtgaaa                                                          >  1:3161784/1‑100 (MQ=255)
                                       gAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGaagtgaaag                                                         <  1:1980070/100‑1 (MQ=255)
                                        aGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGaagtgaaagt                                                        <  1:1069478/100‑1 (MQ=255)
                                         gcgcGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGaagtgaaagtg                                                       <  2:536/100‑1 (MQ=255)
                                          cgcgGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGaagtgaaagtga                                                      <  2:913238/100‑1 (MQ=255)
                                           gcgGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAagtgaaagtgaa                                                     <  1:203177/100‑1 (MQ=255)
                                           gcgGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAagtgaaagtgaa                                                     >  1:1689618/1‑100 (MQ=255)
                                            cgGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAAt                                                    <  1:1452390/100‑1 (MQ=255)
                                             ggTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAAtt                                                   >  2:790128/1‑100 (MQ=255)
                                             ggTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAAtt                                                   >  1:401313/1‑100 (MQ=255)
                                             ggTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAAtt                                                   <  2:3896293/100‑1 (MQ=255)
                                               tAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCt                                                 <  1:2920199/100‑1 (MQ=255)
                                                     cAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTgcg                                           >  1:2445087/1‑100 (MQ=255)
                                                     cAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTgcg                                           >  1:908256/1‑100 (MQ=255)
                                                      aaaCATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTgcgc                                          >  1:3110460/1‑100 (MQ=255)
                                                       aaCATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTAAAAGTGAATTCTATTGCGcc                                         >  1:3823477/1‑100 (MQ=255)
                                                          atatTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATc                                      >  1:194259/1‑100 (MQ=255)
                                                           tatTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCg                                     >  2:280019/1‑100 (MQ=255)
                                                           tatTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCg                                     >  1:1712244/1‑100 (MQ=255)
                                                                 gTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATc                               <  1:3440897/100‑1 (MQ=255)
                                                                 gTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATc                               <  1:1377870/100‑1 (MQ=255)
                                                                 gTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATc                               <  2:1952906/100‑1 (MQ=255)
                                                                 gTATGCTGCAAGCAAAGCGGCACTGGAGAATAGGACCCGCTCGTTTGCCCGCCAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATc                               >  2:564586/1‑100 (MQ=255)
                                                                 atctgctgcAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATc                               <  1:2913193/97‑1 (MQ=255)
                                                                  tATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATcc                              >  1:2790706/1‑100 (MQ=255)
                                                                    tgctgcAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTg                            <  2:1902706/100‑1 (MQ=255)
                                                                     gctgcAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATcagat                           >  2:2869489/1‑96 (MQ=255)
                                                                     gctgcAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGt                           >  1:1788482/1‑100 (MQ=255)
                                                                      ctgcAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGtt                          >  1:1278342/1‑100 (MQ=255)
                                                                      ctgcAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGtt                          >  2:2036616/1‑100 (MQ=255)
                                                                      ctgcAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGtt                          <  1:240265/100‑1 (MQ=255)
                                                                      ctgcAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGtt                          <  1:833379/100‑1 (MQ=255)
                                                                      ctgcAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGtt                          <  1:4039781/100‑1 (MQ=255)
                                                                       tgcAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGttt                         <  2:1898312/100‑1 (MQ=255)
                                                                        gcaagcaaAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTa                        <  1:1907431/100‑1 (MQ=255)
                                                                        gcaagcaaAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTa                        >  1:779858/1‑100 (MQ=255)
                                                                        gcaagcaaAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTa                        <  1:1865703/100‑1 (MQ=255)
                                                                         caagcacAGCGGCAGTAGATAATAGGACCCGCTCGTTTGCCGGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTaa                       <  2:2359053/100‑1 (MQ=255)
                                                                         caagcaaAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGCAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTaa                       >  2:2760490/1‑100 (MQ=255)
                                                                          aagcaaAGCGGCACTGGATAATCGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAAt                      <  1:3716037/100‑1 (MQ=255)
                                                                          aagcaaAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAAt                      <  2:916642/100‑1 (MQ=255)
                                                                          aagcaaAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAAt                      <  2:25382/100‑1 (MQ=255)
                                                                          aagcaaAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAAt                      <  2:3103442/100‑1 (MQ=255)
                                                                            gcaaAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGa                    <  2:2526794/100‑1 (MQ=255)
                                                                            gcaaAGCGGCACGGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGa                    <  1:570774/100‑1 (MQ=255)
                                                                             caaAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGaa                   >  1:1127987/1‑100 (MQ=255)
                                                                              aaaGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAAc                  >  1:2199496/1‑100 (MQ=255)
                                                                              aaaGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAAc                  >  1:1767052/1‑100 (MQ=255)
                                                                                   ggCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgat             <  2:1812378/100‑1 (MQ=255)
                                                                                   ggCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgat             <  1:3433594/100‑1 (MQ=255)
                                                                                   ggCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGAACCGGAAGTGAAAGTGAATTCTATTGCGCCATAGCTGATCCTGTTTAATGAACatgat             >  1:1728514/1‑100 (MQ=255)
                                                                                    gCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgatg            >  1:338192/1‑100 (MQ=255)
                                                                                     cACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgatga           >  1:1121283/1‑100 (MQ=255)
                                                                                     cACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgatga           <  1:705899/100‑1 (MQ=255)
                                                                                     cACTGGATAATAGGACCCGCTCGTTTACACGCAAGCTGGCACCGGAAGTGAATTTGAATTCTATTGCGCCATCGCTGTTCCTGTTTAATGAACatgatga           <  2:2113099/100‑1 (MQ=255)
                                                                                      aCTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgatgat          <  2:3397784/100‑1 (MQ=255)
                                                                                      aCTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgatgat          >  1:3808120/1‑100 (MQ=255)
                                                                                       cTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgatgatg         <  2:1486298/100‑1 (MQ=255)
                                                                                         ggATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGcc       >  2:3894809/1‑100 (MQ=255)
                                                                                         ggATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGcc       >  1:3304182/1‑100 (MQ=255)
                                                                                          gATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCg      >  2:2064614/1‑100 (MQ=255)
                                                                                          gATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCg      <  2:1700168/100‑1 (MQ=255)
                                                                                          gATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCg      <  1:2701433/100‑1 (MQ=255)
                                                                                          gATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCg      <  1:2646933/100‑1 (MQ=255)
                                                                                          gATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCg      >  1:1822654/1‑100 (MQ=255)
                                                                                           ataataGGACCCGCTCGTTTGCCCGCAATCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGa     >  1:3211338/1‑100 (MQ=255)
                                                                                              atagGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAata  >  2:925595/1‑100 (MQ=255)
                                                                                                 |                                                                                                
GCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATA  >  NC_000913/1681363‑1681556

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