Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,114,539 A→G 6.3% intergenic (‑67/+30) fpr ← / ← glpX ferredoxin‑NADP reductase; flavodoxin reductase/fructose 1,6‑bisphosphatase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,114,5390AG 6.3% 422.6 / 13.2 143intergenic (‑67/+30)fpr/glpXferredoxin‑NADP reductase; flavodoxin reductase/fructose 1,6‑bisphosphatase II
Reads supporting (aligned to +/- strand):  ref base A (73/61);  new base G (6/3);  total (79/64)
Fisher's exact test for biased strand distribution p-value = 7.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

Only 100 of 146 total aligned reads displayed.

CGGGGGCGTGAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGATATTGCC  >  NC_000913/4114391‑4114688
                                                                                                                                                    |                                                                                                                                                     
      cgggAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTgg                                                                                                                                                <  2:723309/147‑1 (MQ=255)
        tGAACGGTGAGACTAAACAGGGCGGCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGGTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGaa                                                                                                                                              <  1:1289328/150‑1 (MQ=255)
         gAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAg                                                                                                                                             >  1:760919/1‑150 (MQ=255)
         gAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAg                                                                                                                                             <  1:1268887/150‑1 (MQ=255)
          aaCGGTGAGACTAAACAGGGCGTCTGTCCAGTTCTGCACTTTAGTGACTATTCCTGTTACCCAATCAGCCATCTTTATCTCCTGTATAGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCTTCGGGACTGGaatt                                                                                                                                            <  1:282192/150‑3 (MQ=255)
          aaCGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCTTCGGGACTGGaatt                                                                                                                                            <  1:1209076/150‑3 (MQ=255)
          aaCGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAgg                                                                                                                                            <  1:1366367/150‑1 (MQ=255)
          aaCGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAgg                                                                                                                                            <  1:1159519/150‑1 (MQ=255)
           aCGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAATGACTGGAAGGc                                                                                                                                           >  2:1578714/1‑150 (MQ=255)
             ggTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTc                                                                                                                                         <  1:976825/150‑1 (MQ=255)
             ggTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCCTCGTTCTCCGAACAAAAAAATTCCAGTCCTTCGGGACTGGaattttt                                                                                                                                         <  1:1119755/150‑6 (MQ=255)
              gTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCa                                                                                                                                        >  1:267028/1‑150 (MQ=255)
                  gaCTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCg                                                                                                                                    <  2:786710/150‑1 (MQ=255)
                    cTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGat                                                                                                                                  >  2:1710141/1‑150 (MQ=255)
                     tAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTACGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGatc                                                                                                                                 >  2:1053618/1‑150 (MQ=255)
                      aaaCAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGatca                                                                                                                                <  1:691299/150‑1 (MQ=255)
                                  ggTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCagag                                                                                                                    <  1:201958/150‑1 (MQ=255)
                                   gTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAgg                                                                                                                   >  1:475758/1‑150 (MQ=255)
                                     ccAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGAt                                                                                                                 >  1:10483/1‑150 (MQ=255)
                                      cAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGAtg                                                                                                                <  2:1469967/150‑1 (MQ=255)
                                         ttCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGc                                                                                                             <  2:548959/150‑1 (MQ=255)
                                           cTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCAc                                                                                                           <  2:112334/150‑1 (MQ=255)
                                             gCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCt                                                                                                         >  1:349411/1‑150 (MQ=255)
                                                cTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCa                                                                                                      >  1:848275/1‑150 (MQ=255)
                                                 tttAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCAt                                                                                                     >  1:786501/1‑150 (MQ=255)
                                                     gTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTcc                                                                                                 >  1:1235837/1‑150 (MQ=255)
                                                      tGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCg                                                                                                <  2:1695126/150‑1 (MQ=255)
                                                        aCTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCggg                                                                                              <  2:308518/150‑1 (MQ=255)
                                                        aCTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCggg                                                                                              >  1:886817/1‑150 (MQ=255)
                                                        aCTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCggg                                                                                              <  2:349411/150‑1 (MQ=255)
                                                        aCTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCggg                                                                                              <  2:1424721/150‑1 (MQ=255)
                                                        aCTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCggg                                                                                              <  2:220642/150‑1 (MQ=255)
                                                         cTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGt                                                                                             >  2:648289/1‑150 (MQ=255)
                                                             gCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCttt                                                                                         <  2:50481/150‑1 (MQ=255)
                                                             gCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCttt                                                                                         <  2:1514259/150‑1 (MQ=255)
                                                              ccTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTg                                                                                        >  1:1658799/1‑150 (MQ=255)
                                                              ccTGTTACCCAATCAGCCATGCTTTTCACCTGTTTTGAATGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGGGCACCTGCATTTCCGGGTCTTTg                                                                                        <  2:1215717/150‑1 (MQ=255)
                                                                 gTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGa                                                                                     <  1:1175673/150‑1 (MQ=255)
                                                                 gTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGATCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGa                                                                                     >  2:1277546/1‑150 (MQ=255)
                                                                  ttACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGAt                                                                                    <  2:760919/150‑1 (MQ=255)
                                                                   tACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATc                                                                                   <  2:848275/150‑1 (MQ=255)
                                                                     cccAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGATGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGTGATCCa                                                                                 <  1:860934/150‑1 (MQ=255)
                                                                     cccAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCa                                                                                 <  1:1710141/150‑1 (MQ=255)
                                                                           cAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCACGAAGGACTGGAAGGCTCAATCGATCAAAAAAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCATATAGt                                                                           >  1:1239379/1‑150 (MQ=255)
                                                                             gCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGtgg                                                                         >  2:1571680/1‑150 (MQ=255)
                                                                             gCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGtgg                                                                         >  2:134761/1‑150 (MQ=255)
                                                                              ccATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGtgga                                                                        >  2:1490617/1‑150 (MQ=255)
                                                                               cATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGATCAAAAAAATTCCAGTCCCTAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGtggat                                                                       >  1:1149278/1‑150 (MQ=255)
                                                                                aTGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGtggatg                                                                      <  2:465461/150‑1 (MQ=255)
                                                                                 tGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGtggatgg                                                                     >  1:863821/1‑150 (MQ=255)
                                                                                  gTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGtggatgga                                                                    >  1:910083/1‑150 (MQ=255)
                                                                                       tctcCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGaa                                                               >  1:538952/1‑150 (MQ=255)
                                                                                             gTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATgcggc                                                         <  1:144812/150‑1 (MQ=255)
                                                                                                 tgattgaCTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAAt                                                     >  2:398224/1‑150 (MQ=255)
                                                                                                 tgattgaCTTGACTTATCTACGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAACGCTCAAACGATCAAATCAATCAGAGGAAGTGCACCTGCATATCCGGGACTTTGCGATACAGATAGTGGATGGACTGAATGCGGCGAAt                                                     >  2:1604761/1‑150 (MQ=255)
                                                                                                  gattgaCTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATg                                                    <  2:10483/150‑1 (MQ=255)
                                                                                                       aCTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATgcggc                                                         <  2:689365/140‑1 (MQ=255)
                                                                                                       aCTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATgcggc                                                         >  1:689365/1‑140 (MQ=255)
                                                                                                       aCTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCg                                               <  1:1512586/150‑1 (MQ=255)
                                                                                                       aCTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCg                                               <  2:267028/150‑1 (MQ=255)
                                                                                                        cTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGt                                              >  1:1069640/1‑150 (MQ=255)
                                                                                                          tGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGa                                            <  1:1330148/150‑1 (MQ=255)
                                                                                                                aTCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGcc                                      >  2:455523/1‑150 (MQ=255)
                                                                                                                 tCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCcg                                     <  1:1053618/150‑1 (MQ=255)
                                                                                                                 tCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCcg                                     >  2:70590/1‑150 (MQ=255)
                                                                                                                 tCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCAGCGAATGGTGCGTGACTTGCcg                                     <  2:927133/150‑1 (MQ=255)
                                                                                                                  cTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCcgc                                    >  1:348919/1‑150 (MQ=255)
                                                                                                                   ttCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCcgcg                                   >  1:1224393/1‑150 (MQ=255)
                                                                                                                     cGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGa                                 >  1:1592245/1‑150 (MQ=255)
                                                                                                                      gTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGAt                                >  1:1375849/1‑150 (MQ=255)
                                                                                                                      gTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGAGTGACTTGCCGCGGAt                                >  2:101556/1‑150 (MQ=255)
                                                                                                                       ttCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATgcgg                                                          <  2:373381/123‑1 (MQ=255)
                                                                                                                         ctcCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATcag                             >  2:145704/1‑150 (MQ=255)
                                                                                                                           ccGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATcagca                           <  1:395605/150‑1 (MQ=255)
                                                                                                                           ccGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATcagca                           <  1:1277023/150‑1 (MQ=255)
                                                                                                                           cagAACAAAAAAATTCCAGTCCCGAAGGACTGGAATGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCTGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATcagca                           >  2:589174/3‑150 (MQ=255)
                                                                                                                            cGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCTGATcagcag                          >  2:541701/1‑150 (MQ=255)
                                                                                                                            cGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATcagcag                          <  2:1224393/150‑1 (MQ=255)
                                                                                                                             gAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCagcagc                         <  2:1112452/150‑1 (MQ=255)
                                                                                                                                cAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGtt                      >  2:1162501/1‑150 (MQ=255)
                                                                                                                                  aaaaaaTTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTc                    <  2:863821/150‑1 (MQ=255)
                                                                                                                                  aaaaaaTTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTc                    >  2:744621/1‑150 (MQ=255)
                                                                                                                                    aaaaTTCCAGTCCCGAATGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATTGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCgg                  <  2:1239379/150‑1 (MQ=255)
                                                                                                                                    aaaaTTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCgg                  >  2:1607585/1‑150 (MQ=255)
                                                                                                                                    aaaaTTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCgg                  <  1:1571680/150‑1 (MQ=255)
                                                                                                                                     aaaTTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGt                 <  2:1426258/150‑1 (MQ=255)
                                                                                                                                       cttCCAGTCCTTCGGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAg               >  2:45109/2‑150 (MQ=255)
                                                                                                                                       cttCCAGTCCTTCGGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAg               >  1:966175/2‑150 (MQ=255)
                                                                                                                                       aTTACAGTACCGAAGGACTGGAAGGCTCAATCGATCAAATCATTCAGAGTATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATTGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAg               <  2:1115873/150‑1 (MQ=255)
                                                                                                                                        ttCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGt              <  1:398224/150‑1 (MQ=255)
                                                                                                                                          ccAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATAAGAGGATGTGCACCTGCATTTCCGGGTATTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCAGCGGATCAGCAGCGTTTCGGTAGTcg            >  1:948327/1‑150 (MQ=255)
                                                                                                                                            aGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTcgcg          >  1:1461705/1‑150 (MQ=255)
                                                                                                                                               cccGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGAGata       >  1:110897/1‑150 (MQ=255)
                                                                                                                                                ccGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCCCCTGCATTTCCGGGTCTTTGCGATCCAGATTGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGatat      >  2:243388/1‑150 (MQ=255)
                                                                                                                                                ccGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCTGTAGTCGCGatat      >  2:293328/1‑150 (MQ=255)
                                                                                                                                                ccGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGatat      >  1:1447379/1‑150 (MQ=255)
                                                                                                                                                ccGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGatat      >  2:1415965/1‑150 (MQ=255)
                                                                                                                                                    aGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGATATTGcc  <  2:381026/150‑1 (MQ=255)
                                                                                                                                                    aGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGATATTGcc  <  2:968871/150‑1 (MQ=255)
                                                                                                                                                    aGGACTGGAAGGCTCAAACGATCAAATCAATGAGAGGTTGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGAGGGACTGAATCCGGCGAATGGAGCGTGACTTGCCGCGGATCAGCAGCGGTTCGGTAGTCGCGATATTGcc  <  2:63339/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGGGGGCGTGAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGATATTGCC  >  NC_000913/4114391‑4114688

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