Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 134,166 A→G 5.0% Q851R (CAG→CGG)  acnB → aconitate hydratase 2; aconitase B; 2‑methyl‑cis‑aconitate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913134,1660AG 5.0% 308.6 / 4.7 101Q851R (CAG→CGG) acnBaconitate hydratase 2; aconitase B; 2‑methyl‑cis‑aconitate hydratase
Reads supporting (aligned to +/- strand):  ref base A (49/46);  new base G (3/2);  total (52/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.38e-01

Only 100 of 105 total aligned reads displayed.

GAATGTCTTCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTTCTCC  >  NC_000913/134017‑134315
                                                                                                                                                     |                                                                                                                                                      
gAATGTCTTCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCa                                                                                                                                                        >  1:378734/1‑150 (MQ=255)
  aTGTCTTCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGCCGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACATACCGTTATCTGAACTTCAACCATCTTTCTCAGt                                                                                                                                                      >  2:758980/1‑150 (MQ=255)
   tGTCTTCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTa                                                                                                                                                     >  1:1679836/1‑150 (MQ=255)
    gTCTTCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTac                                                                                                                                                    >  2:1149609/1‑150 (MQ=255)
        tCCTGGCTTCTGCGGAACTGGCGGCTGGTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCg                                                                                                                                                <  1:1450911/150‑1 (MQ=255)
         ccTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCga                                                                                                                                               >  1:1382364/1‑150 (MQ=255)
         ccTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCga                                                                                                                                               >  1:1383372/1‑150 (MQ=255)
          cTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCgag                                                                                                                                              >  1:748707/1‑150 (MQ=255)
             gctgctGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCAGTAAGAGTACCATACATAAGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGaaa                                                                                                                                           <  1:28263/150‑1 (MQ=255)
              cTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAg                                                                                                                                          >  2:1561917/1‑150 (MQ=255)
                tCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCTGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGcc                                                                                                                                        <  2:1676086/150‑1 (MQ=255)
                tCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCGGTGTACTGAGAAAGct                                                                                                                                        <  1:1462913/150‑2 (MQ=255)
                 cTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCGGTGTACTGAGAAAGctg                                                                                                                                       <  1:701284/150‑3 (MQ=255)
                 cTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCg                                                                                                                                       >  2:1718569/1‑150 (MQ=255)
                      gAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATggg                                                                                                                                  <  2:1079966/150‑1 (MQ=255)
                      gAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATggg                                                                                                                                  >  1:810105/1‑150 (MQ=255)
                        aCTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGt                                                                                                                                <  1:1110621/150‑1 (MQ=255)
                         cTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTg                                                                                                                               >  1:1183386/1‑150 (MQ=255)
                         cTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTg                                                                                                                               <  1:891506/150‑1 (MQ=255)
                          tGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACCCCGAGAAAGCCGATGGGGTGa                                                                                                                              >  2:1577093/1‑150 (MQ=255)
                          tGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGa                                                                                                                              >  1:647247/1‑150 (MQ=255)
                           ggcggcTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGAt                                                                                                                             >  1:670900/1‑150 (MQ=255)
                                 tGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCa                                                                                                                       <  2:588421/150‑1 (MQ=255)
                                          gcTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGt                                                                                                              <  2:1374776/150‑1 (MQ=255)
                                          gcTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGt                                                                                                              >  2:927260/1‑150 (MQ=255)
                                           cTGATTGGCAAACTGCCGACCCCGGAAGAGTACCAGACCTACGTGGCGCAGCTAGATAAAACAGCCGTAGATACTTCACGATATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGAAGTGGTGATTTTCCAGAATGCGGtt                                                                                                             >  1:284426/1‑150 (MQ=255)
                                               ttGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGACGGTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTaaa                                                                                                         <  1:97642/150‑1 (MQ=255)
                                                 ggCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTaaaaa                                                                                                       <  1:29524/150‑1 (MQ=255)
                                                 ggCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTaaaaa                                                                                                       <  1:61322/150‑1 (MQ=255)
                                                  gCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGGGGCGCAGGTAGATAAAACAGCCGGTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAg                                                                                                      <  2:826171/150‑1 (MQ=255)
                                                       cTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAgc                                                                                                 >  1:578372/1‑150 (MQ=255)
                                                         gCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAgcgc                                                                                               >  1:1249267/1‑150 (MQ=255)
                                                            gACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCTTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCAcg                                                                                            <  2:161087/150‑1 (MQ=255)
                                                              cGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGAGCAcgag                                                                                          >  1:13162/1‑148 (MQ=255)
                                                                 cGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGACAGCGCAAGCGCTg                                                                                       >  1:93628/1‑150 (MQ=255)
                                                                  ggAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTgc                                                                                      >  1:1052829/1‑150 (MQ=255)
                                                                  ggAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTgc                                                                                      <  1:1208136/150‑1 (MQ=255)
                                                                  ggAAGAGTACCAGACCAACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTATACAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTgc                                                                                      >  1:96980/1‑150 (MQ=255)
                                                                     agagTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCa                                                                                   >  2:1472789/1‑150 (MQ=255)
                                                                      gagTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCCGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCAt                                                                                  <  1:958181/150‑1 (MQ=255)
                                                                      gagTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCAt                                                                                  <  1:1561917/150‑1 (MQ=255)
                                                                      gagTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCAt                                                                                  <  1:913735/150‑1 (MQ=255)
                                                                          aCCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCATCCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAAc                                                                              <  1:17352/150‑1 (MQ=255)
                                                                          aCCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAAc                                                                              >  1:1326567/1‑150 (MQ=255)
                                                                          aCCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACATCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATTGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATATAc                                                                              >  2:927156/1‑150 (MQ=255)
                                                                                ccTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGAcaca                                                                        >  2:171090/1‑150 (MQ=255)
                                                                                    cGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGc                                                                    >  1:381058/1‑150 (MQ=255)
                                                                                       ggCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCg                                                                 <  1:75611/150‑1 (MQ=255)
                                                                                        gcgcAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCgg                                                                <  2:1195466/150‑1 (MQ=255)
                                                                                         cgcAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGt                                                               <  1:111072/150‑1 (MQ=255)
                                                                                          gcAGGTAGATAAAAGAGCCGTTGATAATTACCGTTATGTGAACTTCTACCAGGTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTTAAAATTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTg                                                              <  2:593325/150‑1 (MQ=255)
                                                                                             ggTAGATAAAACAGCCGTTTATACTTACCGTTATCTGAACTTCAAACAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAAt                                                           <  1:374589/150‑1 (MQ=255)
                                                                                              gTAGATAAAACAGCCGTTGATACTTACCGTTATCTGACCTTCAACCAGCATTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATg                                                          <  2:906735/150‑1 (MQ=255)
                                                                                              gTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATg                                                          <  1:14364/150‑1 (MQ=255)
                                                                                                  aTAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGAt                                                      >  1:1463818/1‑150 (MQ=255)
                                                                                                    aaaaCAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATTCCCGGTGAATGAGATTc                                                    >  2:876053/1‑150 (MQ=255)
                                                                                                        cAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACAGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCgg                                               >  1:223837/1‑150 (MQ=255)
                                                                                                         aGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCggg                                              >  2:1488110/1‑150 (MQ=255)
                                                                                                         aGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCggg                                              <  2:985910/150‑1 (MQ=255)
                                                                                                          gCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGc                                             >  2:521381/1‑150 (MQ=255)
                                                                                                            cGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGGCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCAt                                           <  1:1304589/150‑1 (MQ=255)
                                                                                                              tggATACTTACCGTTATCTGAACTTCACCCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAATGGCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGGAATCTCATTCACCGGGCAttg                                         <  1:302483/148‑2 (MQ=255)
                                                                                                               tGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCAtttt                                        >  2:1602064/1‑150 (MQ=255)
                                                                                                               tGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCAtttt                                        <  2:748707/150‑1 (MQ=255)
                                                                                                                taTACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCAttttt                                       <  2:1383372/149‑1 (MQ=255)
                                                                                                                  tACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCAGGGCAttttttt                                     >  1:1074934/1‑150 (MQ=255)
                                                                                                                      tACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTAttt                                 >  1:964381/1‑150 (MQ=255)
                                                                                                                       aCCGTTATCTGAACTTGAACCAGCTTGCTCAGTACCCCGCGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTATAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGATATTCCCGGGCATTTTTTTATTTc                                <  2:1194605/150‑1 (MQ=255)
                                                                                                                       aCCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTc                                >  1:115805/1‑150 (MQ=255)
                                                                                                                       aCCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTc                                <  2:134238/150‑1 (MQ=255)
                                                                                                                       aCCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTc                                <  2:466988/150‑1 (MQ=255)
                                                                                                                        ccGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGGAATCTCATTCACCGGGCAttgtgtcg                                    <  2:1374425/145‑7 (MQ=38)
                                                                                                                        ccGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGGAATCTCATTCACCGGGCAttgtgtcg                                    >  1:1374425/1‑139 (MQ=38)
                                                                                                                          gTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTaa                             <  2:1554272/150‑1 (MQ=255)
                                                                                                                          gTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTaa                             >  2:932825/1‑150 (MQ=255)
                                                                                                                           ttATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTaaa                            <  2:459673/150‑1 (MQ=255)
                                                                                                                             aTCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGAATTTTTTTATTTCTAAAcc                          >  2:1258367/1‑150 (MQ=255)
                                                                                                                                  aaCTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGc                     <  1:171090/150‑1 (MQ=255)
                                                                                                                                   aCTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGcc                    <  2:855937/150‑1 (MQ=255)
                                                                                                                                    cTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCg                   >  1:313670/1‑150 (MQ=255)
                                                                                                                                        aaCCAGCTTTCTCAGTACACCGAGAAAGGCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGATGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTcc               <  2:1183386/150‑1 (MQ=255)
                                                                                                                                        aaCCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTcc               <  2:1249267/150‑1 (MQ=255)
                                                                                                                                          ccAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTcctc             >  1:638461/1‑148 (MQ=255)
                                                                                                                                          ccAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGc             <  1:1472789/150‑1 (MQ=255)
                                                                                                                                          ccAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGAGTTAACGACACTATGAGGGTTGAATGAGGTTACCTGGCATTTTTGTAGTTCTAAAGCATCGCCGTTcggg             >  1:31775/1‑147 (MQ=255)
                                                                                                                                           cggCTTTCTCGGTGTACTGAGAAAGCTGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCt            >  2:713577/3‑150 (MQ=255)
                                                                                                                                            aGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTg           >  1:837182/1‑150 (MQ=255)
                                                                                                                                            aGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTg           <  2:1006791/150‑1 (MQ=255)
                                                                                                                                            aGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTg           >  1:486243/1‑150 (MQ=255)
                                                                                                                                              cTTTCTCGGTGTACTGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACAATCGCCGTTCCGCTGtt         >  2:310634/1‑150 (MQ=255)
                                                                                                                                                ttCTCTGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGTTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGtttt       <  2:1031348/150‑1 (MQ=255)
                                                                                                                                                ttCTCAGTACACCGAGAAAGCTGAAGGGGCGATTTCCCAGACTGCGGTGTAAAAAGTCAGCGCACGCGCTGCGCATAAACGAAACAATGCCCGGTGAATGAGATTCACGTGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGtttt       <  2:906137/150‑1 (MQ=255)
                                                                                                                                                ttCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGtttt       <  2:1380354/150‑1 (MQ=255)
                                                                                                                                                ttCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCAGTTTTTTATTTCTAAACCATCGCCGTTCCGCTGtttt       >  1:1458694/1‑150 (MQ=255)
                                                                                                                                                   tcAGTACACAGAGAAAGACGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGACATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTtct    <  1:758980/150‑1 (MQ=255)
                                                                                                                                                    cAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTtctc   >  2:652169/1‑150 (MQ=255)
                                                                                                                                                    cAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTtctc   >  2:879696/1‑150 (MQ=255)
                                                                                                                                                     aGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTTCTcc  >  2:1626145/1‑150 (MQ=255)
                                                                                                                                                     aGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTTCTcc  >  2:1664294/1‑150 (MQ=255)
                                                                                                                                                     aGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTAAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTTCTcc  <  2:670900/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GAATGTCTTCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTTCTCC  >  NC_000913/134017‑134315

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: