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.

CGGGGGCGTGAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGATATTGCC  >  NC_000913/4114391‑4114688
                                                                                                                                                     |                                                                                                                                                     
   gggCGTGAACGGTGATACTAATCAGGGCGTTGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCTAAC‑AAAAAAAATCCAGTCCCGAAGGAc                                                                                                                                                   <  2:1005033/150‑1 (MQ=255)
   gggCGTGAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGAc                                                                                                                                                   <  1:324187/150‑1 (MQ=255)
       gTGAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCTTCGGGACTGGa                                                                                                                                               <  1:412198/150‑1 (MQ=255)
          aaCGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACACAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCACCGAAC‑AAAAAAATTCCAGTCCTTCGGGACTGGaatt                                                                                                                                            <  2:57362/150‑3 (MQ=255)
           aCGGTGAGACTAAACAGGTCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATTTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGc                                                                                                                                           <  2:1407479/150‑1 (MQ=255)
            cGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCTTCGGGACTGGaatttt                                                                                                                                          <  1:1291195/150‑5 (MQ=255)
            cGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCTTCGGGACTGGaatttt                                                                                                                                          <  1:864536/150‑5 (MQ=255)
             ggTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTc                                                                                                                                         >  1:1724369/1‑150 (MQ=255)
              gTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCTTCGGGACTGGaatttttt                                                                                                                                        <  1:979348/150‑7 (MQ=255)
                  gaCTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCg                                                                                                                                    >  2:706589/1‑150 (MQ=255)
                  gaCTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCg                                                                                                                                    >  2:367630/1‑150 (MQ=255)
                     tAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGatc                                                                                                                                 <  2:56773/150‑1 (MQ=255)
                     tAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGatc                                                                                                                                 <  2:1457612/150‑1 (MQ=255)
                     tAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGatc                                                                                                                                 <  2:1099166/150‑1 (MQ=255)
                      aaaCAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGatca                                                                                                                                <  1:87103/150‑1 (MQ=255)
                         cAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGatcaaat                                                                                                                             <  1:1569390/150‑1 (MQ=255)
                           gggCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGatcaaatca                                                                                                                           <  2:1208612/150‑1 (MQ=255)
                                  ggTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCagag                                                                                                                    >  2:165664/1‑150 (MQ=255)
                                   gTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAgg                                                                                                                   >  1:1269439/1‑150 (MQ=255)
                                   gTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAgg                                                                                                                   >  2:187956/1‑150 (MQ=255)
                                    tCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGa                                                                                                                  <  1:6873/150‑1 (MQ=255)
                                    tCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGa                                                                                                                  <  2:529041/150‑1 (MQ=255)
                                      cAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGAtg                                                                                                                >  2:156190/1‑150 (MQ=255)
                                               aCTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGc                                                                                                       <  2:677165/150‑1 (MQ=255)
                                                cTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCACGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGAACCTGCa                                                                                                      >  1:371848/1‑150 (MQ=255)
                                                              ccTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTg                                                                                        >  2:120302/1‑150 (MQ=255)
                                                              ccTGTTACCCAATCACGCATGTTGATCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTGTCCTAAC‑AAAACAATTCGAGCCCCGAAGGACTGGAAGGCTCAATCGATCACATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTg                                                                                        <  1:940593/150‑1 (MQ=255)
                                                               cTGTTTCACATTCAGCCATGTTTTTCTCCTGTTTTGATTTACTTGCCTTATCTTCGTTCTCCTAAC‑AAAAAAAATCCAGTCACTAAGGACTGGAAGGCTCAATCGTTCAAATCAATCAGAGGATGTGCTCCTGCATTTCCGGGTCTTTGc                                                                                       <  2:265781/150‑1 (MQ=255)
                                                               cTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTATGc                                                                                       >  2:226710/1‑150 (MQ=255)
                                                                 gTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCa                                                                                                      <  1:907943/133‑1 (MQ=255)
                                                                   tACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATc                                                                                   >  1:1060434/1‑150 (MQ=255)
                                                                   tACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGACCCTGCATTTCCGGGTCTTTGAGATc                                                                                   >  2:260301/1‑150 (MQ=255)
                                                                    aCCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATcc                                                                                  <  1:466281/150‑1 (MQ=255)
                                                                     cccAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCa                                                                                 >  1:1454400/1‑150 (MQ=255)
                                                                     cccAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGAATGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCa                                                                                 <  1:270102/150‑1 (MQ=255)
                                                                       cAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAg                                                                                <  2:1604334/150‑1 (MQ=255)
                                                                       cAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGa                                                                               >  1:919988/1‑150 (MQ=255)
                                                                             gCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGtgg                                                                         >  2:1480450/1‑150 (MQ=255)
                                                                             gCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGtgg                                                                         >  2:477527/1‑150 (MQ=255)
                                                                             gCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCACGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGtgg                                                                         >  1:59960/1‑150 (MQ=255)
                                                                             gCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTAGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGtgg                                                                         <  2:1194642/150‑1 (MQ=255)
                                                                               cATGTTTTTCTCGTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGTACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGtggat                                                                       <  1:1664742/150‑1 (MQ=255)
                                                                               cATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGtggat                                                                       >  2:1696143/1‑150 (MQ=255)
                                                                                aTGTTTTTCTCCTGTTTTGTTTGACTTGCCTTATCTTCGTTCTACGAAC‑AAAAAAATTCCAGTCCCGACGGACTGGAAGGCTCATTCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGtggatg                                                                      <  1:1449130/150‑1 (MQ=255)
                                                                                 tGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGtggatgg                                                                     >  1:826892/1‑150 (MQ=255)
                                                                                   tttttCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGAc                                                                   >  1:957689/1‑150 (MQ=255)
                                                                                      ttCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGa                                                                <  2:1656359/150‑1 (MQ=255)
                                                                                       tctcCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGaa                                                               >  2:1180903/1‑150 (MQ=255)
                                                                                         tcCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATg                                                             <  2:14977/150‑1 (MQ=255)
                                                                                         tcCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATg                                                             <  2:1789851/150‑1 (MQ=255)
                                                                                          ccTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAAGCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATgc                                                            <  1:524272/150‑1 (MQ=255)
                                                                                           cTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATgcg                                                           <  1:1337305/150‑1 (MQ=255)
                                                                                            tGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATgcgg                                                          >  1:1113565/1‑150 (MQ=255)
                                                                                             gTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATgcggc                                                         >  1:1593432/1‑150 (MQ=255)
                                                                                              ttttGATAGACTTTCCTTATCTTCGTTCTCCGAAC‑AAAACAATTCCAGGCCCGATGGACTGGAAGGCTCAATCGCTCATATCAATCAGAGGATGTTCACATGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATgcggcg                                                        <  2:1151478/150‑1 (MQ=255)
                                                                                               tttGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGa                                                       <  2:1722434/150‑1 (MQ=255)
                                                                                                 tgattgaCTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAACTCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAAt                                                     <  1:369260/150‑1 (MQ=255)
                                                                                                    ttgaCTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGt                                                  <  2:919988/150‑1 (MQ=255)
                                                                                                         ttGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTg                                             <  1:1811091/150‑1 (MQ=255)
                                                                                                          tGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGa                                            <  2:515399/150‑1 (MQ=255)
                                                                                                           gCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGAc                                           >  2:1902517/1‑150 (MQ=255)
                                                                                                           gCCTTATCTTCGTTCTCCGAAC‑‑AAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACt                                          >  1:1277479/1‑150 (MQ=255)
                                                                                                           gCCTTACCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGAc                                           >  1:1267164/1‑150 (MQ=255)
                                                                                                             cTTATCTTCGTTCTACGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACGTGCATATCCGGGTCTTGGCGATCCAGATAGTGGATGGACTGAATGCGGCGCAGGGTGCGTGACtt                                         >  2:271718/1‑150 (MQ=255)
                                                                                                                aTCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGcc                                      >  1:1412380/1‑150 (MQ=255)
                                                                                                                 tCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCcg                                     >  1:837746/1‑150 (MQ=255)
                                                                                                                  cTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCcgc                                    >  2:1695746/1‑150 (MQ=255)
                                                                                                                  cTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCAAGATAGTGGATGGACTGAAAGCGGCGAATGGTGCGTGACTTGCcgc                                    >  2:1855450/1‑150 (MQ=255)
                                                                                                                   ttCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGAAGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCcgcg                                   >  2:235415/1‑150 (MQ=255)
                                                                                                                    tCGTTCTCCGAAC‑ATAATAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCgg                                  <  1:1132321/150‑1 (MQ=255)
                                                                                                                        tctcCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATca                              <  2:1593432/150‑1 (MQ=255)
                                                                                                                          tcCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATcagc                            >  1:1894005/1‑150 (MQ=255)
                                                                                                                           ccGAAC‑AAAAAAATTCCAGTACCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGCGGAGGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGTATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATcagca                           <  2:1168786/150‑1 (MQ=255)
                                                                                                                            cGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATcaggag                          >  1:514497/1‑150 (MQ=255)
                                                                                                                             gAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCagcagc                         >  1:207585/1‑150 (MQ=255)
                                                                                                                               aC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGt                       >  2:1713941/1‑150 (MQ=255)
                                                                                                                               aC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGt                       <  1:637607/150‑1 (MQ=255)
                                                                                                                                c‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGtt                      >  1:1761835/1‑150 (MQ=255)
                                                                                                                                  aaaaaaaTT‑CAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTTCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGAGACTTGCCGCGGATCAGCAGCGTTTc                    <  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)
                                                                                                                                                     |                                                                                                                                                     
CGGGGGCGTGAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGATATTGCC  >  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: ‑