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,4600TG99.5% 754.3 / ‑4.3 212M162R (ATG→AGG) folMdihydromonapterin reductase, NADPH‑dependent; dihydrofolate reductase isozyme
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (107/104);  total (108/104)
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.02e-01

Only 100 of 214 total aligned reads displayed.

GCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATA  >  NC_000913/1681363‑1681556
                                                                                                 |                                                                                                
gccgccAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGa                                                                                                <  2:25125/100‑1 (MQ=255)
 ccgccAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGAc                                                                                               <  1:3328797/100‑1 (MQ=255)
 ccgccAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGAc                                                                                               <  1:2792699/100‑1 (MQ=255)
  cgccAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGAcc                                                                                              <  1:1507822/100‑1 (MQ=255)
   gccAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGAccc                                                                                             <  1:260705/100‑1 (MQ=255)
   gccAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGAATGCTGCAAGCAAAGCGGCACTGGATAATAGGAccc                                                                                             >  2:4054988/1‑100 (MQ=255)
    ccAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCg                                                                                            >  1:1928744/1‑100 (MQ=255)
     cAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGc                                                                                           <  2:643538/100‑1 (MQ=255)
      aGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCt                                                                                          <  2:2211899/100‑1 (MQ=255)
      aGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCt                                                                                          <  1:2147922/100‑1 (MQ=255)
      aGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCt                                                                                          >  2:4106571/1‑100 (MQ=255)
        cGATATGATTCACTTTACCGATTAGGTGGTTGATCGCGGTAGCGATAAACATATTGCGTCTGCAGCAAGCAAAGCGGCACTGGATAATAGGACCCTCTCg                                                                                        <  2:5167404/100‑1 (MQ=255)
        cGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCg                                                                                        >  1:1971878/1‑100 (MQ=255)
        cGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCg                                                                                        >  1:531863/1‑100 (MQ=255)
         gATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGt                                                                                       >  1:3982661/1‑100 (MQ=255)
         gATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGt                                                                                       <  1:4703216/100‑1 (MQ=255)
            atCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTg                                                                                    <  1:3110173/100‑1 (MQ=255)
            atCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTg                                                                                    <  2:1513424/100‑1 (MQ=255)
              cATTCCCTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGcc                                                                                  >  2:4845951/1‑100 (MQ=255)
              cATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGcc                                                                                  >  1:5140398/1‑100 (MQ=255)
              cATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGcc                                                                                  <  1:3075546/100‑1 (MQ=255)
              cATCCACTTCTCCGATTATGTGGTGGAGCGCGGCAGCGACAAACATATTGCTTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGcc                                                                                  <  2:4629042/100‑1 (MQ=255)
               aTTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAACGCGGAACTGGATAATAGGACACGCTCGTTTGcct                                                                                 >  2:2756112/1‑99 (MQ=255)
               aTTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGccc                                                                                 >  2:1216908/1‑100 (MQ=255)
                ttCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCg                                                                                <  2:3436280/100‑1 (MQ=255)
                 tCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGc                                                                               >  2:3491308/1‑100 (MQ=255)
                  cACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCa                                                                              >  1:2008644/1‑100 (MQ=255)
                    cTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAg                                                                            >  2:831949/1‑100 (MQ=255)
                       tACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTg                                                                         <  2:4562496/100‑1 (MQ=255)
                         ccGATTATGTGGTGGAGCTCGTGAGCGACAATCATATTGCGTATGCTGCAATCAAAGCTGCACTTGATAAGAGTACCCGCTCGTGTGCGCTGGAGCTGGc                                                                       >  2:3024712/1‑100 (MQ=255)
                         ccGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGc                                                                       >  1:2508257/1‑100 (MQ=255)
                          cGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCa                                                                      >  1:1752697/1‑100 (MQ=255)
                           gATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCAc                                                                     >  1:4532068/1‑100 (MQ=255)
                                tgtgGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTAGTTTGCCCGCAAGCTGGCACCGGaa                                                                >  1:2782553/1‑100 (MQ=255)
                                tgtgGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAACGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGaa                                                                <  1:5032283/100‑1 (MQ=255)
                                 gtggtgGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGaag                                                               >  2:357209/1‑100 (MQ=255)
                                  tggtggAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGaagt                                                              >  1:957795/1‑100 (MQ=255)
                                  tggtggAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGaagt                                                              >  2:3862126/1‑100 (MQ=255)
                                  tggtggAGCGCGGTAGCGACAAACATATCGCGTATGCTGCAAGCAAAGCCGCACTCGATAATAGGACCCACTCGTTTGCCCGCAAGCTGGCCGCGGaagt                                                              >  2:1605235/1‑100 (MQ=255)
                                   ggtggAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGaagtg                                                             >  1:209209/1‑100 (MQ=255)
                                   ggtggAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAACAAAGCGACACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGaagtg                                                             >  1:5291068/1‑100 (MQ=255)
                                    gtggAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGaagtga                                                            <  2:5349833/100‑1 (MQ=255)
                                     tggAGCGCGGTAGCGACAAACATACTTCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGaagtgaa                                                           <  1:4572388/100‑1 (MQ=255)
                                       gAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGaagtgaaag                                                         >  2:3977947/1‑100 (MQ=255)
                                         gcgcGGTAGCGACAAACATTTTGCGTATGCTGCAAGCAAAGCGGCACTGGATACTAGGACCTGCTCGCTTGCCCGCAAGCTGGCACCGGaagtgaaagtg                                                       <  1:4954258/100‑1 (MQ=255)
                                         gcgcGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGaagtgaaagtg                                                       <  2:60975/100‑1 (MQ=255)
                                               tAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCt                                                 >  2:1292354/1‑100 (MQ=255)
                                                aGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTa                                                >  1:4495846/1‑100 (MQ=255)
                                                aGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTa                                                <  1:3897760/100‑1 (MQ=255)
                                                 gCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTAt                                               <  1:2569212/100‑1 (MQ=255)
                                                 gCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTAt                                               >  1:2240476/1‑100 (MQ=255)
                                                   ggcAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTg                                             >  2:3660837/3‑100 (MQ=255)
                                                     cAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTgcg                                           >  1:4749343/1‑100 (MQ=255)
                                                       aaCATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGcc                                         <  2:3330754/100‑1 (MQ=255)
                                                       aaCATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGcc                                         <  1:560859/100‑1 (MQ=255)
                                                       aaCATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGcc                                         >  1:3792556/1‑100 (MQ=255)
                                                        aCATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCa                                        >  2:2286361/1‑100 (MQ=255)
                                                         cATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCAt                                       >  1:200475/1‑100 (MQ=255)
                                                         cATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCAt                                       >  1:1653316/1‑100 (MQ=255)
                                                           tatTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCg                                     <  2:3373845/100‑1 (MQ=255)
                                                               gCGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGa                                 >  2:923789/1‑100 (MQ=255)
                                                               gCGTATGCTGCAAGCAAAGCGGCACTGGAAAATAGGACCCGCTCGTTTGCACGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCTCTGa                                 >  2:1113970/1‑100 (MQ=255)
                                                                cGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGAt                                <  2:173235/100‑1 (MQ=255)
                                                                cGTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGAt                                <  1:3292363/100‑1 (MQ=255)
                                                                 gTATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATc                               <  2:2046717/100‑1 (MQ=255)
                                                                  tATGCTGCAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATcc                              >  2:428510/1‑100 (MQ=255)
                                                                  tATGCTGCAAGCAAAGCGGCACCGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATcc                              >  1:4845327/1‑100 (MQ=255)
                                                                    tgctgaAAGCCAAGCGGCACTGGATAACAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTg                            <  1:3053139/100‑1 (MQ=255)
                                                                     gctgcAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGt                           >  2:1034571/1‑100 (MQ=255)
                                                                      ctgcAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGTAAGCTGGCACCGGAAGTGTAAGTGAATTCTATTGCGCCATCGCTGATCCTGtt                          <  2:5069643/100‑1 (MQ=255)
                                                                      ctgcAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGtt                          <  1:22856/100‑1 (MQ=255)
                                                                      ctgcAAGCAAAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGtt                          <  2:3044178/100‑1 (MQ=255)
                                                                       tgcAAGCAAAGGGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGttt                         <  2:525840/100‑1 (MQ=255)
                                                                        gcaagcaaAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTa                        >  1:4615312/1‑100 (MQ=255)
                                                                         caagcaaAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTaa                       <  2:2637091/100‑1 (MQ=255)
                                                                         caagcaaAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTaa                       <  2:2369679/100‑1 (MQ=255)
                                                                          aagcaaAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAAt                      <  1:3581310/100‑1 (MQ=255)
                                                                           agcaaAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATg                     <  1:1684874/100‑1 (MQ=255)
                                                                            gcaaAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAGTGa                    <  1:4282661/100‑1 (MQ=255)
                                                                            gcaaAGCGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGa                    <  1:2733237/100‑1 (MQ=255)
                                                                                  cGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatga              <  1:5372831/100‑1 (MQ=255)
                                                                                  cGGCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatga              <  2:2030178/100‑1 (MQ=255)
                                                                                   ggCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgat             <  1:113457/100‑1 (MQ=255)
                                                                                   ggCACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgat             <  2:1423832/100‑1 (MQ=255)
                                                                                     cACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgatga           >  1:2951265/1‑100 (MQ=255)
                                                                                     cACTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgatga           >  1:5156430/1‑100 (MQ=255)
                                                                                     cACCGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgatga           <  2:4776877/100‑1 (MQ=255)
                                                                                      aCTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgatgat          <  1:4895874/100‑1 (MQ=255)
                                                                                       cTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgatgatg         >  1:5479775/1‑100 (MQ=255)
                                                                                       cTGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgatgatg         >  1:1232518/1‑100 (MQ=255)
                                                                                        tGGATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGc        >  1:2071559/1‑100 (MQ=255)
                                                                                         ngATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGcc       >  1:480533/2‑100 (MQ=255)
                                                                                         ggATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGcc       >  1:5145882/1‑100 (MQ=255)
                                                                                          gATAATAGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCg      <  1:4144003/100‑1 (MQ=255)
                                                                                           ataataGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTAAATTCTATTGAGCCATCGCTTATCCTGTTTAATGAACATGATGATGCCGa     <  2:1955688/100‑1 (MQ=255)
                                                                                            taataGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGaa    <  1:3357237/100‑1 (MQ=255)
                                                                                            taataGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGaa    >  2:396171/1‑100 (MQ=255)
                                                                                            taataGGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGaa    >  1:4040981/1‑100 (MQ=255)
                                                                                              atagGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAata  >  1:3158705/1‑100 (MQ=255)
                                                                                              atagGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTATTTAATGAACATGATGATGCCGAata  >  2:4424341/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: