Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,114,538 A→C 6.8% intergenic (‑66/+31) 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,5380AC 6.8% 379.0 / 14.1 132intergenic (‑66/+31)fpr/glpXferredoxin‑NADP reductase; flavodoxin reductase/fructose 1,6‑bisphosphatase II
Reads supporting (aligned to +/- strand):  ref base A (65/58);  new base C (4/5);  total (69/63)
Fisher's exact test for biased strand distribution p-value = 7.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.13e-01

Only 100 of 133 total aligned reads displayed.

CGGGGGCGTGAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGATATTGC  >  NC_000913/4114391‑4114687
                                                                                                                                                    |                                                                                                                                                     
cGGGGGCGTGAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAg                                                                                                                                                     <  2:519570/150‑1 (MQ=255)
 gggggCGTGAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAgg                                                                                                                                                    <  2:292432/150‑1 (MQ=255)
   gggCGTGAACGGTGATACTAATCAGGGCGTTGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCTAAC‑AAAAAAAATCCAGTCCCGAAGGAc                                                                                                                                                  <  2:1005033/150‑1 (MQ=255)
   gggCGTGAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGAc                                                                                                                                                  <  1:324187/150‑1 (MQ=255)
      cgggAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTgg                                                                                                                                               <  2:1039369/147‑1 (MQ=255)
       gTGAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCTTCGGGACTGGa                                                                                                                                              <  1:412198/150‑1 (MQ=255)
           aCGGTGAGACTAAACAGGTCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATTTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGc                                                                                                                                          <  2:1407479/150‑1 (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)
                     tAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGatc                                                                                                                                <  2:56773/150‑1 (MQ=255)
                     tAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGatc                                                                                                                                <  2:1457612/150‑1 (MQ=255)
                      aaaCAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGatca                                                                                                                               <  1:87103/150‑1 (MQ=255)
                         cAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGatcaaat                                                                                                                            <  1:1569390/150‑1 (MQ=255)
                         cAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCAAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGatcaaat                                                                                                                            <  2:371011/150‑1 (MQ=255)
                              cGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCATCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAaatcaatc                                                                                                                       <  2:117915/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                                                                                                                 <  1:6873/150‑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)
                                               aCTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGc                                                                                                      <  2:677165/150‑1 (MQ=255)
                                                cTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCACGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGAACCTGCa                                                                                                     >  1:371848/1‑150 (MQ=255)
                                                     gTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGATCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTcc                                                                                                >  1:95011/1‑150 (MQ=255)
                                                      tGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCg                                                                                               <  2:1509380/150‑1 (MQ=255)
                                                       gACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCgg                                                                                              >  1:826994/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: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: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)
                                                                           cAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGt                                                                          >  2:1658570/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)
                                                                              ccATGTTTTTCTCCTGTTTTGTTTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGtgga                                                                       <  1:1476013/150‑1 (MQ=255)
                                                                               cATGTTTTTCTCGTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGTACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGtggat                                                                      <  1:1664742/150‑1 (MQ=255)
                                                                               cATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGtggat                                                                      >  2:1696143/1‑150 (MQ=255)
                                                                                 tGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGtggatgg                                                                    >  1:826892/1‑150 (MQ=255)
                                                                                         tcCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATg                                                            <  2:1789851/150‑1 (MQ=255)
                                                                                         tcCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATg                                                            <  2:14977/150‑1 (MQ=255)
                                                                                          ccTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATgc                                                           >  2:1451366/1‑150 (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)
                                                                                                 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)
                                                                                                         ttGCCTTATATTCGTTCTCCTAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTg                                            <  2:1626454/150‑1 (MQ=255)
                                                                                                           gCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGAc                                          >  2:1902517/1‑150 (MQ=255)
                                                                                                           gCCTTATCTTCGTTCTCCGAAC‑‑AAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACt                                         >  1:1277479/1‑150 (MQ=255)
                                                                                                                aTCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGcc                                     >  1:1412380/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)
                                                                                                                         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                  >  2:332006/1‑150 (MQ=255)
                                                                                                                                      aaaTTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTACGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGt                <  2:1269439/150‑1 (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)
                                                                                                                                             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)
                                                                                                                                                cccGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGata      >  1:1905176/1‑150 (MQ=255)
                                                                                                                                                 ccGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGatat     >  2:991045/1‑150 (MQ=255)
                                                                                                                                                 ccGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGatat     >  2:1700908/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)
                                                                                                                                                    |                                                                                                                                                     
CGGGGGCGTGAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAAC‑AAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGATATTGC  >  NC_000913/4114391‑4114687

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