Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,387,553 T→A 7.3% R212R (CGA→CGT aaeB ← aromatic carboxylic acid efflux pump subunit AaeB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,387,5530TA 7.3% 293.5 / 8.1 94R212R (CGA→CGTaaeBaromatic carboxylic acid efflux pump subunit AaeB
Reads supporting (aligned to +/- strand):  ref base T (48/39);  new base A (2/5);  total (50/44)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

ATTCTGAATAAGATAAGTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAAGAGCAAATC  >  NC_000913/3387409‑3387699
                                                                                                                                                |                                                                                                                                                  
aTTCTGAATAAGATAAGTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCAc                                                                                                                                               >  1:188647/1‑150 (MQ=255)
aTTCTGAATAAGATAAGTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCAc                                                                                                                                               <  1:1523422/150‑1 (MQ=255)
  tCTGAATAAGATAAGTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCa                                                                                                                                             >  1:1438906/1‑150 (MQ=255)
   cTGAATAAGATAAGTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAg                                                                                                                                            <  2:1705559/150‑1 (MQ=255)
    tGAATAAGATAAGTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAgg                                                                                                                                           >  2:496454/1‑150 (MQ=255)
     gAATAAGATAAGTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGACGCACCacgg                                                                                                                                          <  2:1182774/150‑4 (MQ=255)
      aaTAAGATAAGTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGACGCACCacggc                                                                                                                                         <  1:1114922/150‑5 (MQ=255)
      aaTAAGATAAGTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGACGCACCacggc                                                                                                                                         <  1:289615/150‑5 (MQ=255)
        taagataagTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGc                                                                                                                                       <  1:1033781/150‑1 (MQ=255)
          agataagTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGccc                                                                                                                                     >  2:975154/1‑150 (MQ=255)
               agTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGACCTGGTGCGACGCACCACGgcgctacaagg                                                                                                                                <  1:1486023/150‑8 (MQ=255)
                gTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCt                                                                                                                               >  1:137008/1‑150 (MQ=255)
                gTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCt                                                                                                                               <  1:420748/150‑1 (MQ=255)
                gTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCt                                                                                                                               >  1:310975/1‑150 (MQ=255)
                    cGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTAt                                                                                                                           <  1:1291050/150‑1 (MQ=255)
                     gCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAgg                                                                                                                                           >  2:1244824/1‑133 (MQ=255)
                     gCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAgg                                                                                                                                           <  1:1244824/133‑1 (MQ=255)
                     gCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATc                                                                                                                          <  2:347596/150‑1 (MQ=255)
                     gCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATc                                                                                                                          <  2:1703116/150‑1 (MQ=255)
                     gCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATc                                                                                                                          <  1:27643/150‑1 (MQ=255)
                     gCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATc                                                                                                                          <  1:1662343/150‑1 (MQ=255)
                        ggATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGAc                                                                                                                       >  1:1325898/1‑150 (MQ=255)
                        ggATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGAc                                                                                                                       <  2:1325898/150‑1 (MQ=255)
                        ggATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGAc                                                                                                                       >  1:784222/1‑150 (MQ=255)
                             gggTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACtt                                                                                                                  >  1:642952/1‑150 (MQ=255)
                             gggTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACtt                                                                                                                  >  1:230019/1‑150 (MQ=255)
                              ggTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTc                                                                                                                 <  2:137008/150‑1 (MQ=255)
                                  aTCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCa                                                                                                             >  1:1949471/1‑150 (MQ=255)
                                    cagcgtcagcgATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATc                                                                                                           >  2:1511973/1‑150 (MQ=255)
                                    cagcgtcagcgATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTAACCAtt                                                                                                           <  1:624825/150‑2 (MQ=255)
                                     agcgtcagcgATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCg                                                                                                          >  1:1130267/1‑150 (MQ=255)
                                      gcgtcagcgATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGc                                                                                                         <  1:1702511/150‑1 (MQ=255)
                                        gtcagcgATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCa                                                                                                       >  2:778539/1‑150 (MQ=255)
                                          cagcgATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTATAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATg                                                                                                     <  1:570153/150‑1 (MQ=255)
                                            gcgATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCt                                                                                                   >  1:1615810/1‑150 (MQ=255)
                                            gcgATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCt                                                                                                   >  2:798479/1‑150 (MQ=255)
                                             cgATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCtt                                                                                                  <  1:1582946/150‑1 (MQ=255)
                                             cgATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCtt                                                                                                  <  1:496454/150‑1 (MQ=255)
                                              gATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTg                                                                                                 <  2:1209710/150‑1 (MQ=255)
                                                    gTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACag                                                                                           <  2:1628346/150‑1 (MQ=255)
                                                         gATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTCCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGTACGCCCCAGCCTTTATCGACAACTACACCAACGCCATGCTAGAAACAGAGTTg                                                                                      >  2:832207/1‑150 (MQ=255)
                                                         gATCGCTTTTAAACGTCCATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCATTGTACCGCCTTGTTAAGTCGCACCAGCTAGCCCCAGGCTTTATAGACAAATTAACCATCGCCATGCATTATACAGAGATg                                                                                      >  1:179547/1‑150 (MQ=255)
                                                          aTCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGc                                                                                     >  2:238131/1‑150 (MQ=255)
                                                                 ttAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAAtt                                                                              <  2:319228/150‑1 (MQ=255)
                                                                  tAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTg                                                                             >  1:1978094/1‑150 (MQ=255)
                                                                      cGTCGATTGGCCCGCGCCCAGCGGGAAGATTCGATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCATTATCGACAACATCACCATCGCCATGCTTGATACACAGTTGCATTAATTGatat                                                                         >  2:1023135/1‑150 (MQ=255)
                                                                        ccGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTg                                                                       <  1:876349/149‑1 (MQ=255)
                                                                          gATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCt                                                                                                   <  1:1857083/120‑1 (MQ=255)
                                                                          gATTGGCCCGCGCCCAGAGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCt                                                                                                   >  2:1857083/1‑120 (MQ=255)
                                                                              ggCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGAc                                                                 >  2:1181418/1‑150 (MQ=255)
                                                                               gCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGAcc                                                                >  2:936113/1‑150 (MQ=255)
                                                                                cccGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCa                                                               >  1:1649084/1‑150 (MQ=255)
                                                                                cccGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCa                                                               >  2:450089/1‑150 (MQ=255)
                                                                                cccGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCa                                                               >  2:1758846/1‑150 (MQ=255)
                                                                                      cccAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAAc                                                         <  2:188647/150‑1 (MQ=255)
                                                                                      cccAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAAc                                                         <  2:1615810/150‑1 (MQ=255)
                                                                                       ccAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACAATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACt                                                        >  2:311978/1‑150 (MQ=255)
                                                                                           cGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTcc                                                    >  1:696211/1‑150 (MQ=255)
                                                                                            gggAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCa                                                   >  2:1908304/1‑150 (MQ=255)
                                                                                            gggAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCa                                                   <  2:642952/150‑1 (MQ=255)
                                                                                            gggAAGATTCCATATTCAGGTTGCCGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCa                                                   >  2:1516444/1‑150 (MQ=255)
                                                                                             ggAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAg                                                  <  2:1949471/150‑1 (MQ=255)
                                                                                                 gATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGctct                                              >  2:1177175/1‑150 (MQ=255)
                                                                                                 gATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGctct                                              >  2:864418/1‑150 (MQ=255)
                                                                                                     ccATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGAt                                          >  2:1001199/1‑150 (MQ=255)
                                                                                                     ccATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTAGCCCCAGGCTTTATCGACAACTTCACCATAGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGAt                                          >  2:494577/1‑150 (MQ=255)
                                                                                                      cATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATc                                         >  1:1255292/1‑150 (MQ=255)
                                                                                                      cATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTACAGCTCTCGATc                                         >  1:201986/1‑150 (MQ=255)
                                                                                                            cAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCActtc                                   >  2:364720/1‑150 (MQ=255)
                                                                                                              ggTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCActtctt                                 >  1:1188117/1‑150 (MQ=255)
                                                                                                                 tgctgcGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGtt                              <  2:310975/150‑1 (MQ=255)
                                                                                                                  gctgcGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTa                                                                                 >  2:383869/1‑98 (MQ=255)
                                                                                                                  gctgcGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTa                                                                                 <  1:383869/98‑1 (MQ=255)
                                                                                                                   ctgcGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTg                            <  2:1438906/150‑1 (MQ=255)
                                                                                                                   ctgcGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTg                            >  2:675898/1‑150 (MQ=255)
                                                                                                                   ctgcGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGGCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTg                            <  2:132926/150‑1 (MQ=255)
                                                                                                                      cgcATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCTCCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTgatc                         <  2:102948/150‑1 (MQ=255)
                                                                                                                      cgcATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTgatc                         <  2:396472/150‑1 (MQ=255)
                                                                                                                      cgcATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTgatc                         <  2:696211/150‑1 (MQ=255)
                                                                                                                         aTGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTgatcgat                      >  2:792023/1‑150 (MQ=255)
                                                                                                                             cTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATcgcg                  <  2:1415546/150‑1 (MQ=255)
                                                                                                                             cTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATcgcg                  >  2:170915/1‑150 (MQ=255)
                                                                                                                              ttGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCgg                 <  1:778539/150‑1 (MQ=255)
                                                                                                                              ttGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCgg                 <  1:975154/150‑1 (MQ=255)
                                                                                                                                   gcgcCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGaaa            >  2:1359146/1‑150 (MQ=255)
                                                                                                                                   gcgcCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGaaa            >  1:1983559/1‑150 (MQ=255)
                                                                                                                                     gcCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAag          <  1:612628/150‑1 (MQ=255)
                                                                                                                                     gcCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAag          <  1:311978/150‑1 (MQ=255)
                                                                                                                                     gcCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAag          <  1:1359146/150‑1 (MQ=255)
                                                                                                                                      acctggtgCGACGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTgat                          <  2:1102980/130‑1 (MQ=255)
                                                                                                                                      acctggtgCGACGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTgat                          >  1:1102980/4‑133 (MQ=255)
                                                                                                                                          ggtgCGACGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAAGAGCaa     >  2:422977/1‑150 (MQ=255)
                                                                                                                                             gCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAAGAGCAAATc  >  1:967301/1‑150 (MQ=255)
                                                                                                                                             gCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAAGAGCAAATc  >  1:1170897/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
ATTCTGAATAAGATAAGTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAAGAGCAAATC  >  NC_000913/3387409‑3387699

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