Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 378,400 G→A 5.4% intergenic (‑32/+62) frmB ← / ← frmA S‑formylglutathione hydrolase/alcohol dehydrogenase class III; glutathione‑dependent formaldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913378,4000GA 5.4% 298.4 / 4.2 93intergenic (‑32/+62)frmB/frmAS‑formylglutathione hydrolase/alcohol dehydrogenase class III; glutathione‑dependent formaldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (45/42);  new base A (2/3);  total (47/46)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAGGCTCATGG  >  NC_000913/378253‑378541
                                                                                                                                                   |                                                                                                                                             
gcgGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACgg                                                                                                                                              >  2:321403/1‑149 (MQ=255)
          gtggGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCggg                                                                                                                                     <  1:321403/148‑1 (MQ=255)
          gtggGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCggg                                                                                                                                     <  3:231265/148‑1 (MQ=255)
           tggGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGGCAGCCGTCCgggg                                                                                                                                    <  3:359693/148‑1 (MQ=255)
           tggGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGga                                                                                                                                   <  1:233975/149‑1 (MQ=255)
           tggGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGga                                                                                                                                   <  1:290047/149‑1 (MQ=255)
             ggAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGGCAGCCGTCCggggga                                                                                                                                  <  6:246986/148‑3 (MQ=255)
             ggAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGgaga                                                                                                                                 <  5:225323/149‑1 (MQ=255)
             ggAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGga                                                                                                                                   >  3:275366/1‑147 (MQ=255)
              gagaTAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGagag                                                                                                                                <  3:121160/149‑1 (MQ=255)
               agaTAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGagag                                                                                                                                >  5:242572/1‑148 (MQ=255)
                gaTAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCggg                                                                                                                                     <  3:346420/142‑1 (MQ=255)
                gaTAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCggg                                                                                                                                     >  4:346420/1‑142 (MQ=255)
                    gACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTAttt                                                                                                                           <  1:343537/148‑2 (MQ=255)
                    gACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTaa                                                                                                                          <  6:242572/149‑1 (MQ=255)
                    gACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTaa                                                                                                                          <  4:387425/149‑1 (MQ=255)
                      cGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAg                                                                                                                         >  5:372272/1‑148 (MQ=255)
                          gACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAAccc                                                                                                                    >  8:309237/1‑149 (MQ=255)
                          gACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCAGGACGGCTGCCCGGGGAGAGTATTAAGAAccc                                                                                                                    >  4:269364/1‑149 (MQ=255)
                           aCATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAAcccc                                                                                                                   >  5:208134/1‑149 (MQ=255)
                             attcattTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCg                                                                                                                 >  8:37770/1‑149 (MQ=255)
                              ttcattTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGc                                                                                                                <  4:275366/149‑1 (MQ=255)
                              ttcattTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGc                                                                                                                <  5:290176/149‑1 (MQ=255)
                                  tttcacatttcaGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTg                                                                                                             >  8:67813/1‑148 (MQ=255)
                                   ttcacatttcaGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTg                                                                                                             >  8:109202/1‑147 (MQ=255)
                                     cacatttcaGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGt                                                                                                         >  2:113505/1‑149 (MQ=255)
                                         attcaGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGAc                                                                                                     >  8:397092/2‑149 (MQ=255)
                                            caGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACggg                                                                                                  >  6:176606/1‑149 (MQ=255)
                                             aGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACgggg                                                                                                 >  3:377244/1‑149 (MQ=255)
                                             aGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACggg                                                                                                  <  1:119842/148‑1 (MQ=255)
                                                 aTTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTaaa                                                                                             >  7:123862/1‑149 (MQ=255)
                                                  ttGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAAc                                                                                            <  7:309237/149‑1 (MQ=255)
                                                  ttGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATATCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGCCGGGGTAAAc                                                                                            >  5:239623/1‑149 (MQ=255)
                                                       aataatGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCg                                                                                       <  7:122653/149‑1 (MQ=255)
                                                       aataatGCAGATACATATTTTGACAGCCGCCAACGATGACATGTTTTTCAATGAGTTCCATCAGTCATTATAGCAGGGTTAATCCCCCGGACTGCTGCCCGGGGAGATTATTAAGACCCCCCGCCTGAAGTGGACGGGGTAAACCTGCg                                                                                       <  6:22538/149‑1 (MQ=255)
                                                         taatGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCggg                                                                                     <  2:421226/149‑1 (MQ=255)
                                                           atGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGga                                                                                                                                   >  4:244754/1‑101 (MQ=255)
                                                           atGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGga                                                                                                                                   <  3:244754/101‑1 (MQ=255)
                                                           atGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAAccccc                                                                                                                  >  5:325189/1‑118 (MQ=255)
                                                           atGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAAccccc                                                                                                                  <  6:325189/118‑1 (MQ=255)
                                                              ccGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAAt                                                                                 >  3:369414/1‑148 (MQ=255)
                                                                 aTACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGt                                                                             <  6:285222/149‑1 (MQ=255)
                                                                 aTACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGt                                                                             <  6:208134/149‑1 (MQ=255)
                                                                     atatTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAAc                                                                          >  3:334916/1‑148 (MQ=255)
                                                                        ttttGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGaa                                                                       <  7:37770/148‑1 (MQ=255)
                                                                              cAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGt                                                                >  7:8732/1‑149 (MQ=255)
                                                                               aGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGtt                                                               >  3:375418/1‑149 (MQ=255)
                                                                               aGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGtt                                                               >  5:239330/1‑149 (MQ=255)
                                                                                gccgccAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTc                                                              <  1:113505/149‑1 (MQ=255)
                                                                                gccgccAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTc                                                              >  2:263274/1‑149 (MQ=255)
                                                                                gccgccAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTc                                                              <  3:269364/149‑1 (MQ=255)
                                                                                gccgccAAAGATGACATTTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCTCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACTGTTc                                                              <  1:223746/149‑1 (MQ=255)
                                                                                    ccAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAAt                                                          >  1:384749/1‑149 (MQ=255)
                                                                                       aaGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGa                                                       >  3:175210/1‑149 (MQ=255)
                                                                                          cTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGAttt                                                    <  6:372272/149‑1 (MQ=255)
                                                                                           tGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTg                                                   >  5:385976/1‑149 (MQ=255)
                                                                                               aTGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGTTTCGAATGGAttt                                                    >  6:279767/1‑144 (MQ=255)
                                                                                               aTGTTTTTCAATGAGTTCCATCAGTAATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGTTTCGAATGGAttt                                                    <  5:279767/144‑1 (MQ=255)
                                                                                                tGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGa                                                       >  7:342775/1‑140 (MQ=255)
                                                                                                tGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGa                                                       <  8:342775/140‑1 (MQ=255)
                                                                                                tGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTc                                              >  3:104650/1‑149 (MQ=255)
                                                                                                 gTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCa                                             >  1:52531/1‑149 (MQ=255)
                                                                                                  tttttCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCAt                                            <  4:369414/149‑1 (MQ=255)
                                                                                                  tttttCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCAt                                            >  6:89516/1‑149 (MQ=255)
                                                                                                   ttttCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCAt                                            <  7:67813/148‑1 (MQ=255)
                                                                                                        aaTGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATc                                       <  4:334916/148‑1 (MQ=255)
                                                                                                         aTGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAg                                     >  7:328585/1‑149 (MQ=255)
                                                                                                          tGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTc                                              <  6:369177/139‑1 (MQ=255)
                                                                                                          tGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTc                                              >  5:369177/1‑139 (MQ=255)
                                                                                                           gAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATc                                       >  7:54000/1‑145 (MQ=255)
                                                                                                           gAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATc                                       <  8:54000/145‑1 (MQ=255)
                                                                                                               tCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCg                                 >  8:20097/1‑147 (MQ=255)
                                                                                                                ccATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGaa                               <  5:89516/148‑1 (MQ=255)
                                                                                                                 cATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAgg                             >  8:41619/1‑149 (MQ=255)
                                                                                                                       tCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCAt                       >  5:302639/1‑149 (MQ=255)
                                                                                                                         aTTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCAtt                      >  1:20153/1‑148 (MQ=255)
                                                                                                                           tATCTCAGGGTTACTCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCGCCGCCTGTAGTGGACGGGGTAATCCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGACGGCGTCATTaa                    <  6:13934/148‑1 (MQ=255)
                                                                                                                           tATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGTGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATg                                           <  8:34509/125‑1 (MQ=255)
                                                                                                                           tATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGTGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATg                                           >  7:34509/1‑125 (MQ=255)
                                                                                                                            aTCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACACCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAAtt                  >  1:27092/1‑149 (MQ=255)
                                                                                                                              ctcAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCCGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGaa                               >  7:293528/1‑134 (MQ=255)
                                                                                                                              ctcAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCCGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGaa                               <  8:293528/134‑1 (MQ=255)
                                                                                                                                 aGGGTTAATCCCCCGGACGGCGGCCCGGGGAGAGTATTAAGAACCCCCGCCTGCAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAAttt                 <  1:338518/145‑1 (MQ=255)
                                                                                                                                   ggTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCa           <  8:123862/149‑1 (MQ=255)
                                                                                                                                   ggTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCa           <  1:263274/149‑1 (MQ=255)
                                                                                                                                      tAATCCCCCGGACGGCTGCCCGGGGAGAGTTTTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAGGc        >  7:105886/1‑149 (MQ=255)
                                                                                                                                        aTCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAGGCt       <  7:41619/148‑1 (MQ=255)
                                                                                                                                         tCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTACCGAATTACGGCTCGAATGGATTTGCCTGCATGCATCAGGTCGAAGGCGGCATTAATTTCATCCAGGCTCa     <  4:173685/149‑1 (MQ=255)
                                                                                                                                          tccccGGGCAGCCGTCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGCAGGCGTCATTAATTTCATCCAGGCTCAt    >  2:329069/2‑149 (MQ=255)
                                                                                                                                           ccccGGGCAGCCGTCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAgg                                    >  4:105904/1‑116 (MQ=255)
                                                                                                                                           ccccGGGCAGCCGTCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAgg                                    <  3:105904/116‑1 (MQ=255)
                                                                                                                                           ccccGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAGGCTCAt    <  2:52531/148‑1 (MQ=255)
                                                                                                                                           ccccGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAGGCTCATg   >  4:27962/1‑149 (MQ=255)
                                                                                                                                             ccGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAGGCTCATgg  >  1:426285/1‑148 (MQ=255)
                                                                                                                                                   ggCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAAtt                  >  5:14346/1‑126 (MQ=255)
                                                                                                                                                   ggCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAAtt                  <  6:14346/126‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
GCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAGGCTCATGG  >  NC_000913/378253‑378541

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