Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,683,472 G→A 5.0% intergenic (‑25/+158) dctA ← / ← yhjK C4‑dicarboxylic acid, orotate and citrate transporter/cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,683,4720GA 5.0% 387.6 / 5.5 121intergenic (‑25/+158)dctA/yhjKC4‑dicarboxylic acid, orotate and citrate transporter/cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (66/46);  new base A (3/3);  total (70/51)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 124 total aligned reads displayed.

GAAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAG  >  NC_000913/3683325‑3683619
                                                                                                                                                   |                                                                                                                                                   
ccgtcgccaagcggtttcatttgctcgcctatttcaggatagaaatggccaaggagaataccaatggctatcgctgtcaggacctgaaagtaaaggcttttaaacagagaggttttcatagggtgtcctttagtaaaaccacagacag                                                                                                                                                     >  8:483062‑M2/145‑148 (MQ=255)
  aGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGt                                                                                                                                                   <  6:421055/148‑1 (MQ=255)
        cGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGt                                                                                                                                              >  8:399203/1‑147 (MQ=255)
        cGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGt                                                                                                                                              >  8:171145/1‑147 (MQ=255)
        cGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGt                                                                                                                                              >  4:85049/1‑147 (MQ=255)
         gCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGCTGTTTTCATAGGGTGTCCTTTAGTAAAACCACATATCTTGTa                                                                                                                                             <  1:122137/147‑1 (MQ=255)
         gCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTaa                                                                                                                                            <  5:297115/148‑1 (MQ=255)
         gCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGATCTTGTa                                                                                                                                             <  7:171469/147‑1 (MQ=255)
             agccgcttCCTTTGCCCGCCTATTTCACGACACAAATCCCCAAGGAGAATACCAATGGCTATCGCGGTCAGGACCTGAAAGTAAACGCTCTTAAACAGAGAGGGTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGtt                                                                                                                                        <  4:355910/142‑1 (MQ=255)
             aGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGt                                                                                                                                         >  2:276414/1‑147 (MQ=255)
               gggTTTCTTTTGGTCGCGTTTTTCTGGATTGTATTGGGCTAGGAGTATTCGTTTGGCTATCGCTGTCAGGACCTGAAAGTTAAGGCTTTTAAACAGTGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTa                                                                                                                                       <  4:432167/146‑1 (MQ=255)
                 gTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGATCTTGTAAGGTTATgg                                                                                                                                    <  1:350993/148‑1 (MQ=255)
                     cATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTaa                                                                                                                                            <  4:130952/136‑1 (MQ=255)
                     cATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTaa                                                                                                                                            >  3:130952/1‑136 (MQ=255)
                       tttGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTAcc                                                                                                                              >  2:466264/1‑148 (MQ=255)
                            tCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAAGTCTTGTAAGGTTATGGGGTACCTGCgg                                                                                                                         >  7:254909/1‑148 (MQ=255)
                                 tATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTAcc                                                                                                                              >  3:23379/1‑138 (MQ=255)
                                 tATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTAcc                                                                                                                              <  4:23379/138‑1 (MQ=255)
                                  aTTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCTCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTa                                                                                                                    >  4:281099/1‑147 (MQ=255)
                                  aTTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTa                                                                                                                    >  5:166665/1‑147 (MQ=255)
                                  aTTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGCAAGGTTATGGGGTACCTGCGGCCTTa                                                                                                                    >  4:369155/1‑147 (MQ=255)
                                  aTTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAAGTCTTGTAAGGtt                                                                                                                                        <  6:437062/127‑1 (MQ=255)
                                  aTTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAAGTCTTGTAAGGtt                                                                                                                                        >  5:437062/1‑127 (MQ=255)
                                   tttCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGTTCTTGTAAGGTTATGGGGTACCTGCGGCCTTaa                                                                                                                   <  8:269161/147‑1 (MQ=255)
                                   tttCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTaaa                                                                                                                  <  5:380139/148‑1 (MQ=255)
                                   tttCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAAGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTaaa                                                                                                                  <  2:236080/148‑1 (MQ=255)
                                       aGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATa                                                                                                               <  6:166665/147‑1 (MQ=255)
                                       aGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATa                                                                                                               <  1:287181/147‑1 (MQ=255)
                                        tgATAGAAATGGAGAAGGAGAATACCAATGGCTTTCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTACTTTAGTAAAACCACAAGTATTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAac                                                                                                             <  8:254909/147‑1 (MQ=255)
                                         taTAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAaca                                                                                                            <  5:32884/147‑1 (MQ=255)
                                            aGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCta                                                                                                                            <  7:144755/129‑2 (MQ=255)
                                            aGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCta                                                                                                                            >  8:144755/1‑128 (MQ=255)
                                             gAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGCAAGGTTATGGGGTACCTGCGGCCTTAAAATAACAccc                                                                                                         <  2:12181/147‑1 (MQ=255)
                                                 tGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGTTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATa                                                                                                               >  1:108117/1‑137 (MQ=255)
                                                 tGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGTTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATa                                                                                                               <  2:108117/137‑1 (MQ=255)
                                                     cAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACAt                                                                                                >  3:272408/1‑148 (MQ=255)
                                                      aaGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACAt                                                                                                >  7:173918/1‑147 (MQ=255)
                                                       aGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAAACACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATca                                                                                              >  5:2687/1‑148 (MQ=255)
                                                         gagaATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATca                                                                                              >  1:255454/1‑146 (MQ=255)
                                                         gagaATACAAATGGATATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTAAGTAAAACCAAAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATcaca                                                                                            >  1:326557/1‑148 (MQ=255)
                                                            aaTACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGaa                                                                                         <  6:261682/148‑1 (MQ=255)
                                                             aTACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGCAAGGTTATGGTGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGaa                                                                                         <  6:78386/147‑1 (MQ=255)
                                                              aatggctatcgctgtcaggacctgaaagtaaaggcttttaaacagagaggttttcatagggtgtcctttagtaaaaccacagACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCa                                                                                                        <  6:316404‑M2/49‑1 (MQ=255)
                                                              aatggctatcgctgtcaggacctgaaagtaaaggcttttaaacagagaggttttcatagggtgtcctttagtaaaaccacagACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCa                                                                                                        >  5:316404‑M2/83‑131 (MQ=255)
                                                               aCCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATg                                                                                      >  7:222940/1‑148 (MQ=255)
                                                                 cAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTa                                                                                    <  2:238317/148‑1 (MQ=255)
                                                                  gctatcgctgtcaggacctgaaagtaaaggcttttaaacagagaggttttcatagggtgtcctttagtaaaaccacagACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTa                                                                                    <  7:483062‑M2/69‑1 (MQ=255)
                                                                    tGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCt                                                                                 >  3:242086/1‑148 (MQ=255)
                                                                    tGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTAATAGGGTGTCCTTTATTAAAACCACAGGTCTTGTAAGGTTATGGGGTACATACGGCCTTAAAATAAAACACAGACAACATCACACAAATGTACCt                                                                                 >  3:129515/1‑148 (MQ=255)
                                                                     ggCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTg                                                                                >  4:354844/1‑148 (MQ=255)
                                                                     ggCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGGCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTg                                                                                >  2:381274/1‑148 (MQ=255)
                                                                         aTCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTAt                                                                                                                                      <  1:218149/90‑1 (MQ=255)
                                                                         aTCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTAt                                                                                                                                      >  2:218149/1‑90 (MQ=255)
                                                                          tCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGacaac                                                                                                  >  2:210163/1‑125 (MQ=255)
                                                                          tCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGAGTCCATTACTAAAACCACAGGTCATGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTCGATCa                                                                           >  6:336138/1‑148 (MQ=255)
                                                                           cGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTg                                                                                                                                               <  2:125525/79‑1 (MQ=255)
                                                                           cGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTg                                                                                                                                               >  1:125525/1‑79 (MQ=255)
                                                                                tCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAAACACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTg                                                                     <  6:2687/148‑1 (MQ=255)
                                                                                  aGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGa                                                                    <  1:276414/147‑1 (MQ=255)
                                                                                      ccTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGCAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACaa                                                               <  3:369155/148‑1 (MQ=255)
                                                                                        tGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACaaa                                                              >  1:128273/1‑147 (MQ=255)
                                                                                        tGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACaaa                                                              >  4:302265/1‑147 (MQ=255)
                                                                                         gAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACaaaa                                                             >  5:195025/1‑147 (MQ=255)
                                                                                          aaaGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAg                                                           >  4:183585/1‑148 (MQ=255)
                                                                                          aaaGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACaaa                                                              <  6:195025/145‑1 (MQ=255)
                                                                                               aaaGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTa                                                                                    <  7:100873/118‑1 (MQ=255)
                                                                                               aaaGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTa                                                       >  3:404678/1‑147 (MQ=255)
                                                                                               aaaGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGCAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTa                                                       >  6:282164/1‑147 (MQ=255)
                                                                                               aaaGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACAACAGGTCTTGTAAGGTTTTGGGGTACCTACGGCCTTAAAATAAAACCCAGACAAAATCACAGAAATGTACCTGGAAAATAATTGAAACAAAAAGGATa                                                       >  3:411059/1‑147 (MQ=255)
                                                                                                aaGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTaa                                                      >  8:469165/1‑147 (MQ=255)
                                                                                                    cTTTTAAACTGAGAGGTTTTCATCTGGTGTCATTTAGTAAATCCACTGTTATTGTAAGGTTTTGGGGTAACTGAGGAATTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAATAAGTGTAAAAAGt                                                 <  5:308352/148‑1 (MQ=255)
                                                                                                      tttAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGCAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGttt                                               >  8:36525/1‑148 (MQ=255)
                                                                                                       ttAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAAAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGttt                                               >  7:117483/1‑147 (MQ=255)
                                                                                                       ttAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAAGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAAGGTACCTGGATCATAATTGAAACAAAAAGGTAAAAAAGTTTg                                              >  6:66625/1‑148 (MQ=255)
                                                                                                        tAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTg                                              <  7:171145/147‑1 (MQ=255)
                                                                                                        tAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTg                                              <  7:469165/147‑1 (MQ=255)
                                                                                                            cAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCGCAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTg                                         >  8:156996/1‑148 (MQ=255)
                                                                                                              gagagGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGTTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGc                                       >  1:255630/1‑148 (MQ=255)
                                                                                                               agagGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATcac                                                                                             <  5:106851/93‑1 (MQ=255)
                                                                                                               agagGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATcac                                                                                             >  6:106851/1‑93 (MQ=255)
                                                                                                                  ggTTTTCATAGGGTGTCCTTTAGTAAAACCAAAGGTCTTGTAAGGTTATCGGGTACCTGCGGCCTTAAAATAACACCCAGACATCATCACAGAAATGAACATGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTTAGCTGGCTCg                                    >  4:259509/1‑147 (MQ=255)
                                                                                                                     tttCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGa                                             >  8:467164/1‑135 (MQ=255)
                                                                                                                        cATAGGGTGTCCTTTAGTAAAACCACAGATCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACaaa                              <  5:444010/147‑1 (MQ=255)
                                                                                                                         aTAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAAcc                            <  7:66098/148‑1 (MQ=255)
                                                                                                                                gtCCTTTAGTAAAACCACAGATCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACtt                     >  7:343677/1‑148 (MQ=255)
                                                                                                                                 tCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACtt                     <  6:342501/147‑1 (MQ=255)
                                                                                                                                  ccTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACtttt                   >  4:200913/1‑148 (MQ=255)
                                                                                                                                  ccTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACtttt                   >  2:425492/1‑148 (MQ=255)
                                                                                                                                  ccTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACtttt                   >  5:456443/1‑148 (MQ=255)
                                                                                                                                    tttAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACttttt                  >  7:343009/1‑147 (MQ=255)
                                                                                                                                    tttAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACttttt                  >  7:470022/1‑147 (MQ=255)
                                                                                                                                    tttAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACttttt                  >  2:101108/1‑147 (MQ=255)
                                                                                                                                     ttAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTGCCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTa                 >  8:286977/1‑147 (MQ=255)
                                                                                                                                     ttAGTAAAACCACAGGTCCTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTa                 >  4:281374/1‑147 (MQ=255)
                                                                                                                                      tAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTa                 <  3:354844/146‑1 (MQ=255)
                                                                                                                                         tAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGtt            >  6:10259/1‑148 (MQ=255)
                                                                                                                                            aaCCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTg         >  1:20974/1‑148 (MQ=255)
                                                                                                                                              ccACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTa       <  6:456443/148‑1 (MQ=255)
                                                                                                                                                   ggTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAg  >  8:188993/1‑148 (MQ=255)
                                                                                                                                                   ggTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAg  >  1:393716/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GAAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAG  >  NC_000913/3683325‑3683619

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