Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 911,782 G→A 27.9% V123V (GTC→GTT)  hcr ← HCP oxidoreductase, NADH‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913911,7820GA27.9% 110.7 / 47.3 68V123V (GTC→GTT) hcrHCP oxidoreductase, NADH‑dependent
Reads supporting (aligned to +/- strand):  ref base G (28/21);  new base A (9/10);  total (37/31)
Fisher's exact test for biased strand distribution p-value = 5.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

GTTATTTTCTGCCACCAGCGTTACCGGATAGTTACGCCACTCATCGGCGAAAATAACATCCTGCGGCGTACGCACGTTGTAGATCACCCGCACATCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTGACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCGAAAGCCAGAGATAATCACCGCGTTTTACATCGCGCGTCAGCCACTGGGAGCCGACACCGTCATCAATCCG  >  NC_000913/911635‑911928
                                                                                                                                                   |                                                                                                                                                  
gTTATTTTCTGCCACCAGCGTTACCGGATAGTTACGCCACTCATCGGCGAAAATAACATCCTGCGGCGTACGCACGTTGTAGATCACCCGCACATCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTg                                                                                                                                                    <  3:244297/148‑1 (MQ=255)
  tATTTTCTGCCACCAGCGTTACCGGATAGTTACGCCACTCATCGGCGAAAATAACATCCTGCGGCGTACGCACGTTGTAGATCACCCGCACATCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTGAc                                                                                                                                                  >  1:263349/1‑148 (MQ=255)
        cTGCCACCAGCGTTACCGGATAGTTACGCCACTCATCGGCGAAAATAACATCCTGCGGCGTACGCACGTTGTAGATCACCCGCACATCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTGAcgccgc                                                                                                                                             <  6:32400/147‑1 (MQ=255)
        cTGCCACCAGCGTTACCGGATAGTTACGCCACTCATCGGCGAAAATAACATCCTGCGGCGTACGCACGTTGTAGATCACCCGCACATCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTGAcgccgc                                                                                                                                             <  2:263349/147‑1 (MQ=255)
               caGCGTTACCGGATAGTTACGCCACTCATCGGCGAAAATAACATCCTGCGGCGTACGCACGTTGTAGATCACCCGAACCTCGGCCTGTGGCCGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTGACGCCGCAGCCTgcc                                                                                                                                     >  7:248949/1‑148 (MQ=255)
                 gCGTTACCGGATAGTTACGCCACTCATCGGCGAAAATAACATCCTGCGGCGTACGCACGTTGTAGATCACCCGCACATCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTGACGCCGCAGCCTgccgc                                                                                                                                   >  6:99816/1‑148 (MQ=255)
                     tACCGGATAGTTACGCCACTCATCGGCGAAAATAACATCCTGCGGCGTACGCACGTTGTAGATCACCCGCACATCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTg                                                                                                                                                    <  3:58742/127‑1 (MQ=255)
                     tACCGGATAGTTACGCCACTCATCGGCGAAAATAACATCCTGCGGCGTACGCACGTTGTAGATCACCCGCACATCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTg                                                                                                                                                    >  4:58742/1‑127 (MQ=255)
                        cGGATAGTTACGCCACTCATCGGCGAAAATAACATCCTGCGGCGTACGCACGTTGTAGATCACCCGCACATCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTGACGCCGCAGCCTGCCGCCAGCaa                                                                                                                             >  3:185021/1‑147 (MQ=255)
                         ggATAGTTACGCCACTCATCGGCGAAAATAACATCCTGCGGCGTACGCACGTTGTAGATCACCCGCACATCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTGACGCCGCAGCCTGCCGCCAGCaa                                                                                                                             >  4:224339/1‑146 (MQ=255)
                         ggATAGTTACGCCACTCATCGGCGAAAATAACATCCTGCGGCGTACGCACGTTGTAGATCACCCGCACATCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTGACGCCGCAGCCTGCCGCCAGCAACa                                                                                                                           <  2:196305/148‑1 (MQ=255)
                                      aCTCATCGGCGAAAATAACATCCTGCGGCGTACGCACGTTGTAGATCACCCGCACATCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTAACGCCGCAGCCTGCCGCCAGCAACAGGAATTTAtcttc                                                                                                              >  2:155513/1‑148 (MQ=255)
                                       cTCATCGGCGAAAATAACATCCTGCGGCGTACGCACGTTGTAGATCACCCGCACATCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTGACGCCGCAGCCTGCCGCCAGCAACAGGAATTTAtcttct                                                                                                             <  4:185021/148‑1 (MQ=255)
                                            cGGCGAAAATAACATCCTGCGGCGTACGCACGTTGTAGATCACCCGCACATCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTGACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCtt                                                                                                         >  3:151504/1‑147 (MQ=255)
                                            cGGCGAAAATAACATCCTGCGGCGTACGCACGTTGTAGATCACCCGCACATCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTAACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCtt                                                                                                         >  4:178404/1‑147 (MQ=255)
                                              gCGAAAATAACATCCTGCGGCGTACGCACGTTGTAGATCACCCGCACATCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTGACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTAt                                                                                                      >  5:83952/1‑148 (MQ=255)
                                                 aaaaTAACATCCTGCGGCGTACGCACGTTGTAGATCACCCGCACATCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTGACGCCGCAGCCTGCCGCCAGCAAc                                                                                                                            <  2:104828/123‑1 (MQ=255)
                                                 aaaaTAACATCCTGCGGCGTACGCACGTTGTAGATCACCCGCACATCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTGACGCCGCAGCCTGCCGCCAGCAAc                                                                                                                            >  1:104828/1‑123 (MQ=255)
                                                  aaaTAACATCCTGCGGCGTACGCACGTTGTAGATCACCCGCACATCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTGACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTAtcgtc                                                                                                  <  6:315147/148‑1 (MQ=255)
                                                     tAACATCCTGCGGCGTACGCACGTTGTAGATCACCCGCACATCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTGACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCa                                                                                               <  5:99816/148‑1 (MQ=255)
                                                       aCATCCTGCGGCGTACGCACGTTGTAGATCACCCGCACATCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTGACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAgg                                                                                             >  3:90144/1‑148 (MQ=255)
                                                          tCCTGCGGCGTACGCACGTTGTAGATCACCCGCACATCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTGACGCCGCAGCCTGCCGCCAGCCACAGGAAGTTATCTTCTGCTTTATCGTCGCAGGTa                                                                                           >  5:123619/1‑147 (MQ=255)
                                                               cggcgTACGCACGTTGTAGATCACCCGCACATCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTGACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTcc                                                                                     <  2:208797/148‑1 (MQ=255)
                                                                 gcgTACGCACGTTGTAGATCACCCGCACATCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTg                                                                                                                                                    >  8:77165/1‑83 (MQ=255)
                                                                 gcgTACGCACGTTGTAGATCACCCGCACATCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTg                                                                                                                                                    <  7:77165/83‑1 (MQ=255)
                                                                     actcACGTTGTAGATCACCCGCACATCGGCCTGTGGACGGTTCTTCCCCATCCAGCGACGCATCGACATAATCGGCGTGACGCCGCAGACTGACGCCAGCAACAGGAATTTAACTTCTGATTTATCGTCGCAGTTAAATTCACCCATc                                                                               >  4:2/4‑148 (MQ=255)
                                                                       gCACGTTGTAGATCACCCGCACATCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTGACGCCGCAGCCTGCCGCCAGCAACAGTAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATcgc                                                                             <  5:303643/148‑1 (MQ=255)
                                                                         aCGTTGTAGATCACCCGCACATCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTGACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGt                                                                           <  6:83952/148‑1 (MQ=255)
                                                                             tGTAGATCACCCGCACATCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTGACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCg                                                                        >  7:228796/1‑147 (MQ=255)
                                                                               tAGATCACCCGCACATCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTGACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCGaa                                                                      <  1:908/147‑1 (MQ=255)
                                                                                aGATCACCCGCACATCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTGACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCGAAAg                                                                    >  2:219197/1‑148 (MQ=255)
                                                                                  aTCACCCGCACATCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGCCGCATCGACATAATCTGCGTGAATCCGCAGCCTGCCGCCAGAAACAGGAATTTATATTCTGATTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCGAAAga                                                                   >  4:165493/1‑146 (MQ=255)
                                                                                    cACCCGCACATCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTACCGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCGAAAg                                                                    >  8:130867/1‑144 (MQ=255)
                                                                                    cACCCGCACATCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTAACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCGAAAg                                                                    <  7:130867/144‑1 (MQ=255)
                                                                                      cccGCACATCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTAACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCGAAAGCCagag                                                              <  3:178404/148‑1 (MQ=255)
                                                                                        cGCACATCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTGACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCGAAAGCCAgaga                                                             >  7:26438/1‑147 (MQ=255)
                                                                                        cGCACATCGGCCTGTGGACCGTTCTTCGCAACCCAGCGACGCATCGACATAATCGGCGTAACGCCGCAGCCTGCCGCCATCAACACGAATTTAACTTCTCCTTTATCGACGCAGGTAAATTCCCCCAACGCGTCCGAAAGCCAgaga                                                             >  8:239803/1‑147 (MQ=255)
                                                                                          cacaTCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTGACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCGAAAGCCAGAGATaa                                                          >  2:1530/1‑148 (MQ=255)
                                                                                            caTCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTAACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCATGTAAATTCCCCCATCGCGTCCGAAAGCCAGAGATAATc                                                        <  7:239803/148‑1 (MQ=255)
                                                                                             aTCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTAACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCGAAAGCCAGAGATAATCa                                                       <  1:43697/148‑1 (MQ=255)
                                                                                                     gtgGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTAACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCGAAAGCCAGAGATAATCACCGCGttt                                               >  6:36606/1‑148 (MQ=255)
                                                                                                       ggACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTGACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCGAAAGCCAGAGATAATCACCGCGtttt                                              >  7:35773/1‑147 (MQ=255)
                                                                                                         aCGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTAACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGt                                                                           <  3:215372/116‑1 (MQ=255)
                                                                                                         aCGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTAACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGt                                                                           >  4:215372/1‑116 (MQ=255)
                                                                                                         aCGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTAACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCGAAAGCCAGAGATAATCACCGCGTTTTACa                                           >  6:141691/1‑148 (MQ=255)
                                                                                                          cGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTGACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCGAAAGCCAGAGATAATCACCGCGTTTTACAt                                          >  6:185434/1‑148 (MQ=255)
                                                                                                          cGGTTCTTAGCAAGCCAGCGACGATTCGACATACTCGGCGTAACGCAGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCTAAAGCCTGAGATAATCACCGCGTTTTACAt                                          <  4:292539/148‑1 (MQ=255)
                                                                                                              tcttcGCAAGCCAGCGACGCATCGACATAATCGGCGTAACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCGAAAGCCAGAGATAATCACCGCGTTTTACATcgcg                                      >  5:185980/1‑148 (MQ=255)
                                                                                                               cttcGCAAGCCAGCGACGCATCGACATAATCGGCGTAACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCGAAAGCCAGAGATAATCACCGCGTTTTACATcgcg                                      <  1:155513/147‑1 (MQ=255)
                                                                                                               cttcGCAAGCCAGCGACGCATAGACTTATTACGCGTGACGCCTCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCGCATCGCGGCCGAAAGCCAGAGTTAATCACCGCGTTTTACATcgcg                                      <  3:2/147‑1 (MQ=255)
                                                                                                                   gCAAGCCAGCGACGCATCGACATAATCGGCGTGACGCCGCAGCCTGCAGCCAGCAACACGAATTTATCTTCTGCTTTATCGGCGCAGGTAAATTCCCCCATCGCGTCCGAAAGCCAGAGATAATCACCGCGTTTTACATCGCGCGTc                                  <  8:248949/147‑1 (MQ=255)
                                                                                                                    cAAGCCAGCGACGCATCGACATAATCGGCGTGACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCGAACGCCAGAGATAATCACCGCGTTTTACATCGCGCGTCAg                                >  6:208090/1‑148 (MQ=255)
                                                                                                                    cAAGCCAGCGACGCATCGACATAATCGGCGTGACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCGAAAGCCAGAGATAATCACCGCGTTTTACATCGCGCGTCAg                                <  8:26438/148‑1 (MQ=255)
                                                                                                                     aaGCCAGCGACGCATCGACATAATCGGCGTGACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCGAAAGCCAGAGATAATCACCGCGTTTTACATCGCGCGTCAGc                               >  5:84479/1‑148 (MQ=255)
                                                                                                                     aaGCCAGCGACGCATCGACATAATCGGCGTGACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCGAAAGCCAGAGATAATCACCGCGTTTTACATCGCGCGTCAGc                               >  7:62045/1‑148 (MQ=255)
                                                                                                                        ccAGCGACGCATCGACATAATCGGCGTGACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCGAAAGCCAGAGATAATCACCGCGTTTTACATCGCGCGTCAGCCAc                            <  8:35773/148‑1 (MQ=255)
                                                                                                                          aGCGACGCATCGACATAATCGGCGTGACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCGAAAGCCAGAGATAATCACCGCGTTTTACATCGCGCGTCAGCCACTg                          <  1:219197/148‑1 (MQ=255)
                                                                                                                          aGCGACGCATCGACATAATCGGCGTAACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCGAAAGCCAGAGATAATCACCGCGTTTTACATCGCGCGTCAGCCACTg                          <  5:141691/148‑1 (MQ=255)
                                                                                                                            cGACGCATCGACATAATCGGCGTAACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCGAAAGCCAGAGATAATCACCGCGTTTTACATCGCGCGTCAGCCACTggg                        <  7:315521/148‑1 (MQ=255)
                                                                                                                                 cATCGACATAATCGGCGTGACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCGAAAGCCAGAGATAATCACCGCGTTTTACATCGCGCGTCAGCCACTGGGAGCCg                   <  1:1530/148‑1 (MQ=255)
                                                                                                                                      aCATAATCGGCGTGACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCGAAAGCCAGAGATAATCACCGCGTTTTACATCGCGCGTCAGCCACTGGGAGCCGACAcc              >  4:236208/1‑148 (MQ=255)
                                                                                                                                        aTAATCGGCTTGAAGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGGCCGAAAGCCAGAGATAATCACCGCGTTTTACATCGCGCGTCAGCCACTGGGAGCCGACACCGt            <  6:123619/148‑1 (MQ=255)
                                                                                                                                        aTAATCGGCGTGACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCGAAAGCCAGAGATAATCACCGCGTTTTACATCGCGCGTCAGCCACTGGGAGCCGACACCGt            >  7:309361/1‑148 (MQ=255)
                                                                                                                                         tAATCGGCGTGACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCGAAAGCCAGAGATAATCACCGCGTTTTACATCGCGCGTCAGCCACTGGGAGCCGACACCGtc           <  1:220105/148‑1 (MQ=255)
                                                                                                                                         tAATCGGCGTGACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCGAAAGCCAGAGATAATCACCGCGTTTTACATCGCGCGTCAGCCACTGGGAGCCGACACCGt            >  1:247727/1‑147 (MQ=255)
                                                                                                                                         tAATCGGCGTGACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCGAAAGCCAGAGATAATCACCGCGTTTTACATCGCGCGTCAGCCACTGGGAGCCGACACCGt            >  8:100418/1‑147 (MQ=255)
                                                                                                                                         tAATCGGCGTGACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTATGCTTTATCGTCGACGGTAAATTCCCCCATCGCGTCCGAAAGCCAGAGATAATCCCCGCGTTTTACATCGCGCGTCAGCCACTGGGAGCCGACACCGt            >  5:74787/1‑147 (MQ=255)
                                                                                                                                         tAATCGGCGGGACGCCGCAGCCGGCCGCCAGCCACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGGCCGAAAGCCAGAGATAATCACCGCGTTTTACATCGCGCGTCAGCCACTGGGAGCCGACACCGtc           <  6:74787/148‑1 (MQ=255)
                                                                                                                                          aaTCGGCGTGACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCGAAAGCCAGAGATAATCACCGCGTTTTACATCGCGCGTCAg                                <  3:108264/126‑1 (MQ=255)
                                                                                                                                          aaTCGGCGTGACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCGAAAGCCAGAGATAATCACCGCGTTTTACATCGCGCGTCAg                                >  4:108264/1‑126 (MQ=255)
                                                                                                                                          aaTCGGCGTAACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCGAAAGCCAGAGATAATCACCGCGTTTTACATCGCGCGTCAGCCACTGGGAGCCGACACCGtca          >  5:72970/1‑148 (MQ=255)
                                                                                                                                               gCGTAACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCGAAAGCCAGAGATAATCACCGCGTTTTACATCGCGCGTCAGCCACTGGGAGCCGACACCGTCATCAAt     <  6:185980/148‑1 (MQ=255)
                                                                                                                                                  tGACGCCGCCGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCGAACGCCAGAGATAATCACCGCGTTTTACATCGCGCGTCAGCCACTGGGAGCCGACACCGTCATCAATccg  >  2:81778/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTTATTTTCTGCCACCAGCGTTACCGGATAGTTACGCCACTCATCGGCGAAAATAACATCCTGCGGCGTACGCACGTTGTAGATCACCCGCACATCGGCCTGTGGACGGTTCTTCGCAAGCCAGCGACGCATCGACATAATCGGCGTGACGCCGCAGCCTGCCGCCAGCAACAGGAATTTATCTTCTGCTTTATCGTCGCAGGTAAATTCCCCCATCGCGTCCGAAAGCCAGAGATAATCACCGCGTTTTACATCGCGCGTCAGCCACTGGGAGCCGACACCGTCATCAATCCG  >  NC_000913/911635‑911928

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