Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,114,537 G→T 7.5% intergenic (‑65/+32) 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,5370GT 7.5% 229.6 / 7.4 67intergenic (‑65/+32)fpr/glpXferredoxin‑NADP reductase; flavodoxin reductase/fructose 1,6‑bisphosphatase II
Reads supporting (aligned to +/- strand):  ref base G (29/33);  new base T (2/3);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.86e-01

ACGGGGGCGTGAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGAT  >  NC_000913/4114390‑4114682
                                                                                                                                                   |                                                                                                                                                 
aCGGGGGCGTGAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGa                                                                                                                                                  >  2:353516/1‑149 (MQ=255)
  gggggCGTGAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAg                                                                                                                                                <  3:140063/149‑1 (MQ=255)
   ggggCGTGAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAgg                                                                                                                                               <  4:241985/149‑1 (MQ=255)
          gAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCTTCGGGACTGGaa                                                                                                                                        <  1:331329/149‑1 (MQ=255)
           aaCGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAg                                                                                                                                       <  1:137841/149‑1 (MQ=255)
              ggTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCt                                                                                                                                    <  3:400467/149‑1 (MQ=255)
                      tAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGat                                                                                                                            <  2:99462/149‑1 (MQ=255)
                        aaCAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCTTCGGGACTGGaat                                                                                                                                       >  2:30555/1‑135 (MQ=255)
                        aaCAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCTTCGGGACTGGaat                                                                                                                                       <  1:30555/136‑2 (MQ=255)
                        aaCAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGatca                                                                                                                          >  2:58721/1‑149 (MQ=255)
                        aaCAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGat                                                                                                                            <  1:293999/147‑1 (MQ=255)
                        aaCAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGat                                                                                                                            >  2:293999/1‑147 (MQ=255)
                            gggCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGatcaaatc                                                                                                                      <  4:444695/149‑1 (MQ=255)
                                          ttCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGAtgtg                                                                                                        <  4:282081/149‑1 (MQ=255)
                                          ttCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGAtgtg                                                                                                        <  1:52911/149‑1 (MQ=255)
                                                aCTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTg                                                                                                  >  1:210480/1‑149 (MQ=255)
                                                aCTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTg                                                                                                  >  2:263455/1‑149 (MQ=255)
                                                      gTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAatcaatca                                                                                                                 <  2:61745/128‑1 (MQ=255)
                                                      gTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAatcaatca                                                                                                                 >  1:61745/1‑128 (MQ=255)
                                                       tGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTcc                                                                                           >  2:175245/1‑149 (MQ=255)
                                                           tttGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAaatcaat                                                                                                                   >  4:251829/1‑121 (MQ=255)
                                                           tttGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAaatcaat                                                                                                                   <  3:251829/121‑1 (MQ=255)
                                                           tttGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGt                                                                                       <  1:263455/149‑1 (MQ=255)
                                                           tttGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTAATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGt                                                                                       >  1:402908/1‑149 (MQ=255)
                                                           tttGCCTGTTACCCAATCAGCCATGTTTGTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGt                                                                                       >  3:31822/1‑149 (MQ=255)
                                                                 tGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGc                                                                                 <  1:68908/149‑1 (MQ=255)
                                                                  gTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCgg                                                                                         <  2:162565/140‑1 (MQ=255)
                                                                  gTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCgg                                                                                         >  1:162565/1‑140 (MQ=255)
                                                                   ttACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCAt                                                                                               >  2:292290/1‑133 (MQ=255)
                                                                   ttACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCAt                                                                                               <  1:292290/133‑1 (MQ=255)
                                                                       ccAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATcc                                                                            >  3:30622/1‑148 (MQ=255)
                                                                       ccAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATcc                                                                            <  4:30622/148‑1 (MQ=255)
                                                                         aaTCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAatcaatca                                                                                                                 <  3:128157/109‑1 (MQ=255)
                                                                         aaTCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAatcaatca                                                                                                                 >  4:128157/1‑109 (MQ=255)
                                                                         aaTCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGa                                                                         <  3:120686/149‑1 (MQ=255)
                                                                              gCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGtg                                                                    <  1:175245/149‑1 (MQ=255)
                                                                                  tGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCaga                                                                                                               >  2:94062/1‑102 (MQ=255)
                                                                                  tGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCaga                                                                                                               <  1:94062/102‑1 (MQ=255)
                                                                                  tGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGtg                                                                    >  4:337549/1‑145 (MQ=255)
                                                                                  tGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGtg                                                                    <  3:337549/145‑1 (MQ=255)
                                                                                     ttttCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGAc                                                             <  2:210480/149‑1 (MQ=255)
                                                                                              gTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATgcgg                                                    >  4:271909/1‑149 (MQ=255)
                                                                                               ttttGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATgcggc                                                   >  2:358206/1‑149 (MQ=255)
                                                                                                    attgaCCTGACTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCCGCATTTCCGGGTCTTTGCGATCCAGCTAGTGGATGGACTGAATGCGGCGAATg                                              <  2:402908/149‑1 (MQ=255)
                                                                                                      tgaCTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGt                                            <  1:58721/149‑1 (MQ=255)
                                                                                                          ttGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGt                                        >  2:159230/1‑149 (MQ=255)
                                                                                                               ttATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATg                                                       >  2:114705/1‑129 (MQ=255)
                                                                                                               ttATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATg                                                       <  1:114705/129‑1 (MQ=255)
                                                                                                               ttATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACtt                                   >  4:439870/1‑149 (MQ=255)
                                                                                                                   cTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAAt                                                        <  3:210792/124‑1 (MQ=255)
                                                                                                                   cTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAAt                                                        >  4:210792/1‑124 (MQ=255)
                                                                                                                       gTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGa                           >  1:377779/1‑149 (MQ=255)
                                                                                                                       gTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGa                           <  1:202999/149‑1 (MQ=255)
                                                                                                                          ctcCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATca                        <  3:439870/149‑1 (MQ=255)
                                                                                                                           tcCGAACAAAAAAATGCCAGTCCCGAAGGAGTGGAGGGCTCAATCGATCAAATCAATCAGAGGATGTGCGCCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGTATGGTGCGTGACTTGCCGCGGATcag                       >  4:295113/1‑149 (MQ=255)
                                                                                                                            ccGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATcagc                      >  2:201069/1‑149 (MQ=255)
                                                                                                                                aCAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCg                  >  4:70068/1‑149 (MQ=255)
                                                                                                                                  aaaaaaaTTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCgg                            >  2:401385/1‑137 (MQ=255)
                                                                                                                                  aaaaaaaTTCCAGTCCCGAAGGACTGGAAGGCTCAAACGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCgg                            <  1:401385/137‑1 (MQ=255)
                                                                                                                                     aaaaTTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCGGCAGCGTTTCg             <  1:353516/149‑1 (MQ=255)
                                                                                                                                        cttCCAGTCCTTCGGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGtgga                                                                  <  4:322774/92‑1 (MQ=255)
                                                                                                                                         ttCCAGTCCTTCGGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAg         >  4:144671/1‑149 (MQ=255)
                                                                                                                                            cAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGGGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTcg      >  4:172462/1‑149 (MQ=255)
                                                                                                                                            cAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGa                                      >  4:239744/1‑117 (MQ=255)
                                                                                                                                            cAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGa                                      <  3:239744/117‑1 (MQ=255)
                                                                                                                                               tCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGa   <  2:377779/149‑1 (MQ=255)
                                                                                                                                                cccGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGat  <  1:358206/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
ACGGGGGCGTGAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGAT  >  NC_000913/4114390‑4114682

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