Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,114,539 A→G 7.6% intergenic (‑67/+30) fpr ← / ← glpX ferredoxin‑NADP reductase; flavodoxin reductase/fructose 1,6‑bisphosphatase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,114,5390AG 7.6% 204.1 / 6.3 66intergenic (‑67/+30)fpr/glpXferredoxin‑NADP reductase; flavodoxin reductase/fructose 1,6‑bisphosphatase II
Reads supporting (aligned to +/- strand):  ref base A (28/33);  new base G (2/3);  total (30/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 = 6.99e-01

GGGGGCGTGAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGAT  >  NC_000913/4114392‑4114682
                                                                                                                                                   |                                                                                                                                               
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                                                                                                                                       <  1:30555/136‑2 (MQ=255)
                      aaCAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCTTCGGGACTGGaat                                                                                                                                       >  2:30555/1‑135 (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                                                                                                                 >  1:61745/1‑128 (MQ=255)
                                                    gTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAatcaatca                                                                                                                 <  2:61745/128‑1 (MQ=255)
                                                     tGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTcc                                                                                           >  2:175245/1‑149 (MQ=255)
                                                         tttGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAaatcaat                                                                                                                   <  3:251829/121‑1 (MQ=255)
                                                         tttGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAaatcaat                                                                                                                   >  4:251829/1‑121 (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                                                                                               <  1:292290/133‑1 (MQ=255)
                                                                 ttACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCAt                                                                                               >  2:292290/1‑133 (MQ=255)
                                                                     ccAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATcc                                                                            <  4:30622/148‑1 (MQ=255)
                                                                     ccAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATcc                                                                            >  3:30622/1‑148 (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                                                                                                               <  1:94062/102‑1 (MQ=255)
                                                                                tGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCaga                                                                                                               >  2:94062/1‑102 (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                                                       <  1:114705/129‑1 (MQ=255)
                                                                                                             ttATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATg                                                       >  2:114705/1‑129 (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)
                                                                                                                                                   |                                                                                                                                               
GGGGGCGTGAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGAT  >  NC_000913/4114392‑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: