Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,683,497 G→A 43.6% intergenic (‑50/+133) 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,4970GA43.6% 65.5 / 163.1 133intergenic (‑50/+133)dctA/yhjKC4‑dicarboxylic acid, orotate and citrate transporter/cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (40/35);  new base A (34/24);  total (74/59)
Fisher's exact test for biased strand distribution p-value = 5.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.53e-02

Only 100 of 136 total aligned reads displayed.

TCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAG  >  NC_000913/3683353‑3683644
                                                                                                                                                |                                                                                                                                                   
tCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAAGTCTTGTAAGGTTATGGGGTACCTGCgg                                                                                                                                                  >  7:254909/1‑148 (MQ=255)
      aTTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCTCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTa                                                                                                                                             >  4:281099/1‑147 (MQ=255)
      aTTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTa                                                                                                                                             >  5:342501/1‑147 (MQ=255)
      aTTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTa                                                                                                                                             >  5:166665/1‑147 (MQ=255)
      aTTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGCAAGGTTATGGGGTACCTGCGGCCTTa                                                                                                                                             >  4:369155/1‑147 (MQ=255)
       tttCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGTTCTTGTAAGGTTATGGGGTACCTGCGGCCTTaa                                                                                                                                            <  8:269161/147‑1 (MQ=255)
       tttCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTaaa                                                                                                                                           <  5:380139/148‑1 (MQ=255)
       tttCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAAGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTaaa                                                                                                                                           <  2:236080/148‑1 (MQ=255)
        ttCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTaaaa                                                                                                                                          >  7:108099/1‑148 (MQ=255)
          cAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAAt                                                                                                                                         <  6:19656/147‑1 (MQ=255)
           aGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATa                                                                                                                                        <  4:244295/147‑1 (MQ=255)
            tgATAGAAATGGAGAAGGAGAATACCAATGGCTTTCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTACTTTAGTAAAACCACAAGTATTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAac                                                                                                                                      <  8:254909/147‑1 (MQ=255)
             taTAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAaca                                                                                                                                     <  5:32884/147‑1 (MQ=255)
                 gAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGCAAGGTTATGGGGTACCTGCGGCCTTAAAATAACAccc                                                                                                                                  <  2:12181/147‑1 (MQ=255)
                     tGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGTTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATa                                                                                                                                        >  1:108117/1‑137 (MQ=255)
                          aaGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACAt                                                                                                                         >  7:173918/1‑147 (MQ=255)
                           aGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAAACACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATca                                                                                                                       >  5:2687/1‑148 (MQ=255)
                             gagaATACAAATGGATATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTAAGTAAAACCAAAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATcaca                                                                                                                     >  1:326557/1‑148 (MQ=255)
                                 aTACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGCAAGGTTATGGTGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGaa                                                                                                                  <  6:78386/147‑1 (MQ=255)
                                  aatggctatcgctgtcaggacctgaaagtaaaggcttttaaacagagaggttttcatagggtgtcctttagtaaaaccacagACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCa                                                                                                                                 <  6:316404‑M2/49‑1 (MQ=255)
                                   aCCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGATCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATg                                                                                                               >  6:444010/1‑148 (MQ=255)
                                      gctatcgctgtcaggacctgaaagtaaaggcttttaaacagagaggttttcatagggtgtcctttagtaaaaccacagACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTa                                                                                                             <  7:483062‑M2/69‑1 (MQ=255)
                                      aaTGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAAACACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTAc                                                                                                            >  2:185255/1‑148 (MQ=255)
                                        tGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCt                                                                                                          >  3:242086/1‑148 (MQ=255)
                                         ggCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTg                                                                                                         >  4:354844/1‑148 (MQ=255)
                                              tCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGacaac                                                                                                                           >  2:210163/1‑125 (MQ=255)
                                              tCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGacaac                                                                                                                           <  1:210163/125‑1 (MQ=255)
                                                    tCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAAACACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTg                                                                                              <  6:2687/148‑1 (MQ=255)
                                                      aGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGa                                                                                             <  1:276414/147‑1 (MQ=255)
                                                       ggACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGaaa                                                                                           >  7:200183/1‑148 (MQ=255)
                                                            tGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACaaa                                                                                       >  1:128273/1‑147 (MQ=255)
                                                             gAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACaaaa                                                                                      >  5:195025/1‑147 (MQ=255)
                                                              aaaGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAg                                                                                    >  4:183585/1‑148 (MQ=255)
                                                              aaaGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGAGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACaaaaa                                                                                     <  3:281099/147‑1 (MQ=255)
                                                              aaaGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACaaa                                                                                       <  6:195025/145‑1 (MQ=255)
                                                                 gTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGt                                                                                  >  4:333306/1‑147 (MQ=255)
                                                                   aaaGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTa                                                                                                             <  7:100873/118‑1 (MQ=255)
                                                                   aaaGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTa                                                                                                             >  8:100873/1‑118 (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)
                                                                              aaCAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAg                                                                     <  8:222940/147‑1 (MQ=255)
                                                                                cAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCGCAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTg                                                                  >  8:156996/1‑148 (MQ=255)
                                                                                  gagagGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGTTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGc                                                                >  1:255630/1‑148 (MQ=255)
                                                                                   agagGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATcac                                                                                                                      >  6:106851/1‑93 (MQ=255)
                                                                                         tttCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGa                                                                      >  8:467164/1‑135 (MQ=255)
                                                                                         tttCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGa                                                                      <  7:467164/135‑1 (MQ=255)
                                                                                                    gtCCTTTAGTAAAACCACAGATCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACtt                                              >  7:343677/1‑148 (MQ=255)
                                                                                                     tCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACtt                                              <  6:342501/147‑1 (MQ=255)
                                                                                                      ccTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACtttt                                            >  2:425492/1‑148 (MQ=255)
                                                                                                      ccTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACtttt                                            >  4:200913/1‑148 (MQ=255)
                                                                                                      ccTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACtttt                                            >  5:456443/1‑148 (MQ=255)
                                                                                                        tttAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACttttt                                           >  1:60667/1‑147 (MQ=255)
                                                                                                        tttAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACttttt                                           >  7:343009/1‑147 (MQ=255)
                                                                                                        tttAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACttttt                                           >  7:470022/1‑147 (MQ=255)
                                                                                                        tttAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACttttt                                           >  2:101108/1‑147 (MQ=255)
                                                                                                         ttAGTAAAACCACAGGTCCTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTa                                          >  4:281374/1‑147 (MQ=255)
                                                                                                          tAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTa                                          <  3:354844/146‑1 (MQ=255)
                                                                                                          aaaaccacagACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTaaa                                        >  2:79089‑M2/11‑148 (MQ=255)
                                                                                                             tAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGtt                                     >  6:10259/1‑148 (MQ=255)
                                                                                                                aaCCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTg                                  >  1:28324/1‑148 (MQ=255)
                                                                                                                aaCCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTg                                  >  1:20974/1‑148 (MQ=255)
                                                                                                                  ccACAGGTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTa                                <  6:456443/148‑1 (MQ=255)
                                                                                                                   cacaGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTaa                               >  4:66763/1‑148 (MQ=255)
                                                                                                                       ggTCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAg                           >  8:188993/1‑148 (MQ=255)
                                                                                                                        gTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGt                          <  1:425492/148‑1 (MQ=255)
                                                                                                                        gTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGt                          <  7:399203/148‑1 (MQ=255)
                                                                                                                        atCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGt                          >  7:118353/2‑148 (MQ=255)
                                                                                                                         tCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGtt                         <  1:185255/148‑1 (MQ=255)
                                                                                                                         tCTTGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGtt                         <  3:333306/148‑1 (MQ=255)
                                                                                                                         tCTTGCAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGtt                         <  5:282164/148‑1 (MQ=255)
                                                                                                                          cTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACAACCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGATTGAGCTGGCTCGAACAAACCAGCACTTTTTAAAGTTTTTTAATCAGttt                        >  4:294540/1‑148 (MQ=255)
                                                                                                                            tGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTgg                      >  2:205721/1‑148 (MQ=255)
                                                                                                                            tGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCACACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTg                       >  1:407252/1‑147 (MQ=255)
                                                                                                                            tGTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTgg                      >  8:1163/1‑148 (MQ=255)
                                                                                                                             gTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTgg                      <  1:381274/147‑1 (MQ=255)
                                                                                                                             gTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTgg                      <  8:343677/147‑1 (MQ=255)
                                                                                                                             gTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCACACAAACCAGCACTTTTTAAAGTTTTGTAATCAGttt                        <  2:407252/145‑1 (MQ=255)
                                                                                                                             gTAAGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTg                       <  2:255454/146‑1 (MQ=255)
                                                                                                                              tAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTggg                     <  3:200913/147‑1 (MQ=255)
                                                                                                                               aaGGTTATGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTa                                <  8:386053/135‑1 (MQ=255)
                                                                                                                                 ggTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAg                 >  2:288181/1‑148 (MQ=255)
                                                                                                                                     aTGGGGTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTAc             >  7:480594/1‑148 (MQ=255)
                                                                                                                                          gTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAATAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTttc        <  2:326557/148‑1 (MQ=255)
                                                                                                                                          gTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTttc        <  2:255630/148‑1 (MQ=255)
                                                                                                                                          gTACCTACGGCCTTAAAATAACACCCATACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGATGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTttc        >  3:246810/1‑148 (MQ=255)
                                                                                                                                          gTACCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTttc        >  5:15771/1‑148 (MQ=255)
                                                                                                                                            aCCTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGt                          <  7:221155/128‑1 (MQ=255)
                                                                                                                                             ccTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTttcttc     <  8:161940/148‑1 (MQ=255)
                                                                                                                                             ccTACGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTttctt      <  5:10259/147‑1 (MQ=255)
                                                                                                                                              cTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTcc    <  2:20974/148‑1 (MQ=255)
                                                                                                                                              cTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTcc    <  2:28324/148‑1 (MQ=255)
                                                                                                                                              cTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTAACTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTcc    <  1:79089/148‑1 (MQ=255)
                                                                                                                                               tGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTttcttc     >  1:340675/1‑146 (MQ=255)
                                                                                                                                                gCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCt                                                               >  7:395978/1‑87 (MQ=255)
                                                                                                                                                gCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAg  >  1:252284/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAG  >  NC_000913/3683353‑3683644

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