Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,114,539 A→G 6.9% 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.9% 383.7 / 14.5 132intergenic (‑67/+30)fpr/glpXferredoxin‑NADP reductase; flavodoxin reductase/fructose 1,6‑bisphosphatase II
Reads supporting (aligned to +/- strand):  ref base A (65/55);  new base G (4/5);  total (69/63)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.41e-01

Only 100 of 135 total aligned reads displayed.

CGGGGGCGTGAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGATATTGCC  >  NC_000913/4114391‑4114688
                                                                                                                                                     |                                                                                                                                                     
   gggCGTGAACGGTGATACTAATCAGGGCGTTGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCTAACAAAAAAAATCCAGTCCCGAAGGAc                                                                                                                                                   <  2:1005033/150‑1 (MQ=255)
   gggCGTGAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGAc                                                                                                                                                   <  1:324187/150‑1 (MQ=255)
       gTGAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCTTCGGGACTGGa                                                                                                                                               <  1:412198/150‑1 (MQ=255)
          aaCGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACACAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCACCGAACAAAAAAATTCCAGTCCTTCGGGACTGGaatt                                                                                                                                            <  2:57362/150‑3 (MQ=255)
           aCGGTGAGACTAAACAGGTCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATTTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGc                                                                                                                                           <  2:1407479/150‑1 (MQ=255)
            cGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCTTCGGGACTGGaatttt                                                                                                                                          <  1:1291195/150‑5 (MQ=255)
            cGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCTTCGGGACTGGaatttt                                                                                                                                          <  1:864536/150‑5 (MQ=255)
             ggTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTc                                                                                                                                         >  1:1724369/1‑150 (MQ=255)
              gTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCTTCGGGACTGGaatttttt                                                                                                                                        <  1:979348/150‑7 (MQ=255)
                  gaCTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCg                                                                                                                                    >  2:706589/1‑150 (MQ=255)
                  gaCTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCg                                                                                                                                    >  2:367630/1‑150 (MQ=255)
                     tAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGatc                                                                                                                                 <  2:56773/150‑1 (MQ=255)
                     tAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGatc                                                                                                                                 <  2:1457612/150‑1 (MQ=255)
                     tAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGatc                                                                                                                                 <  2:1099166/150‑1 (MQ=255)
                      aaaCAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGatca                                                                                                                                <  1:87103/150‑1 (MQ=255)
                         cAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGatcaaat                                                                                                                             <  1:1569390/150‑1 (MQ=255)
                           gggCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGatcaaatca                                                                                                                           <  2:1208612/150‑1 (MQ=255)
                                  ggTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCagag                                                                                                                    >  2:165664/1‑150 (MQ=255)
                                   gTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAgg                                                                                                                   >  1:1269439/1‑150 (MQ=255)
                                   gTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAgg                                                                                                                   >  2:187956/1‑150 (MQ=255)
                                    tCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGa                                                                                                                  <  1:6873/150‑1 (MQ=255)
                                    tCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGa                                                                                                                  <  2:529041/150‑1 (MQ=255)
                                      cAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGAtg                                                                                                                >  2:156190/1‑150 (MQ=255)
                                               aCTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGc                                                                                                       <  2:677165/150‑1 (MQ=255)
                                                cTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCACGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGAACCTGCa                                                                                                      >  1:371848/1‑150 (MQ=255)
                                                              ccTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTg                                                                                        >  2:120302/1‑150 (MQ=255)
                                                              ccTGTTACCCAATCACGCATGTTGATCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTGTCCTAACAAAACAATTCGAGCCCCGAAGGACTGGAAGGCTCAATCGATCACATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTg                                                                                        <  1:940593/150‑1 (MQ=255)
                                                               cTGTTTCACATTCAGCCATGTTTTTCTCCTGTTTTGATTTACTTGCCTTATCTTCGTTCTCCTAACAAAAAAAATCCAGTCACTAAGGACTGGAAGGCTCAATCGTTCAAATCAATCAGAGGATGTGCTCCTGCATTTCCGGGTCTTTGc                                                                                       <  2:265781/150‑1 (MQ=255)
                                                               cTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTATGc                                                                                       >  2:226710/1‑150 (MQ=255)
                                                                 gTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCa                                                                                                      <  1:907943/133‑1 (MQ=255)
                                                                   tACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATc                                                                                   >  1:1060434/1‑150 (MQ=255)
                                                                   tACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGACCCTGCATTTCCGGGTCTTTGAGATc                                                                                   >  2:260301/1‑150 (MQ=255)
                                                                    aCCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATcc                                                                                  <  1:466281/150‑1 (MQ=255)
                                                                     cccAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCa                                                                                 >  1:1454400/1‑150 (MQ=255)
                                                                     cccAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGAATGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCa                                                                                 <  1:270102/150‑1 (MQ=255)
                                                                       cAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAg                                                                                <  2:1604334/150‑1 (MQ=255)
                                                                       cAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGa                                                                               >  1:919988/1‑150 (MQ=255)
                                                                             gCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGtgg                                                                         >  2:1480450/1‑150 (MQ=255)
                                                                             gCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGtgg                                                                         >  2:477527/1‑150 (MQ=255)
                                                                             gCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCACGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGtgg                                                                         >  1:59960/1‑150 (MQ=255)
                                                                             gCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTAGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGtgg                                                                         <  2:1194642/150‑1 (MQ=255)
                                                                               cATGTTTTTCTCGTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGTACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGtggat                                                                       <  1:1664742/150‑1 (MQ=255)
                                                                               cATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGtggat                                                                       >  2:1696143/1‑150 (MQ=255)
                                                                                aTGTTTTTCTCCTGTTTTGTTTGACTTGCCTTATCTTCGTTCTACGAACAAAAAAATTCCAGTCCCGACGGACTGGAAGGCTCATTCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGtggatg                                                                      <  1:1449130/150‑1 (MQ=255)
                                                                                 tGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGtggatgg                                                                     >  1:826892/1‑150 (MQ=255)
                                                                                   tttttCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGAc                                                                   >  1:957689/1‑150 (MQ=255)
                                                                                      ttCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGa                                                                <  2:1656359/150‑1 (MQ=255)
                                                                                       tctcCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGaa                                                               >  2:1180903/1‑150 (MQ=255)
                                                                                         tcCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATg                                                             <  2:14977/150‑1 (MQ=255)
                                                                                         tcCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATg                                                             <  2:1789851/150‑1 (MQ=255)
                                                                                          ccTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAAGCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATgc                                                            <  1:524272/150‑1 (MQ=255)
                                                                                           cTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATgcg                                                           <  1:1337305/150‑1 (MQ=255)
                                                                                            tGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATgcgg                                                          >  1:1113565/1‑150 (MQ=255)
                                                                                             gTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATgcggc                                                         >  1:1593432/1‑150 (MQ=255)
                                                                                              ttttGATAGACTTTCCTTATCTTCGTTCTCCGAACAAAACAATTCCAGGCCCGATGGACTGGAAGGCTCAATCGCTCATATCAATCAGAGGATGTTCACATGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATgcggcg                                                        <  2:1151478/150‑1 (MQ=255)
                                                                                               tttGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGa                                                       <  2:1722434/150‑1 (MQ=255)
                                                                                                 tgattgaCTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAACTCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAAt                                                     <  1:369260/150‑1 (MQ=255)
                                                                                                    ttgaCTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGt                                                  <  2:919988/150‑1 (MQ=255)
                                                                                                         ttGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTg                                             <  1:1811091/150‑1 (MQ=255)
                                                                                                          tGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGa                                            <  2:515399/150‑1 (MQ=255)
                                                                                                           gCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGAc                                           >  2:1902517/1‑150 (MQ=255)
                                                                                                           gCCTTATCTTCGTTCTCCGAACAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACt                                          >  1:1277479/1‑150 (MQ=255)
                                                                                                           gCCTTACCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGAc                                           >  1:1267164/1‑150 (MQ=255)
                                                                                                             cTTATCTTCGTTCTACGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACGTGCATATCCGGGTCTTGGCGATCCAGATAGTGGATGGACTGAATGCGGCGCAGGGTGCGTGACtt                                         >  2:271718/1‑150 (MQ=255)
                                                                                                                aTCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGcc                                      >  1:1412380/1‑150 (MQ=255)
                                                                                                                 tCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCcg                                     >  1:837746/1‑150 (MQ=255)
                                                                                                                  cTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCcgc                                    >  2:1695746/1‑150 (MQ=255)
                                                                                                                  cTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCAAGATAGTGGATGGACTGAAAGCGGCGAATGGTGCGTGACTTGCcgc                                    >  2:1855450/1‑150 (MQ=255)
                                                                                                                   ttCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGAAGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCcgcg                                   >  2:235415/1‑150 (MQ=255)
                                                                                                                    tCGTTCTCCGAACATAATAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCgg                                  <  1:1132321/150‑1 (MQ=255)
                                                                                                                        tctcCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATca                              <  2:1593432/150‑1 (MQ=255)
                                                                                                                          tcCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATcagc                            >  1:1894005/1‑150 (MQ=255)
                                                                                                                           ccGAACAAAAAAATTCCAGTACCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGCGGAGGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGTATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATcagca                           <  2:1168786/150‑1 (MQ=255)
                                                                                                                            cGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATcaggag                          >  1:514497/1‑150 (MQ=255)
                                                                                                                             gAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCagcagc                         >  1:207585/1‑150 (MQ=255)
                                                                                                                               aCAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGt                       >  2:1713941/1‑150 (MQ=255)
                                                                                                                               aCAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGt                       <  1:637607/150‑1 (MQ=255)
                                                                                                                                cAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGtt                      >  1:1761835/1‑150 (MQ=255)
                                                                                                                                  aaaaaaaTTCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTTCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGAGACTTGCCGCGGATCAGCAGCGTTTc                    <  2:47485/150‑1 (MQ=255)
                                                                                                                                   aaaaaaTTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCAGTTCCGGGTCTTTGCGATCCAGCTAGTGGATGGACTGAATGCGGCGAATGGTGCGGGACTTGCCGCGGATCAGCAGCGGTTc                    >  1:1226574/1‑150 (MQ=255)
                                                                                                                                    aaaaaTTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCg                   >  2:332006/1‑150 (MQ=255)
                                                                                                                                     aaactaccACTCCCCACGGACTGGAAGGCTCAATCGATCAAATCAATCCGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGGGCGTGACTTGCCGCGGATCAGCAGCGTTTCgg                  <  2:133096/144‑1 (MQ=255)
                                                                                                                                       ccttCCAGTCCTTCGGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTa                >  1:1572485/3‑150 (MQ=255)
                                                                                                                                       ccttCCAGTCCTTCGGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGAGGCGATTGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGtt                >  2:1252649/3‑149 (MQ=255)
                                                                                                                                       ccttCCAGTCCTTCGGGACTGGAAGGCTAAATCGATCAAATCAATCAGAGGATGTGCACCTGGATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATAAGCAGCGTTTCGGTa                >  2:1258892/3‑150 (MQ=255)
                                                                                                                                        cttCCAGTCCTTCGGGAGTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCCCCTCCATTTCCGGGTCTTTGCGATCCAGATAGTGCATGGACTGAATGCCGCGAATGGTGCGTGACATGCCGCGGATCAGCAGCGTTTCGGTAg               >  2:1634491/2‑150 (MQ=255)
                                                                                                                                           ccAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCCCCTGCATTTCCGGGTCTTTGCGCTCCAGATAGGGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCGGCGTTTCGGTAGTcg            >  1:327182/1‑150 (MQ=255)
                                                                                                                                           ccAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTcg            >  2:747919/1‑150 (MQ=255)
                                                                                                                                             aGTCCGGAAGGAGGGGAAGGCTGTATCGAGGAAATCAATCAGAGGATGTGCTCTGGCATTTCCGGGTCTTGGCGTGCCAGCTAGTGGATGGACGGAATGCGGCGACGGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTcgcg          <  1:120305/150‑1 (MQ=255)
                                                                                                                                             aGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTcgcg          >  2:325044/1‑150 (MQ=255)
                                                                                                                                             aGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCAGTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTcgcg          <  1:120302/150‑1 (MQ=255)
                                                                                                                                              gTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGa         >  1:75910/1‑150 (MQ=255)
                                                                                                                                              gTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGa         >  2:894031/1‑150 (MQ=255)
                                                                                                                                                cccGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGata       >  1:1905176/1‑150 (MQ=255)
                                                                                                                                                 ccGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGatat      >  2:991045/1‑150 (MQ=255)
                                                                                                                                                   gAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGATATTg    <  2:818634/150‑1 (MQ=255)
                                                                                                                                                   gAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGGCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGATATTg    <  2:826892/150‑1 (MQ=255)
                                                                                                                                                    aaGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGCCTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGATATTGc   >  2:1683671/1‑150 (MQ=255)
                                                                                                                                                     aGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCCCCTGCATTTCTGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGATATTGcc  >  1:1332108/1‑150 (MQ=255)
                                                                                                                                                     aGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGATATTGcc  <  1:1696143/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 < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: