Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,114,536 C→T 7.2% intergenic (‑64/+33) 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,5360CT 7.2% 383.0 / 12.3 126intergenic (‑64/+33)fpr/glpXferredoxin‑NADP reductase; flavodoxin reductase/fructose 1,6‑bisphosphatase II
Reads supporting (aligned to +/- strand):  ref base C (61/55);  new base T (4/5);  total (65/61)
Fisher's exact test for biased strand distribution p-value = 7.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 130 total aligned reads displayed.

CGGGGGCGTGAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGATAT  >  NC_000913/4114391‑4114684
                                                                                                                                                  |                                                                                                                                                    
   gggCGTGAACGGTGATACTAATCAGGGCGTTGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCTAAC‑AAAAAAAATCCAGTCCCGAAGGAc                                                                                                                                               <  2:1005033/150‑1 (MQ=255)
      cgggAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTgg                                                                                                                                            <  2:1039369/147‑1 (MQ=255)
       gTGAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCTTCGGGACTGGa                                                                                                                                           <  1:412198/150‑1 (MQ=255)
          aaCGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACACAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCACCGAAC‑AAAAAAATTCCAGTCCTTCGGGACTGGaatt                                                                                                                                        <  2:57362/150‑3 (MQ=255)
            cGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCTTCGGGACTGGaatttt                                                                                                                                      <  1:864536/150‑5 (MQ=255)
              gTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCTTCGGGACTGGaatttttt                                                                                                                                    <  1:979348/150‑7 (MQ=255)
                  gaCTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCg                                                                                                                                >  2:367630/1‑150 (MQ=255)
                  gaCTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCg                                                                                                                                >  2:706589/1‑150 (MQ=255)
                     tAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGatc                                                                                                                             <  2:1457612/150‑1 (MQ=255)
                     tAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGatc                                                                                                                             <  2:56773/150‑1 (MQ=255)
                     tAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGatc                                                                                                                             <  1:757377/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)
                             gCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAaatcaat                                                                                                                     >  1:1830029/1‑150 (MQ=255)
                              cGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCATCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAaatcaatc                                                                                                                    <  2:117915/150‑1 (MQ=255)
                                  ggTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTGAATCGATCAAATCAATCagag                                                                                                                >  1:1724059/1‑150 (MQ=255)
                                  ggTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCagag                                                                                                                <  1:423645/150‑1 (MQ=255)
                                  ggTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCagag                                                                                                                >  2:165664/1‑150 (MQ=255)
                                  ggTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCagag                                                                                                                >  1:1261998/1‑150 (MQ=255)
                                   gTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAgg                                                                                                               >  1:1269439/1‑150 (MQ=255)
                                    tacAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCTGAGGa                                                                                                              <  1:1287078/148‑1 (MQ=255)
                                    tCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGa                                                                                                              <  2:529041/150‑1 (MQ=255)
                                      cAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGAtg                                                                                                            >  2:156190/1‑150 (MQ=255)
                                            tGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCAcc                                                                                                      >  1:883975/1‑150 (MQ=255)
                                                     gTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGATCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTcc                                                                                             >  1:95011/1‑150 (MQ=255)
                                                            tGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCtt                                                                                      >  2:1811091/1‑150 (MQ=255)
                                                              ccTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTg                                                                                    >  2:120302/1‑150 (MQ=255)
                                                              ccTGTTACCCAATCACGCATGTTGATCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTGTCCTAAC‑AAAACAATTCGAGCCCCGAAGGACTGGAAGGCTCAATCGATCACATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTg                                                                                    <  1:940593/150‑1 (MQ=255)
                                                              acTGTTACCCAATCACCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTAAGAAC‑AAAATAATTCCAGTCCCCAAGGACTGTAAGGCTCAATCTATCAACTCAAGCAGAGGCTGGGCACCTGCAGTTCCAGGTCttag                                                                                    >  2:120305/2‑148 (MQ=255)
                                                               cTGTTTCACATTCAGCCATGTTTTTCTCCTGTTTTGATTTACTTGCCTTATCTTCGTTCTCCTAAC‑AAAAAAAATCCAGTCACTAAGGACTGGAAGGCTCAATCGTTCAAATCAATCAGAGGATGTGCTCCTGCATTTCCGGGTCTTTGc                                                                                   <  2:265781/150‑1 (MQ=255)
                                                               cTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTATGc                                                                                   >  2:226710/1‑150 (MQ=255)
                                                                 gTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCa                                                                                                  >  2:907943/1‑133 (MQ=255)
                                                                   tACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATc                                                                               >  1:248436/1‑150 (MQ=255)
                                                                   tACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATc                                                                               >  1:1194642/1‑150 (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:367630/150‑1 (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)
                                                                              ccATGTTTTTCTCCTGTTTTGTTTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGtgga                                                                    <  1:1476013/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)
                                                                                      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‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATgc                                                        >  2:1451366/1‑150 (MQ=255)
                                                                                           cTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATgcg                                                       <  1:1337305/150‑1 (MQ=255)
                                                                                            tGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATgcgg                                                      >  1:1113565/1‑150 (MQ=255)
                                                                                              ttttGATAGACTTTCCTTATCTTCGTTCTCCGAAC‑AAAACAATTCCAGGCCCGATGGACTGGAAGGCTCAATCGCTCATATCAATCAGAGGATGTTCACATGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATgcggcg                                                    <  2:1151478/150‑1 (MQ=255)
                                                                                                 tgattgaCTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAACTCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAAt                                                 <  1:369260/150‑1 (MQ=255)
                                                                                                    ttgaCTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGt                                              <  2:919988/150‑1 (MQ=255)
                                                                                                      gaCTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGc                                            >  2:215065/1‑150 (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)
                                                                                                             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)
                                                                                                                    tCGTTCTCCGAAC‑ATAATAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCgg                              <  1:1132321/150‑1 (MQ=255)
                                                                                                                       ttCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATc                           <  2:248436/150‑1 (MQ=255)
                                                                                                                        tctcCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATca                          <  2:1593432/150‑1 (MQ=255)
                                                                                                                         ctcCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATcag                         >  2:811384/1‑150 (MQ=255)
                                                                                                                         ctcCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGGCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATcag                         <  1:165664/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:1405812/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               >  1:587831/1‑150 (MQ=255)
                                                                                                                                    aaaaaTTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCg               >  2:332006/1‑150 (MQ=255)
                                                                                                                                     aaactaccACTCCCCACGGACTGGAAGGCTCAATCGATCAAATCAATCCGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGGGCGTGACTTGCCGCGGATCAGCAGCGTTTCgg              <  2:133096/144‑1 (MQ=255)
                                                                                                                                      aaaTTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTACGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGt             <  2:1269439/150‑1 (MQ=255)
                                                                                                                                       ccttCCAGTCCTTCGGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTa            >  1:1572485/3‑150 (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)
                                                                                                                                              gTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGa     >  1:75910/1‑150 (MQ=255)
                                                                                                                                                cccGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGata   >  2:64004/1‑150 (MQ=255)
                                                                                                                                                 ccGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGatat  >  2:1700908/1‑150 (MQ=255)
                                                                                                                                                 ccGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGatat  >  2:991045/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CGGGGGCGTGAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGATAT  >  NC_000913/4114391‑4114684

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: ‑